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Abstract 

This  thesis  extends  the  AFIT  research  directed  towards 
replacing  a  standard  correlation  tracker  with  a  Kalnan 
filter  bank/enhanced  correlation  tracker  in  a  high  energy 
laser  weapon  system.  Airborne  targets  are  tracked  by  a 
Bayesian  multiple  model  adaptive  filtering  (HHAF)  algorithm, 
which  utilizes  an  array  of  infrared  sensing  detectors  as  the 
measurement  information  for  two-dimensional  position  data. 
TVo  different  target  dynamics  models  are  exercised:  a 
linear,  Gauss-Markov  accleration  model,  and  a  nonlinear, 
constant  turn-rate  model.  Performance  analyses  are 
accomplished  via  Monte  Carlo  simulation  techniques. 

Extending  the  adaptive  potential  of  the  tracking  algorithm 
is  of  primary  emphasis.  The  effects  of  bending  and 
vibration  of  a  large  space  structure  on  the  PLIR's  ability 
to  resolve  target  position  is  analyzed.  Also,  a  performance 
compar ison/simulation  time  tradeoff  is  conducted  with  the 
tracking  algorithm  operating  ot  both  30  Hz  and  50  Hz. 
Sensitivity  studies  of  adaptive  responsiveness  to  varying 
target  trajectories,  various  filter-assumed  correlation 
times,  range  to  pixel  size  relationships,  and  pixel  size  to 
filter  driving  white  noise  strength  relationships  are 
performed.  The  robustness  of  the  multiple  model  algorithm 
is  demonstrated  by  its  ability  to  adapt  to  scenarios  which 
it  had  not  been  previously  tuned. 


A  MULTIPLE  MODEL  ADAPTIVE  TRACKING 


ALGORITHM  AGAINST  AIRBORNE  TARGETS 

I .  Introduction 

In  view  of  the  Strategic  Defense  Initiative,  increased 
interest  has  been  placed  on  the  laser  as  a  prime  candidate 
for  a  potential  weapon  in  a  space-based  defensive  system. 
Because  of  a  laser’s  ability  to  transmit  energy  from  the 
weapon  to  the  target  at  the  speed  of  light,  the  need  to 
compute  a  lead  angle,  a  necessity  with  ballistic  projec¬ 
tiles,  is  eliminated.  However,  several  factors  affect  a 
laser's  effectiveness  on  a  given  target:  the  medium  in  which 
the  beam  travels,  the  power  or  intensity  of  the  beam,  the 
capability  to  acquire  the  target,  and  the  ability  to  track 
the  target  for  a  sufficient  time  to  neutralize  it. 

The  precision  pointing  and  accurate  tracking  capabili¬ 
ties  are  crucial  to  the  development  of  a  laser  weapon 
system.  It  would  not  suffice  to  "paint"  the  target  with 
laser  energy  nor  is  it  feasible  to  have  a  space-based  laser 
powerful  enough  that  would  neutralize  a  target 
instantaneously.  Limited  energy  is  available  In  the  laser 
beam,  thus  motivating  the  research  into  highly  accurate 
pointing  and  tracking  systems  for  space  applications. 


The  Air  Force  Weapons  Laboratory  (AFWL)  at  Klrtland 
AFB,  New  Mexico,  is  currently  developing  and  testing  a  high 
energy  laser  weapon  for  use  against  airborne  targets.  The 
current  tracking  method  employs  a  Forward  Looking  Infrared 
Sensor  (FLIR)  to  detect  the  target  passively  and  maintain  a 
low  probability  of  detection  by  the  target.  The  FLIR  detec¬ 
tor  plane  comprises  a  300  x  500  array  of  pixels,  or  picture 
elements,  where  each  pixel  can  focus  detected  energy  through 
an  angle  of  20  micro-radians  in  two  orthogonal  directions. 
The  tracking  algorithm  extracts  a  "tracking  window"  from  the 
larger  array  and  processes  the  detector  outputs  to  align  the 
tracker's  f teld-of -view  (FOV).  This  algorithm  detects  tar¬ 
get  offsets  from  the  center  of  the  FOV. 

Currently,  AFWL  employs  a  correlation  tracker  to  pro¬ 
cess  pointing  and  tracking  sensor  information.  This  tracker 
compares  the  previous  sample  of  FLIR  information  to  the 
present  information.  Relative  position  offsets  from  one 
data  frame  to  the  next  are  cross  correlated  and  this  infor- 
oiation  update  commands  the  control  system  to  center  the 
target  in  the  f le Id-of-view.  This  enables  the  laser  to 
maintain  lock  on  a  given  target  since  the  FLIR  and  laser 
share  the  same  optics. 

Since  no  a  priori  knowledge  of  target  characteristics 
is  required,  the  correlation  tracking  algorithm  can  be 
applied  to  a  wide  variety  of  targets;  however,  it  exhibits 
several  limitations.  First,  the  correlation  algorithm  does 


not  take  advantage  of  target  dynamic  characteristics,  which 
have  been  thoroughly  documented.  In  addition,  it  is  highly 
susceptible  to  noise  and  is  not  capable  of  distinguishing 
between  apparent  target  motion  due  to  signal  corruption,  as 
a  result  of  atmospheric  distortion,  and  true  target  motion. 
This  "jitter**  effect  (18)  could  cause  a  translation  in  the 
FLIR  image  plane  and  thus  result  in  the  laser  beam  pointing 
in  a  direction  other  than  towards  the  target.  Mirror  vibra¬ 
tion  effects  and  system  component  bending  effects  also  con¬ 
tribute  to  the  image  plane  translation  phenomena.  Another 
limitation  is  the  time  difference  betvfeen  computing  the 
image  correlation  and  the  actual  time  for  the  gimbal  system 
to  engage  for  pointing.  A  tracking  algorithm  that  antici¬ 
pates  future  target  position  based  upon  past  and  present 
information  could  compensate  for  these  limitations. 

Since,  under  the  appropriate  conditions,  it  is  an 
optimal  linear  estimator,  the  Kalman  filter  (9)  is  ideally 
suited  for  overcoming  these  limitations.  By  modeling 
different  types  of  targets,  various  parameters  such  as  size, 
shape,  and  acceleration  characteristics  of  each  target  type 
would  either  be  known  or  could  be  estimated.  If  statistical 
characteristics  of  measurement  errors  and  atmospheric  jitter 
are  included,  the  filter  can  predict  target  position  more 
accurately.  This  prediction,  or  optimal  estimate,  reduces 
tracking  errors  due  to  time  delays  or  pointing  system  lags 
by  allowing  the  PLIR  to  anticipate  target  motion. 


1.2 


Previous  AFIT  Research 

Over  the  past  eight  years,  the  Air  Force  Institute  of 
Technology  students  and  staff  have  generated  numerous  papers 
and  theses  investigating  the  feasibility  and  performance 
benefits  of  both  enhanced  correlation/linear  Kalman  filter 
tracking  and  extended  Kalman  filtering  tracking  algorithms, 
all  incorporating  online  adaptations,  with  the  high  energy 
laser  pointing  and  tracking  system.  Both  Netzer  (14)  and 
Tobin  (19)  had  suimnarized  their  predecessor's  endeavors  and 
a  modification  of  their  observations  follows. 

Initially,  the  study  by  Herder  (12)  demonstrated  that 
the  extended  Kalman  filter  (EKF)  algorithm  outperformed  the 
standard  correlation  tracker  algorithm  in  the  ability  to 
track  long  range  targets  modeled  as  infrared  radiation  point 
sources.  The  FLIR  plane  image  of  the  target  was  assumed  to 
have  a  bivariate  Gaussian  distribution  and  was  modeled  with 
equal-intensity,  circular  contours.  The  four-state  filter 
was  predicated  on  a  benign  target  dynamics  model  with  a 
first  order,  zero-mean  Gauss-Markov  (GM)  position  model. 

The  intensity  distribution  due  to  atmospheric  disturbances 
was  initially  approximated  by  a  third  order  shaping  filter, 
and  then  this  was  replaced  by  a  first  order  shaping  filter 
driven  by  white,  Gaussian  noise  through  reduced  order 
modeling.  Filter  measurement  noise  was  considered  to  be 
uncorrelated  in  both  time  and  space.  This  tracking 
algorithm  enhanced  tracking  performance  an  order  of 


magnitude  better  than  observed  with  the  correlation  tracker 
in  benign  scenarios. 

Harnly  and  Jensen  (4)  incorporated  velocity  and  accel¬ 
eration  estimates  to  accommodate  for  the  tracking  of  more 
maneuverable  targets.  They  extended  the  FLIR  plane  image 
constant-intensity  contours  to  represent  an  elliptical  con¬ 
figuration  oriented  at  specified  angles  such  that  the  major 
axis  was  aligned  with  the  estimated  velocity  vector  (versus 
circular  contours),  in  addition  to  incorporating  a  spatially 
correlated  Gaussian  noise  model  to  represent  FLIR  sensor  and 
background  noises.  They  also  Incorporated  a  maneuver  detec¬ 
tion  algorithm  and  a  means  of  appropriately  responding  to 
detected  maneuvers  by  gain  changing  and  modification  of  the 
dynamics  models.  Finally,  they  modified  the  algorithm  to 
estimate  the  target's  true  size  and  shape  adaptively. 

In  previous  research  efforts,  the  extended  Kalman  fil¬ 
ters  were  given  target  intensity  function  information, 
specifically  that  bivariate  Gaussian  shapes  adequately 
depicted  the  targets,  perhaps  requiring  some  estimation  of 
parameters  associated  with  that  function.  Research  by 
Singletery  (17)  and  Rogers  (16)  implemented  algorithms  which 
made  no  a  priori  assumption  about  target  shape  and  tested 
the  algorithms  against  targets  with  multiple  hot-spots  and 
several  dynamic  orientations.  Rogers  developed  an  alterna¬ 
tive  tracker  which  used  the  target  shape  function  as  a 
template  for  an  enhanced  correlator.  The  correlator 


produced  oieasurement  offsets  In  the  two  orthogonal  direc¬ 
tions  from  the  center  of  the  FOV  and  fed  these  measurements 
as  inputs  to  a  linear  Kalman  filler.  A  linear  Kalman  filter 
rather  than  a  non-linear  extended  Kalman  filter  could  be 
utilized  since  the  measuremer ts,  or  offset  distances,  were 
linear  functions  of  the  filter  states.  When  analogous 
performance  is  achievable,  the  linear  Kalman  filter  is 
preferred  over  the  extended  Kalman  filter  due  to  the  lower 
level  of  computational  loading. 

Kozemchak  (7)  continued  the  research  by  implementing 
the  digital  signal  processing  techniques  of  Rogers'  (16) 
shape  Identification  algorithm  for  the  case  of  highly 
dynamic  targets.  He  substantiated  that  filters  based  on 
both  the  Gauss-Harkov  target  acceleration  and  constant  turn- 
rate  target  dynamics  exhibited  good  tracking  performance  in 
the  presence  of  dynamic  images;  these  dynamic  images 
corresponded  to  realistic  projections  of  multiple  hot-spot 
targets  onto  the  FLIR  image  plane.  Tracking  performance  was 
also  favorable  when  experiencing  drastic  image  changes 
inherent  in  constant  g  and  constant  roll-rate  maneuvers. 

Further  research  by  Millner  (13)  included  a  data 
processing  algorithm  to  generate  an  estimated  intensity 
function  relationship  to  be  utilized  as  a  template  in  the 
correlation  algorithm,  which  is  the  same  as  Rogers'  (16) 
alternative  filter  idea,  although  with  a  modified  dynamics 
model  to  address  the  highly  maneuvering  target  scenario. 

This  method  showed  that  nearly  identical  tracking 


performance  was  achieved  In  both  the  single  and  multiple 
hot-spot  scenarios.  Ho%raver,  inadequate  target  tracking 
performance  occurred  for  targets  exhibiting  maneuvers  in 
excess  of  five  g's. 

In  an  attempt  to  thwart  this  limitation,  Flynn  (3) 
investigated  the  possibility  of  multiple  model  adaptive 
filtering  (MHAF)  techniques  for  filter  implementation. 

Suizu  (18)  pursued  this  effort  and  successfully  implemented 
the  HHAF  in  addition  to  demonstrating  that  it  performed  well 
against  a  wide  dynamic  range  of  targets.  The  model  con¬ 
tained  a  bank  of  two  filters,  one  tuned  for  highly 
maneuvering  targets  with  a  POV  of  24  X  24  pixels  and  the 
other  tuned  for  benign  targets  with  a  FOV  of  8  X  8  pixels. 
Based  upon  probabilistic  weighting,  the  filter  adaptively 
changed  the  assumed  target  dynamics  and  corresponding  POV 
array,  which  resulted  in  increasing  filter  performance  to 
allow  the  tracking  of  targets  varying  from  benign  straight- 
line  trajectories  to  pulling  20  g's  at  20  kilometers.  Both 
the  linear  Kalman  filter/correlation  algorithm  and  the  BKP 
were  tested,  and  both  algorithms  demonstrated  the  capability 
to  track  targets  with  trajectories  that  maneuvered  signifi¬ 
cantly  throughout  the  simulation  scenario. 

Loving  (8)  continued  the  rec:earch  by  adding  a  filter  to 
the  HMAF  bank  based  on  an  intermediate  level  of  target 
dynamics,  to  aid  in  tracking  highly  maneuvering  targets.  She 
developed  a  Maximum  a  Posteriori  (MAP)  algorithm  and  com- 


pared  it  to  the  Bayesian  estimator.  The  MAP  algorithm  took 
advantage  of  the  existing  MMAF  structure  and  generated  state 
estimates  from  the  one  elemental  filter  with  the  highest 
probability  of  validity,  versus  probabilistically  weighting 
all  elemental  filter  estimates.  The  addition  of  the  third 
elemental  filter  to  the  HHAF  bank  significantly  enhanced  the 
tracking  ability  of  highly  maneuvering  targets.  Both  the 
MAP  and  Bayesian  technigues  displayed  favorable  tracking 
performance  against  several  target  maneuvers. 

Pollow-on  research  by  Netzer  (14)  expanded  the  results 
produced  by  Loving  (8)  and  employed  the  linear  filter/corre¬ 
lator  algorithm  developed  by  Rogers  (16)  but  modified  for 
harsher  target  dynamics.  He  investigated  the  steady  state 
bias  errors  as  a  result  of  the  algorithm  tracking  a  high-g, 
constant  turn  rate  maneuvering  target.  It  was  determined 
that  the  predominant  cause  of  the  bias  was  due  to  mis¬ 
matching  of  the  FLIR  azimuth-channel  dynamics  model  and  FOV 
size  when  a  true  target  maneuver  in  the  elevation  direction 
forced  the  MMAF  to  select  a  harsh  dynamics  model  and  wide 
f ield-of-view  in  order  to  maintain  lock  on  the  target.  This 
motivated  the  investigation  of  a  MMAF  based  on  separate 
elemental  filters  that  are  tuned  specifically  for  harsh 
maneuvers  in  either  the  x-  or  y-channels,  which  would  enable 
the  tracker  to  distinguish  between  x-direction  and  y- 
direction  maneuvers.  This  technique  would  provide  the  capa¬ 
bility  for  the  tracker  to  expand  the  FOV  in  the  critical 


direction  and  maintain  lock  on  a  maneuvering  target  while 
retaining  maximum  resolution  In  the  non-crltlcal  direction. 

Most  recently,  Tobin  (19)  Implemented  a  constant-turn- 
rate  dynamics  model  since  the  zero-mean,  first  order  Gauss- 
Markov  acceleration  processes  did  not  adequately  describe 
target  dynamics  In  some  cases.  Although  he  demonstrated 
that  the  steady-state  standard  deviation  errors  were 
typically  sotaller  in  the  CTR  model,  the  GM  HHAF  outperformed 
the  CTR  MHAF  with  respect  to  the  transient  characteristics 
of  target  maneuvers.  Tobin  also  included  rectangular  POV 
elemental  filters  tuned  specifically  for  target  maneuvers  in 
both  the  X-  and  y-directions  and  ascertained  that  the 
tracker  maintained  lock  on  the  target  at  the  onset  of 
y-direction  maneuvers  while  retaining  maximum  resolution  in 
the  more  benign  x-direction. 

1 . 3  Objectives 

During  the  past  eight  years,  substantial  development 
has  been  accomplished  on  a  tracker  capable  of  handling 
multiple  hot-spot  targets  in  which  digital  signal  processing 
techniques  have  been  employed  on  FLIR  sensor  information  to 
Identify  the  underlying  target  shape.  Host  recently, 
adaptation  of  the  field  of  view  to  maintain  track  on  harshly 
maneuvering  targets  at  close  ranges  has  been  investigated 
through  multiple  model  adaptive  filtering  techniques.  Both 
the  correlator  tracker  and  the  multiple  model  adaptive 
filter  tracker  have  been  assessed  with  respect  to  several 


tracking  scenarios.  This  research  Is  motivated  by  the  fact 
that  a  further  assessment  of  the  capabilities  of  these  two 
designs  warrants  continued  investigation  In  an  attempt  to 
extend  their  applicability.  Therefore,  the  Intent  Is  to 
continue  the  research  effort  In  developing  a  viable  tracking 
algorithm  by  accomplishing  the  following  objectives. 

1.3.1  Bending  and  Vibration  Phenomena .  Previous 
research  had  considered  the  tracker  to  be  ground-based  and 
the  effects  of  structural  bending  and  vibration  were  felt  to 
be  negligible.  This  research  shall  consider  the  application 
of  the  tracker  to  a  space  structure  and  model  the  effects  of 
bending  and  vibration.  It  shall  analyze  the  resulting 
tracker  performance  when  the  filter  does  not  account  for 
these  effects  through  additional  states  to  be  estimated. 
Implementation  of  bendlng/vibratlon  phenomena  will  only  be 
incorporated  explicitly  into  the  truth  model,  while  filter 
retuning  may  be  performed  to  enhance  filter  performance 
without  increasing  its  state  dimension. 

1.3.2  Implementation  at  50  Hz.  The  current  tracking 
simulation  Is  Implemented  digitally  at  a  30  Hz  rate  and  the 
potential  utilization  of  a  SO  Hz  tracker  Implementation  Is 
presently  being  studied  by  certain  industry  sectors  for 
space  application  (5).  Therefore,  this  effort  will  pursue 
the  Implementation  of  the  current  simulation  software  at  a 
50  Hz  rate  and  analyze  the  potential  benefits. 

1.3.3  Scenario  Sensitivity  Parameters .  This  research 
shall  perform  a  sensitivity  analysis  to  assess  the 


applicability  o£  this  tracking  algorithm  to  a  possible  range 
of  application  scenarios.  This  will  encompass  developing 
the  relationships  to  vary  pixel  size,  target  range,  target 
type,  and  noise  character Istics .  Previous  research  by  both 
Netzer  (14)  and  Tobin  (19)  had  addressed  the  sensitivity  of 
the  filters  to  range;  however,  the  Intent  of  this  objective 
Is  to  develop  a  relationship  with  the  parameters  of  pixel 
size,  range,  noise  characteristics,  and  target  type  and 
evaluate  tracker  limits  as  a  function  of  these  variables. 

1.3.4  Target  Trajectory  Sensitivity.  A  comparison  of 
the  response  of  the  filters  to  target  trajectories  other 
than  those  for  which  the  filters  are  tuned  is  conducted. 

The  five  elemental  filters  in  Tobin's  multiple  model  filter 
algorithm  (19)  had  experienced  favorable  results  against  a 
straight  and  level  trajectory,  a  constant  2-g,  10-g  or  20-g 
pull-up  maneuver  with  the  maneuver  Initiating  in  an  inertial 
X“/y“  plane  parallel  to  the  FLIP  plane,  a  similar  constant-g 
maneuver  but  ending  in  a  straight  trajectory,  and  a  constant 
turn-rate  maneuver  Initiating  In  an  inertial  x-,z-  plane 
perpendicular  to  the  FLIP  plane.  This  effort  will  consider 
the  implementation  of  a  constant  turn-rate  maneuver  similar 
to  the  constant  pull-up  maneuver  initiating  in  the  x-,y- 
plane  and  rotate  the  trajectory  by  45*^.  This  maneuver  is 


expected  to  show  the  effects  of  tracking  a  target  whose 
initial  trajectory  is  as  dynamic  in  the  x-dlrectlon  as  it  is 
in  the  y-direction,  ratther  than  singling  out  one  of  the 


directions  for  initial  dynamic  behavior.  Another  trajectory 
modification  that  will  be  considered  is  to  represent  the 
target's  acceleration  at  the  onset  of  the  turning  maneuver 
as  an  increasing  function  of  time  rather  than  as  a 
physically  unrealistic  step  input  (as  had  been  done  in  the 
past).  The  target  acceleration  rate  will  be  modeled  as 
a  ramp  function  with  the  maximum  attainable  acceleration 
occurring  after  five  sample  periods. 

1.3.5  Rotating  Rectangular  FOV.  An  investigation  into 
the  feasibility  of  implementing  a  rotating  rectangular  field 
of  view  into  the  multiple  model  adaptive  filter  structure  is 
performed.  This  field  of  view  is  considered  to  rotate  in 
such  a  way  as  to  align  the  "elongated"  side  with  the  current 
estimate  of  the  target's  acceleration.  The  benefits  of 
maintaining  lock  on  a  target  that  can  accelerate  in 
arbitrary  directions  forms  the  basis  of  this  evaluation. 

1 . 4  Overview 

This  chapter  has  presented  a  synopsis  of  the  effort 
performed  to  date  in  developing  a  viable  tracking  algorithm 
as  well  as  described  the  intentions  of  the  areas  requiring 
continued  study.  Chapter  II  introduces  the  filter  theory 
concepts  that  will  be  needed  to  understand  the  algorithm's 
theoretical  foundation.  Chapter  III  provides  the  develop¬ 
ment  of  the  various  truth  models  and  Chapter  IV  discusses 
the  aspects  of  the  two  possible  filter-assumed  dynamics 
models,  a  Gauss-Harkov  acceleration  model  and  a  constant 


turn-rate  dynamics  model.  In  addition  to  the  measurement 
model  employed.  Chapter  V  discusses  the  tracking  algorithms 
along  with  discussing  the  tools  used  for  gathering 
statistics  and  describing  £ield-o£-view  processing 
characteristics.  Chapter  VI  presents  the  results  of  the 
various  analyses  conducted,  and  Chapter  VII  presents  the 
summary  of  results  and  suggests  areas  for  further  research. 
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II .  Filter  Theory 


The  Intent  of  this  chapter  is  to  highlight  the 
nathenatical  forns  of  the  Extended  Kaliean  Filter  (EKF)  and 
the  Multiple  Model  Adaptive  Filter  (HMAP) .  A  basic 
understanding  and  knowledge  of  linear  Kalman  filtering 
techniques  Is  assumed  (9,10). 

2.1  Extended  Kalman  Filter 

The  extended  Kalman  filter  (EKF),  for  which  a  complete 
derivation  appears  In  (10)  and  is  summarized  in  (19),  is  an 
algorithm  that  allows  the  states  of  a  nonlinear  stochastic 
system  model  to  be  estimated.  Unlike  the  linear  Kalman 
filter  (LKP),  the  extended  Kalman  filter  requires  the  Taylor 
seles  expansion  of  nonlinear  system  dynamics  and  measurement 
equations  at  the  desired  sample  intervals.  These  Taylor 
series  expansions  neglect  the  effect  of  second  and  higher 
order  terms,  which  results  in  the  EKF  not  complying  with  the 
optimality  criteria  inherent  in  the  linear  Kalsian  filter. 

The  system  state  relationship  for  the  EKF  takes  the  form  of 
the  following  nonlinear  stochastic  differential  equation: 

x(t)  =  f lx(t),u(t),t]  +  G(t)  K(t)  (2-1) 

where : 

2e(t)  a  n-dimensi onal  state  vector 
ll(t)  «  r-dlmenslonal  control  input  vector 
t  time 

f ( • ]  ■  n-dimensional  nonlinear  system  plant  dynamics 

vector  function 


£(t)  ■  n  X  s  noise  distribution  matrix 

j£(t)  s  s-vector  zero-mean,  white  Gaussian  noise 

process  of  strength  ft(t);  independent  of  2c(to) 

As  compared  to  the  linear  Kalstan  filter,  l(x(t),j2.(t),t] 
replaces  I£(t)x(t)  ♦  6(t)ji(t)},  but  similar  to  the  LKP  is 
the  assumption  that  the  initial  condition  is  assumed  to  be  a 

A 

(Gaussian)  random  n-vector  with  Man  jcq  and  covariance 
The  discrete-time  measurements  can  be  represented  in  the 
form  of  the  following  nonlinear  vector  function: 

i(ti)  »  h(x(ti),ti)  ♦  v(ti)  (2-2) 

where : 

2(t|)  =  m-dimensional  measurement  vector 

h I • J  -  m  X  n  nonlinear  vector  function  relating 
the  measurements  to  the  states 

v(tj)  a  m-vector  process  of  discrete-time  zero-mean, 
white  Gaussian  noise  with  covariance  R(t£); 
independent  of  both  and  w(t) 

The  extended  Kalman  filter  measurement  update  incorporates 
the  measurements  via  the  following  relationships: 

!<(ti)  =  E(ti”)li*'(ti){H(ti)E(ti")!i’’(ti)  +  R(ti))"^  (2-3) 

x(ti*)  *  X(ti‘)  ♦  lL(ti){ii  -  l»(x(ti"),ti  ))  (2-4) 

E(ti'^)  a  E(ti")  -  IS.(ti)H(ti)E(ti')  (2-5) 

where : 

|((tj^)  a  Kalman  filter  gain  at  time  t^ 

A 

2i.(t£)  ■  estimate  of  2.  tiM  tj^ 

E(tj^)  a  n  x  n  state  error  covariance  matrix 


(t£  )  ■  instant  Inaiediately  prior  to  incorporating 
measureswnts  at  tiM  t^ 

(t^^)  =  instant  iaaMdiately  after  incorporating 
sieasureBients  at  tiaie  t^ 

The  ii(ti^)  Mtrix  is  defined  as  an  at  x  n  matrix  of  partial 
derivatives  of  ^  with  respect  to  the  first  arguaent, 
evaluated  along  the  noainal  trajectory  (10:41): 

-  dti(x,  tj) 

tt(ti)  -  tt(x(ti~),  ti)  «  -  -  _  (2-6) 

d)c  S.*x(ti  ) 

Utilizing  the  results  of  the  aeasureaent  update  equa¬ 
tions  for  X  and  as  defined  in  Equations  (2-4)  and  (2-5)  as 
the  initial  conditions  for  the  propagation  cycle  froa  t^  to 

tf^Xf  the  following  relationships  apply: 

« 

A  A 

}t(t/ti)  «  i(s.(t/ti),a(t),t)  (2-7) 

E(t/ti)  -  E(t)E(t/ti)  t  E(t/ti)E’’(t) 

+  G(t)g(t)£'’^(t)  (2-8) 

where  : 

(t/t^)  •=  estiaate  at  tiae  t  given  aeasureaents  through 
time  ti 

with  the  initial  conditions: 

A  A 

X(ti/ti)  =  Xiti"*^)  (2-9) 

P(ti/ti)  -  (2-10) 

The  E<t)  aatrix  is  defined  as  the  n  x  n  matrix  of  partial 
derivatives  of  f.  with  respect  to  its  first  arguaent, 
evaluated  along  the  noainal  trajectory  (10:41): 


.  ,% 

v'* 


i: 


E(t)  -  E(2L(t/ti  ),ti  1  - 


dx 


3t-X(t/ti  ) 


(2-11) 


Note  that  one  difference  between  the  linear  Kalnan  filter 
and  the  extended  Kalun  filter  la  that  the  equations  for 
propagating  and  updating  the  eatiution  error  covariance 
Mtrix  are  coupled  to  the  state  estiaiate  relations.  This 
relationship  precludes  the  precoskputat ion  of  the  covariance 
and  gain  okatrices  until  the  state  estimates  and  sieasureMnt 
values  become  knovm  at  specified  sample  times. 


2 . 2  Multiple  Model  Adaptive  Filter 

In  any  given  real  world  application,  one  can  predict 
only  to  some  degree  of  accuracy  all  of  the  possible  scena¬ 
rios  of  a  deterministic  target  model,  thus  allowing  the  idea 
of  uncertainties  in  the  model.  The  Kalman  filter's  ability 
to  achieve  a  high  level  of  performance  against  a  wide 
variety  of  situations  would  require  the  switching  of  the 
uncertain  parasieters  of  the  dynamics  model  to  the  dynamics 
of  the  target.  However,  some  of  these  paraswters  are 
changing  in  time  and  the  system  designer  does  not  always 
have  a  priori  Icnowledge  of  the  values  of  parameters  that 
provide  optimal  perforswince.  One  swans  of  overcoming  this 
limitation  is  to  incorporate  Multiple  Model  Adaptive  Filter 
(MMAF)  techniques  as  developed  by  (8,14,18,19,10). 

The  MMAF  structure  can  best  be  developed  (10:129,136) 
by  considering  a  first-order,  linear,  stochastic 
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dlffecentlal  equation  for  a  given  tl»e-invar lant  system 
model  of  the  form: 

x(t)  »  E(a)x(t)  ♦  t(a)a(t)  +  G(a)w(t)  (2-12) 

and  noise  corrupted,  discrete  tiaie  measureotents  as: 

l(ti)  •  H(a)x(ti)  +  y(ti)  (2-13) 

where : 

x(t)  «  n-diamnsional  system  state  vector 

y.(t)  «  r-diamnsional  deterministic  control  vector 

w(t)  «  s-dimenslonal  discrete-t ism  white  Gaussian, 
zero-otean  noise  vector  process  of  strength 

Q(a) 

2.(t|)  m-distensional  measurement  vector 

vitj)  *  m-dimensional  discrete-tism  white 

Gaussian,  zero-sman  noise  vector  process 
of  covariance  g.(i.) 

a  =  uncertain  and/or  dynamic  parameter  vector 

£(a)  =  n  X  n  system  plant  matrix 

B(a)  «  n  x  r  input  distribution  smtrix 

S(a)  -  n  X  s  noise  distribution  statrix 

H(a)  >  m  x  n  vector  relating  the  measuresmnts  to 
the  states 

To  represent  a  target  with  K  significantly  different  sets  of 
dynamics  characteristics,  it  is  necessary  to  discretize  a 
into  a  set  of  K  finite  vector  values,  a^r  •••/  SK'  The 

HHAP  consists  of  a  bank  of  K  independent  Kalsian  filters 
which  are  processed  In  parallel,  in  which  each  filter  Is 
tuned  for  a  specific  target  dynamics  characteristic 


I 

n| 

iS 


i  •: 


determined  by  the  appropriate  41^,  for  k-  1,2,  K.  At 

discrete  sample  tlsws,  the  residuals  o£  each  elesMntal 
filter  are  used  to  calculate  the  conditional  probability 
that  4  assuMS  the  value  of  4|(  associated  with  that 
particular  elemental  filter,  conditioned  on  the  observed 
measurement  history.  This  conditional  probability,  called 
the  hypothesis  conditional  probability  and  denoted  P|((t^), 
Identifies  the  elemental  filter  with  the  greatest 
probability  of  best  performance  at  a  given  time. 

The  state  estloiate  of  the  HHAF,  t » .  Is  the 

probabilistically  %/eighted  average  of  the  elesiental  bank  of 


filters  (19:14): 


Jtj„*af(ti^>  =  ^  Pk<ti)  JCKCti"*^) 

k«l 


(2-14) 


This  form  of  filter  structure,  as  depicted  in  Figure  2.1,  is 
referred  to  as  the  Bayesian  HMAF  (10:132),  where  the  state 
estimate  consists  of  weighting  all  k  filters  using  the 
hypothesis  conditional  probability  function: 


Pli(ti)  = 


f  (S-i  l^/li-l^'Pk^^i-l) 

4<ti) l4/Z(ti-i) _ 

K  f  <11  laj/Zi-i)* Pj(ti-i) 

E  l(ti)  l4/Z.(ti_i) 

j-1 


(2-15) 


where : 

i  (aiiak^Zi-x) 

4(ti)|a,2(ti_i) 


exp  {  } 


, _  .m/2 , .  , .  .  ,  1/2 

(2k)  IHk^^l ) I 


(2-16) 


Figure  2.1.  Bayesian  Multiple  Model  Filtering  Algorithm 


As  is  evident  from  Equation  (2-15),  the  nu«erator  is 
the  filter's  product  of  the  previous  hypothesis  probabi¬ 

lity  and  the  conditional  probability  density  of  the  current 
MasuresMnt  2.(^1)  provided  that  the  k^^  filter's  assusied 
paraneter  value  Sk  previous  Masure»ent  history 

2.(ti.]^)  are  know.  The  denoainator  is  the  sum  of  slallar 
products  for  all  K  filters  in  the  bank;  this  implies  that 
the  filter  out  of  the  K  filters  that  produces  the  smallest 
squared  residual  relative  to  the  filter-computed  residual 
covariance  best  aiatches  the  real  world  situation  (10:17). 

The  filter  structure  of  Figure  2.1  conveys  that  the  K 
filters  in  the  bank  process  their  estimates  and  residuals  in 
parallel.  It  is  isiportant  that  the  residuals  from  the 
"best-statched"  filter  be  distinguishable  from  those  of  the 
mismatched  filters.  Failure  to  obtain  such  distinction 
can  cause  inappropriate  assignment  of  large  probabilities 
to  filters  based  on  incorrect  models,  leading  to  poor 
perforMnce.  To  preclude  such  possible  perforstance  degrada¬ 
tion,  each  filter  in  the  bank  should  be  tuned  for  a  unique 
target  trajectory  to  statch  its  Internal  dynamics  model, 
since  the  filter  which  represents  the  closest  "satch*  to  the 
true  target  dynamics  will  produce  much  smaller  residuals 
relative  to  the  f ilter-cosiputed  covariance  than  the 
misstatched  filters.  In  addition,  to  prevent  the  possibility 
of  masking  the  distinction  between  the  residuals  of  the 
filters  based  on  different  target  dynamics  models,  the  addi- 


tlon  of  excessive  asK>unts  of  pseudonoise  to  compensate  for 
Model  Inadequacies  should  be  avoided  ( 10 : 133; 14 : 18 ) . 

The  computation  of  the  MHAP  conditional  covariance 
matrix  takes  the  following  form  (19:17): 

K  ^ 

E«Maaf<^l*)  •  E  Pk(ti)IEk(t/)  +  ^  (ti*  )  1  (2-20) 

k*l 

where : 

Yk(ti^)  -  -  Xmmaf<t/) 

P)((t^)  =  k-th  filter's  conditional  hypothesis 
probabi 1 Ity 

Pj{(ti^)  =  k-th  filter's  state  error  covariance 

smtrlx  after  Incorporating  the  measuresient 

The  values  of  Eosaaf^^l^)  cannot  be  computed  a  priori  as  can 
be  accomplished  for  the  the  case  of  a  linear  Kalman  filter 
(as,  for  Instance,  each  elemental  filter)  because  the  values 
of  both  P)((t^)  and  x.^p^^( t ^  )  depend  upon  the  measurements 
taken  through  tlsm  t|^. 

Finally,  to  prevent  a  mismatched  filter's  Pt((ti)  value 
from  converging  to  zero,  an  artificial  lower  bound  Is 
Imposed  on  each  of  the  filter's  hypothesis  conditional 
probabilities  ( 10 : 135; 14 : 18  ) .  Without  the  lower  bound,  once 
a  conditional  probability  reaches  zero,  or  essentially 
reaches  zero.  It  resalns  zero  (or  essentially  zero)  for  all 
time,  since  it  is  a  function  of  the  conditional  probability 
at  the  previous  time,  as  seen  In  equation  (2-lS).  This 
effectively  Impedes  the  contributions  of  that  filter's 
estlskates  for  all  future  tlsms  In  the  MHAP  structure  and 


can  degrade  the  ability  o£  that  filter  to  respond  to  future 
changes  In  the  true  paraneter  values.  The  loss  of  a  filter 
due  to  the  oversight  of  not  placing  a  lower  bound  on  Its 
probability  value  could  degrade  the  HHAP  performance  if  the 
target  dynamics  model  at  a  later  time  matched  the  model  for 
which  the  probability  was  allotted  to  reach  zero,  thus 
preventing  It  being  weighted  in  the  filter  structure.  The 
establishment  of  a  lower  bound  of  P)c(ti)Bin  ‘  0.001  (8:19) 
in  past  efforts  is  continued  in  this  study. 


2. 4  Suiwaarv 


This  chapter  has  presented  discussions  of  both  the 
extended  Kalshan  filter  and  the  concepts  underlying  multiple 
model  adaptive  filtering  techniques.  The  Intent  of  this 
discussion  is  to  provide  a  basic  understanding  of  the  theory 
that  motivates  both  the  tracking  algorithm  and  the  filter 
BK)del  in  the  ensuing  chapters. 


23 


'/.w  III.  Hath  Mffdsl  p^velopwent 

3.1  Introduction 

The  truth  Model  o£  any  real  world  filtering  and/or 
control  application  ia  the  Most  coMplete  MathOMatical 
representation  available  to  the  systeM  designer.  This 
truth  Model  becoMes  the  standard  with  which  to  evaluate  the 
overall  perforMance  of  the  filter,  which  is  discussed  in  the 
next  chapter.  This  chapter  describes  the  characteristics  of 
the  truth  Model  as  utilized  in  previous  ground-based  tracker 
studies  and  Most  recently  in  (19).  Ho%[/ever,  the  Modeling  of 
hard%Mire  bending  and  vibration  effects  is  Included  to  repre¬ 
sent  additional  detector  perturbations  inherent  in  airborne 
and  other  Moving-base  tracking  systeMS. 

For  a  given  tracking  scenario,  the  physical  phenoMena 
of  target  dynaMics,  PLIR  systeM  vibrations,  and  atsK>spheric 
jitter  can  cause  apparent  Motion  between  the  target  and  the 
tracking  sensor,  where  apparent  iMplies  soMe  offset  froM  the 
actual  position.  Therefore,  the  location  of  the  centroid  of 
the  sensed  target  iMage  on  the  FLIR  detector  plane  is  a 
coMbinatlon  of  true  target  smtlon,  a  corruption  due  to  opti¬ 
cal  hardware  bending  and  vibration,  and  atMospheric  jitter 
caused  by  infrared  wavefront  distortion.  If  Xq  represents 
the  avlMUthal  offset  distance  of  the  apparent  target 
centroid  froM  the  center  of  the  PLIR  plane  POV,  then  the 
apparent  location  of  the  centroid  as  detected  by  the  PLIR 
array  is  as  follows: 
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Xc  -  Xt  ♦  Xa  +  Xfa 


(3-1) 


where : 

xt  ■  Xc  component  due  to  actual  target  dynamics 
Xa  *  Xq  component  due  to  atmospheric  litter 
xp  ~  Xc  component  due  to  bending/vibration 
and  all  components  are  measured  In  pixels. 

Since  the  PLIR  plane  Is  a  two-dimensional  array,  an 
equivalent  y-dlrection  relationship  to  Equation  (3-1) 
applies.  The  truth  model  that  describes  this  simulatloi 
corresponds  to  a  two-state  position  model  augmented  by  both 
a  six-state  atmospheric  jitter  model  and  a  four-state 
bending/vibration  model.  Subsequent  sections  will  describe 
the  dynamics  models,  follovrad  by  the  measurement  model, 
target  model,  and  the  simulation  model  which  allows  computer 
simulation  of  the  actual  tracking  scenario. 

3. 2  Dynamics  Model 

The  twelve-state  system  dynamics  model  can  be 
described  by  the  following  linear,  stochastic  differential 
equation: 

sr(t)  «  (3-2) 

where : 

Et  -  12  X  12  time-invariant  truth  model  system 
plant  matrix 

x<i>(t)  -  12-dlamnsional  truth  model  state  vector 
comprised  of  2  target  position  states, 

6  atmospheric  states,  and  4  bending/vibration 
states 
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^•P  ■  12  X  2  t iM-invar lant  truth  OKsdel  Input 
distribution  natrix 

y.<p(t)  *  2-dlBensional  truth  ntodel  deternlnlst ic 
input  vector 

j£>p(t)  ■  12-diMenslonal  truth  Model,  zero-mean, 

white  Gaussian  noise  vector  process  with 
autocorrelation  function: 


E{»r(t)wi.  (t  +  r))  - 


(3-3) 


The  solution  to  the  stochastic  differential  equation  in  a 
discrete-time  sense  takes  the  fora  (Notice  that  the  *T" 
subscript  denoting  truth  model  has  been  dropped  for 
convenience  on  the  following  developments): 


X(ti4.i)  -  ♦x(ti)  +  ^a<j(ti)  +  ifc|(ti) 


(3-4) 


where : 


x(t^)  <=  12-dimensional  discrete-time  state  vector 

±  «  state  transition  matrix  computed  from  fNp  over 
the  sample  period  ^t 

=  ti  +  i  -  ti 


ti+i 


l(ti+i,T)  irdT 


^(tj^)  2-diMnslonal  discrete-time  input  vector 

\j^(t]^)  =  12-diinenslonal  discrete-tinie,  zero-mean 

white  Gaussian  noise  process  of  covariance: 


tl4l 

/ 


tl+l/T)  ft  1  (ti+l,T)dT 


(3-5) 


S 
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Recall  that  the  twelve-state  discrete-tlM  truth  aodel 
vector  is  the  result  o£  augmenting  a  two-state  target 
dynamics  vector,  ^  six-state  atmospheric  jitter 

vector,  x^(tj^),  and  a  four-state  bending  and  vibration 
vector,  X5(ti).  The  partitioning  of  the  truth  model 
dynamics  system  sKkdel  solution  the  takes  the  following  form: 


*t(ti+l> 

it  2x2  1  ^2x6  1  ^2x4 

Xt(ti) 

a 

^x2  1  ia  6x6  1  ^x4 

0-4x2  1  0-4x6  1  ib  4x4_ 

Xb<*^i> 

Sdt  2x2 

02x1 

06x2 

Udt<ti>  + 

Sdat^l >6x1 

04x2 

Sidb<ti)4xl 

where : 

=  2-dimensional  target  dynamics  state  vector 

Xgi  «  6-dimenBional  atmospheric  jitter  state  vector 

^  E  4-dimen8ional  bending/vibration  state  vector 

ydt(^l)  ”  2-dimen8ional  discrete-time  deterministic 
input  vector 

iMa^^i^  “  6-dintenslonaI  discrete- time  white  Gaussian 
noise  related  to  atsK?spheric  states 


Wdb^^i)  *  4-dimensional  discrete-time  white  Gaussian 
noise  related  to  bending/vibration  states 


Realizing  that  the  block  diagonal  form  of  the  state  transi¬ 
tion  BMtrix  allows  independent  evaluation  of  each  particular 
Hiodel,  the  following  subsections  develop  the  propagation 
relationships  for  the  target,  ataospheric  jitter,  and  the 
bending/vibration  models  which  are  ultisately  incorporated 
into  Equation  (3-6)  for  the  entire  system  propagation 
relationship. 

3.2.1  Target  States .  The  truth  model  for  target 
dynamics  takes  the  form  of  a  continuous-time  deterministic 
model  representative  of  the  entire  range  of  target 
trajectories,  from  highly  sbaneuvering  to  very  benign.  Since 
a  PLIR  sensor  is  of  priaiary  interest  in  tracking  the  target, 
the  elevation,  a(t),  and  azimuth,  a(t),  angles  with  respect 
to  the  center  of  the  POV  on  the  two-dimensional  PLIR  image 
plane  become  the  measurement  parameters  by  which  true  target 
location  is  determined.  The  target  states,  therefore,  are 
propagated  forward  in  time  with  respect  to  the  flat  surface 
approximation  a-B  plane,  where  the  surface  is  actually  part 
of  a  huge  sphere. 

The  a-B  plane,  as  described  by  (7:36;  18:11-16;  19:26) 
consists  of  an  array  of  IR  sensing  pixels  and  is  perpendicu¬ 
lar  to  the  sensor-to-target  llne-of-sight  vector.  If  the 
target  is  far  from  the  PLIR  sensor,  the  PLIR  azimuth  and 
elevation  angles  become  proportional  to  the  linear  transla¬ 
tional  coordinates  Xf^  and  y^  of  the  target.  Pigure  3.1 
depicts  the  relationship  between  the  PLIR  angles  and  the 


•  :  TARGET 
Q  :  AZIMUTH 
3  :  ELEVATION 


Figure  3.1.  The  a-«  Plane 

target  centroid.  It  a  and  B  are  Measured  in  Micro-radians 

and  and  are  Measured  in  pixels,  then  the  proportions 

lity  constant  (kp)  is  the  angular  FOV  of  a  single  pixel 

which  has  been  20  Micro-radians/pixel  in  previous  efforts 

and  developed  by  Harnly  and  Jensen  (4:33). 

•  • 

If  a  and  B  are  considered  constant  over  the  filter's 
saMple  period  t  seconds  in  length,  then  the  change  in  the 
target  position  froM  one  sanple  period  to  the  next  can  be 
described  in  discrete  forM  by: 


Xt(tl^l) 

Ti  o1 

j  I 

xt(tj) 

kp^t  0 

«(ti) 

a  ' 

♦ 

* 

yt<ti^i ) 

1°  1  1 

yt(ti ) 

0  kp  lit 

B(ti  ) 

where : 


a(ti)  E  da/dt,  Meaaured  In  micro-radlans/second 
* 

S(tj^)  >  da/dt,  Measured  in  aiicro-radians/second 
kp  ■>  0.05  pixels/Micro-radian 

'  ^1  +  1  ' 


Equation  (3-7)  represents  the  discrete  fora  solution  of 
the  stochastic  differential  equation  for  the  truth  SKxlel 
target  states,  and  thus  describes  the  upper  partition  of  the 
truth  aodel  identified  in  Section  3.2.  This  type  of  model 
is  chosen  to  fit  the  state  space  model  form  similar  to  the 
rest  of  the  dynamics,  and  it  also  allows  for  easy  addition 
of  stochastic  driving  tersa,  if  desired. 

3.2.2  Atmospheric  States .  The  atmospheric  effects  as 
developed  by  Herder  (12:73,79)  for  both  the  x-  and  y- 
directions  of  the  PLIR  can  be  modeled  as  the  ouput  of  third 
order  shaping  filters  driven  by  white  Gaussian  noise  and 
represented  by  the  transfer  function: 

Xa(s)  ^2 


Wj|(s)  (s*tJi)  (Sfi^2) 


(3-8) 


where : 


Xa  ^  atmospheric  state  shaping  filter  output 
Wa  >  zero-MMan,  unit  strength,  white  Gaussian  noise 
Ka  *  atmospheric  RMS  adjustable  gain 
e  break  frequency,  14.14  rad. /sec. 
cjj  ^  double-pole  break  frequency,  659.5  rad. /sec. 


Phase  front  distortion  of  the  IR  radiation  waves  due  to 
atMspherlc  disturbances  causes  translational  shifts  on  the 
PLIR  isiage  plane,  called  "jitter”.  This  atmospheric  jitter 
perturbation  is  devoid  of  directional  dependencies; 
therefore,  the  modelinq  of  the  PLIR  plane  effects  in  both 
the  X-  and  y-  directions  can  be  treated  separately  and  iden¬ 
tically,  and  then  augmented  together  to  form  the  atmospheric 
system  description.  The  effects  can  be  described  by  the 
stochastic  differential  equation: 


(3-9) 


where : 


)^(t)  =  six-diMnsional  atmospheric  state  vector 

Ea  *  atmospheric  time-invariant  plant  matrix 

)fa(t)  *  six-dimensional,  zero-mean,  white  Gaussian 
noise  process  of  strength 


When  the  state  vector  and  the  plant  matrix  are 
expressed  in  Jordan  canonical  form  (12:212),  the  siatrix 
appears  as  (19:31): 


Qa 


Qi^ 

Q1O2 

01Q3 

Q1Q2 

Q2^ 

02Q3 

Q1Q3 

Q2Q3 

03^ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Qi^ 

Q1Q2 

Q1Q3 

Q1Q2 

02^ 

Q2Q3 

O1Q3 

Q2Q3 

03^ 

(3-10) 


where : 


Ol  ■  iKU}^CJ2^  )  / 
Ql  ®  ~0i 
Q3  * 


The  truth  aiodel's  center  partition  in  Equation  (3-6) 
is  £orned  by  representing  Equation  (3-9)  in  discrete  form. 
The  statistics  of  the  discrete-tiM  noise  covariance  can  be 
described  by  the  following  integral: 

t 

flda  * 


/ 


tat-rJOata 


(3-11) 


£^3  =  6X6  atnospheric  discrete-tiste  noise 
covariance  statrix 

^^(t)  =6X6  atmospheric  state  transition  matrix 
associated  with  of  Eq.  (3-9) 


and  the  6X6  atmospheric  state  transition  siatrix  forms  the 
center  partition  of  the  12  X  12  truth  model  state  transition 
aatrix.  The  atmospheric  state  transition  matrix  has  been 
developed  and  can  be  shown  to  be  of  the  form  (19:32): 


ia(4t) 


-CJl^t 

e  0  0 

-u>2At  -Ci>2-4t 
0  e  Ate 

-oj2At 

0  0  e 

0  0  0 

0  0  0 

0  0  0 


0 

0 


0 

0 


0 

0 

0 


0  0 

-u>lAt 

e  0  0 

-0J2At  ~<^At 
0  e  /4te 

-102  At 

0  0  e 


(3-12) 
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3.2.3 


States .  The  aodeling  of  the 


mechanical  bending/vibration  effects  is  based  on  a  report 
by  R  a  D  Associates  (15)  that  wias  prepared  for  APWL 
concerning  dynamic  SKMlelling  for  space  based  structures. 
Figure  3.2  depicts  the  response  spectrum  for  the  x-directlon 
line-of-sight  for  a  particular  satellite  structure  and  the 
approximated  curve  to  isiplement  a  reduced-order  model  that 
is  representative  of  this  type  of  response.  The  effects  due 
to  the  bending  and  vibration  of  the  Mchanical  structure  of 
the  optical  equipment  are  developed  in  Appendix  A.  Since 
the  effects  are  similar  in  both  the  x-  and  y-directions  of 
the  PLIR,  both  directions  can  be  modeled  identically  by 
swans  of  a  second  order  shaping  filter  driven  by  white 


Gaussian  noise  and  represented  by  the  transfer  function: 


Xb(s) 


where ; 


Wb(8)  8^  ♦ 


(3-13) 


Xb  «  swchanical  state  shaping  filter  output 

Wb  zero-mean,  unit  strength,  white  Gaussian  noise 
with  an  autocorrelation: 

B(Wb(t)Wb(t-T)  J  *  Qb  <)(t-T);  Qb  •  1  (3-14) 

Kb  *  gain  adjusted  to  obtain  the  desired  root  swan 
square  (RMS)  bend ing/vi brat  ion  output 

^  a  damping  coefficient,  >.15 

•  natural  mechanical  frequency,  >  n  rad. /sec. 

A  closer  approxiswt ion  to  the  low  frequency  spectra  of 
Figure  3.2  can  be  obtained  by  representing  each  direction  as 
a  sixth  order  system  versus  a  second  order;  however,  twelve 


Figure  2.2.  Response  Spectrum  for  LOS  x-directlon. 


states  versus  four  for  the  two  axes  was  deemed  as  beyond  the 
point  of  diminishing  returns  for  a  first-cut  analysis,  and 
the  modeling  of  the  fundamental  frequency  effect  of  bending 
%fas  determined  to  be  sufficient  enough  to  describe  the 
effect.  The  modeling  of  the  effects  in  both  the  x-  and  y- 


directions  can  be  treated  Independently  and  augnented 
together  to  forai  the  aechanlcal  systea  description.  These 
effects  can  be  described  by  the  stochastic  differential 
equation  of  the  following  fora: 

Xb<t)  “  Eb*b<t)  ♦  GfaWbCt)  (3-14) 

where : 

jc^lt)  s  f our -dlaens Iona  1  aechanlcal  state  vector 

Eb  "  tiate  invariant  aechanlcal  plant  aatrix 

w^(t)  s  two-diaenslonal,  2ero-aean,  white  (Gaussian 
noise  process  of  strength  Q5  I. 

a  4  X  2  noise  distribution  aatrix 

As  developed  in  Appendix  A,  the  Q5  Mtrix  is  shown  as: 


Qb  = 


(3-15) 


The  truth  aodel's  lower  partition  in  Equation  (3-6)  is 
foraed  by  representing  Equation  (3-14)  in  discrete  fora. 

The  value  of  the  d iscrete-t iae  Mchanical  noise  covariance 
can  be  described  by  the  following  integral  (9:171): 

Jit 

Obd  *  I lb(^>abftbSb'*^lb^(^)«lT  (3-16) 

b 


where : 

a  4  X  4  discrete-tiae  aechanlcal  noise  covariance 
aatrix 

f2}  4  X  4  Mchanical  vibration  state  transition 

aatrix  associated  with  £5  of  Bq.  (3-14) 


^t  -  saaple  period,  ti^i  -  t^ 


and  the  4X4  Mchanical  vibration  state  transition  SMtrix 


foras  the  lower  partition  of  the  12  X  12  truth  aodel  state 
transition  satrix.  The  Mchanical  state  transition  Matrix 
as  developed  in  Appendix  A  results  in  a  4  X  4  block  diagonal 
Matrix  of  two  identical  2X2  blocks  with  the  following  forn 


♦b  ■= 


where : 


♦bl 

♦b2 

0 

0 

♦b3 

♦b4 

0 

0 

0 

0 

^1 

♦b2 

0 

0 

^b3 

♦b4 

(3-17) 


=  exp(-a54t)(cos(Wb4t)  ♦  IsinCcj^^t )  ) 

#1,2  “  exp(-ab4t )  I  (IAl/j,  )8ln(i^b^^  )  J 

♦b3  “  oxp(-a5/lt)  {-(1  +  (054Jb>  ) ) 

♦b4  *  exp( -Ob^t)  Icoadt^b^t)  -  (0^/605) s in (w^^t )  J 

Ob  »  real  part  to  the  characteristic  equation  of 
Equation  (3-13) 

c^b  =  iaaglnary  part  to  the  characteristic  equation  of 
Equation  (3-13) 


3 . 3  Measurewwnt  Model 

The  BMasureaents  accessible  to  the  tracking  algorithm 
represent  the  average  intensity  of  target  and  background 
radiation  incident  on  the  PLIR  array  of  detector  elesMnts. 
The  measured  IR  iMge,  or  intensity  function,  consists  of  a 
collection  of  inforsMtion  from  the  target  IR  radiation. 


background  noises,  and  FLIR  sensor  noises.  Previous 
research  (16:11-20)  has  determined  that  single  target 
hotspot  patterns  could  be  modeled  as  a  bivariate  Oaussian 


function  with  constant-intensity  contours  when  describing 
distant  targets.  This  bivariate  Gaussian  intensity  function 
can  be  described  by  : 

1 1  X, y,  Xpg^l(  ( t )  ,ypQ3|(  ( t )  ]  =  lBkaix®*P^  ^ ^"’'peak  ^  ^y“ypeak  ^  ^ 

(Eafll  f  ^ *^“*peak ^  ^ y“ypeak ^  ^  ^  (3-18) 

Inhere : 

I^ax  ‘  naxiaiua  intensity  of  the  hotspot 

^ ^peak rypeak )  *  coordinates  of  the  centroid  of  the 

apparent  target  intensity  profile  on  the 
image  plane,  aieasured  in  pixels 

X/y  =  spatial  coordinates  calculated  relative  to 
the  center  of  the  tracker  POV. 

2 

Pflifl  *  2  X22  satrix  whose  eigenvalues  (Oy  and 
Opy  )  define  the  dispersion  of  the 
elliptical  constant  intensity  contours  in 
the  9-a  frasMt.  The  eigenvectors  of 
define  the  principal  axes  of  the  bivariate 
distribution. 

Although  single  hotspot  IR  sources  can  be  represented 
by  bivariate  Gaussian  distributions,  targets  with  multiple 
hotspots  cannot.  Figure  3.3  shows  the  PLIR  plane  intensity 
profile  for  a  three  hotspot  target  (14:35).  The  centroid  of 
the  apparent  target  intensity  profile  cannot  be  represented 
as  a  (Saussian  distribution  even  though  the  individual 
hotspots  are  modeled  as  Gaussian  with  elliptical  constant- 
intensity  contours.  Furthermore,  the  intensity  swasurement 
provided  by  each  pixel  corresponds  to  the  average  intensity 
on  that  pixel  resulting  from  the  effects  of  the  target 
intensity  function,  the  spatially  correlated  background 
noise,  and  the  FLIR  sensor  noise. 


Eq I 
Intensity 
C'oiitoix's 


ctxd  Armi^ 


'  'ei .  t  I’o  i  q  of 
Apjxd  I  tn  ,  L  1  j  1  c]e  I 
rit  ens  i  t  ^  In  cn'  i  1  e 


Figure  3.3.  FLIR  Isage  Plane  Intensity  Function 
for  a  Three  Hotspot  Target 


The  Measured  output  for  the  pixel  in  the  j-th  row  and  k-th 
column  at  sample  time  t^  is  described  by  (13:33,34): 


*jk^^i* * *  ”  (Il/ApI  Xpeakm^  ^i  ^ '^peakm^  ^i  ^  ^  *** 

m“l 

pixel  jk 


%#here  : 


+  nj)5  (tj  )  +  bjij  (tj  )  (  3- 

Zjk^^i)  *  output  of  pixel  jk 
Ap  -  area  of  one  pixel 
M  *  total  nusiber  of  hotspots 
I„(  . 1  ■  Intensity  function  due  tT  m-th  hotspot 
(x.y)  =  coordinates  of  a  point  within  pixel  jk 


(3-19) 


(^peakm^ypeakm)  =  coordinates  of  maximun  intensity  point 

of  B-th  hotspot 

~  effect  of  sensor  errors  on  pixel  jk 
bjkitf)  effect  of  background  noise  on  pixel  jk 

Temporally  and  spatially  uncorrelated  thermal  noise  and 
dark  current  in  the  IR  sensitive  pixels  comprise  the  sensor 
error,  nj|((t£)  (19:35);  whereas  the  background  noise, 
bjk(ti),  can  be  represented  as  a  spatially  correlated  noise. 
This  spatially  correlated  noise  is  modeled  as  an  exponen¬ 
tially  decaying  correlation  pattern  with  radial  symmetry  and 
characterized  by  a  correlation  distance  of  approximately  two 
pixels  (4:37,40).  This  is  simulated  by  otalntaining  non-zero 
correlation  coefficients  between  pixels  separated  from  each 
other  by  two  pixels  or  leas  in  all  directions  (19:35). 

3 . 4  Target  Model 

A  three  hotspot  planform,  with  a  coordinate  frame  that 
is  discussed  in  Section  3.5.1,  is  depicted  in  Figure  3.4. 

It  is  considered  to  represent  the  target  of  interest  for 
this  analysis.  The  target's  angle  of  attack  and  sideslip 
angle  are  assumed  to  be  zero  in  order  to  maintain  the  align¬ 
ment  of  the  semi-major  axes  of  each  of  the  three  elliptical 
IR  constant-intensity  contours  with  the  target's  velocity 
vector.  This  enables  a  simplistic  simulation  of  the  space 
geometry  without  substantially  degrading  tracker  performance 
analysis  (19:38). 
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Dlsplaceiftents  of  Hotspots  Proa  the  Mass  Centroid  Along 
©v/  Spv»  Sppv  Directions 


Hotspot 

fiv  ) 

Spv 

£ppv  ^  ) 

1 

1.0 

0.0 

0.0 

2 

0.5 

0.0 

3 

0.0 

-0.5 

0.0 

Figure  3.4.  Target  Hotspot  Distribution 

3.4.1  Target  Trajectories .  The  target  trajectories 
avaialble  in  this  study  include  the  four  trajectories  used 
by  Tobin  (19:38,41),  and  two  additional  target  scenarios  as 
shown  in  Figure  3.5.  Bach  is  described  as  follows: 

Trajectory  One.  This  trajectory  is  a  constant 
velocity,  straight  and  level  flight  with  a  siaulation  tiae 
of  five  seconds.  The  inertial  velocity  reaains  in  a  plane 
parallel  to  the  Cx-e^  plane.  Note  that  it  is  not  constant 
velocity  when  projected  onto  the  FLIR  iaage  plane  because 


the  tracker  is  a  rotating  fraae. 


Trajectory  Tvo .  This  trajectory  Involves  a  constant 
10-g  or  20-g  pull-up  suineuver  that  is  represented  by  a  step 
change  In  the  acceleration  at  t-2.0  seconds  and  lasts  for 
the  resiainder  of  the  5  second  sisnilation.  The  target's  turn 
is  portrayed  by  a  constant-smgnitude  angular  rate  vector 
oriented  parallel  to  the  negative  ^  axis. 

Trajectory  Three .  This  trajectory  is  similar  to 
Trajectory  Two  with  the  exception  that  at  t-3.5  seconds,  the 
constant-g  pull-up  Maneuver  abruptly  terninates  and  the 
pitch  rate  is  set  to  zero  for  the  remainder  of  the  6.67 
second  siHulatlon.  This  longer  simulation  interval  allows 
the  tracker  time  to  settle  down  and  corresponds  to  200 
sample  periods  of  duration. 

Trajectory  Four .  This  trajectory  is  similar  to 
Trajectory  Two  with  the  exception  that  the  angular  rate 
vector  is  parallel  to  the  Cy  axis,  which  simulates  the 
target  turning  towards  the  tracker. 

Trajectory  Five .  This  is  a  Trajectory  Two  rotated  at 
an  angle  of  45^  with  respect  to  the  FLIR  plane.  The  intent 
of  this  scenario  is  to  evaluate  how  well  the  tracking 
algorithm  performs  against  a  maneuvering  target  for  which 
the  tracking  algorithm  is  not  tuned  (the  appropriate 
direction  of  elongated  f ields-of-view,  to  be  discussed  in 
Section  5.4,  %rould  be  neither  the  x-  nor  the  y-  direction  in 
the  FLIR  image  plane). 

Trajectory  Six.  This  trajectory  is  Identical  to 


Trajectory  Two  with  the  exception  that  the  acceleration  to 


Figure  3.5.  Target  Trajectories 


start  the  pull-up  saneuver  is  allo%red  to  develop  as  a 
positively  increasing  raaip  rather  than  as  a  step  change. 
Previous  studies  had  BK>deled  the  accelerated  turns  with  a 
step  change  in  acceleration,  a  rather  harsh  uneuver  for  the 
tracking  algorithsi  to  staintain  lock  and  rather  unrealistic 
for  physical  targets  to  perform. 

3.5  Simulation  Space  Model 

The  sisKilation  space  model  allows  the  simulation  of  the 
PLIR  tracker  operation  on  a  digital  computer.  This  sK>del 
allows  for  the  realistic  target  trajectory  propagation  in 
three  dimensional  space  and  provides  a  siathematical  means  of 
describing  the  target's  image  and  velocity  vector  with 
respect  to  the  PLIR  image  plane  (19:36).  Both  the  PLIR 
plane  velocity  projection  and  the  PLIR  plane  target  image 
projection  are  discussed  following  the  description  of  the 
various  coordinate  fraaws. 

3.5.1  Coordinate  Prasies .  The  coordinate  frames  used 
for  the  PLIR  tracker  digital  simulation  include  the  target 
frame,  the  inertial  frame,  the  a  -  a  -  r  frame  and  the  a  -  a 
(PLIR  image)  frame.  The  origin  of  each  frame  and  the 
orientation  of  the  reference  axes  are  as  follows: 

Target  Praam ; 

Origin:  target  center  of  smss 

Axes:  g.v  -  along  velocity  vector 

-  out  right  side  of  target, 
perpendicular  to  e^ 


Epv 


Sppy  -  vector  completing  right  hand 

coordinate  set,  perpendicular  to 
both  Oy  and  e^y 

Inertial  Frame; 

Origin:  location  of  PLIR  sensor 

Axes:  -  due  north,  tangent  to  earth's 

surface,  defines  zero  azimuth 

ey  -  inertial  'up*  with  respect  to 
a  flat  earth  approximation 

Sjg  -  vector  completing  right  hagd 
coordinate  set,  defines  90” 
azimuth 

Note:  Azimuth  (a)  is  measured  east%«ard  from 

e^.  Elevation  (B)  is  Masured  from 
the  plane  defined  by  and  e^ . 

a  -  z.  L  Frame: 

Origin:  target  center  of  siass 

Axes:  ejr  -  coincident  with  the  true 

sensor-to-target  LOS  vector 

^,e^  -  define  a  plane  perpendicular 
to  gf,  rotated  from  the 
inertial  e^^  and  By  by  the 
azimuth  and  elevation  angles 

a  ^  a  (FLIR  Image )  Plane : 

The  PLIR  iaage  plane  by  which  the  sensor 
sieasureMnts  are  taken  comprise  this  plane. 

The  azimuth  and  elevation  angles  can  be 
considered  the  linear  translational  coordinates 
X  and  y.  Observing  the  PLIR  plane  along  the 
LOS  vector,  x  is  chosen  to  be  positive  to  the 
right  and  y  is  positive  do%m.  This  choice  of 
coordinates  maintains  a  right-handed  coordinate 
set  with  the  target's  range  measured  positive 
away  from  the  sensor. 


Both  the  inertial  and  a-B  fraMS  are  depicted  in  Figure  3.1, 
whereas  the  target  frasw  is  illustrated  in  Figure  3.3. 


3.5.2  PLIB  Plane  Velocit 


The  target 


inertial  velocities  are  transformed  into  the  deterministic 

•  • 
azimuth  velocity,  a(t),  and  elevation  velocity,  fi(t),  in  the 

truth  model  differential  equation.  Equation  (3-2),  by 

projecting  the  Inertial  velocities  onto  the  FLIR  image  plane 

based  on  the  geometry  of  Figure  3.6  which  %mis  utilized  by 

Loving  (8:27,29).  Prom  Figure  3.6(b),  a  can  be  derived  by 

the  relationship: 

a(t)  «=  tan"^  (2(t)/x(t))  (3-20) 


and  i9(t)  can  be  derived  from  Figure  3.6(c)  by: 

fl(t)  «  tan'^  (y(t)/rh(t))  (3-21) 

where : 

2  2  1/2 

r))(t)  >  horizontal  range;  >:  ( x  *  z  ) 

However,  the  variables  of  interest  for  velocity  projections 
•  • 

are  a(t)  and  a(t).  Therefore,  talcing  the  time  derivative  of 
Equations  (3-20)  and  (3-21)  yields: 


«(t ) 


X(t)V2(t)  -  Z(t)Vx(t) 
x^(t)  +  z^(t) 


and 


fl(t) 


rh(t)Vy(t)  -  y(t)rh(t) 


r^(t) 


(3-22) 


(3-23) 


where : 

Vx,  Vy,  V2  target  velocity  components  in  the  inertial 
«xf  ®y/  E*  directions,  respectively 


(a)  .  Target/Inertial  Frame  Geometry 


Azimuth  Geometry 


^  I  Target 

e  ,  e  :  inertial  axes  _  _  _ 

^  r  !  sensor-to-target  range;  =  (X  ^  y  +  z  ) 
rj^  :  horizontal  range;  =  (x^  +  z^)  '5 
V  :  inertial  velocity  of  target 
a  t  azimuth  angle 
3  j  elevation  angle 


KEY 


t  Sensor 


rh(t) 


x(t)Vx(t)  ♦  Z(t)V2(t) 
rh(t) 


(3-24) 


Equations  (3-22)  and  (3-23)  define  t)ie  detersiinist ic  input 

vector  in  the  truth  Model  difference  equation,  shown  in 

*  ’  T 

Equation  (3-6),  where  U^t^^i)  *  («( t ^ ) , A( t ^ ) ]  . 

3.5.3  FLIR  Plane  Target  I  sage  Projection .  The 
projection  of  the  target's  hotspots  onto  the  PLIR  isage 
plane  corresponds  to  the  the  detection  of  the  target  by  the 
array  of  pixels.  Although  the  hotspots  are  assusted  to  lie 
in  the  plane  forsad  by  the  wings  of  the  target,  the 
orientation  and  location  of  the  hotspots  on  the  PLIR  isage 
plane  changes  as  the  target  rotates  and  translates  with 
respect  to  the  IR  sensor.  Pigure  3.7  illustrates  the 
geoswtry  of  the  target  with  inertial  velocity  vector  assumed 


Pigure  3.7.  Target  Isage  Projection  Geosatry 


-j.  -  .  -  . .N  V  ..N 


N  S. 


to  project  through  the  nose  of  the  target,  the  «-B  plane 
perpendicular  to  the  true  LOS  fro*  the  tracker  to  the 
target,  and  the  origin  of  the  a-B  plane  coinciding  with  the 
detected  target  center  of  *ass .  The  inage  of  the  target  as 
perceived  by  the  IR  sensor  changes  as  the  target  changes 
angular  orientation  relative  to  the  tracker  or  noves  farther 
away  frost  or  closer  to  the  PLIR.  The  current  image  size  of 
the  target  with  respect  to  the  largest  planform  isiage,  or 
reference  istage,  is  given  by  the  following  (12:11-24): 

OTpy  “  ^pvo  f  *■  o^ J  (3-25) 

and 

Oy  ~  ( Zq /Z  1  t UpyQ  ♦  ( UyQ  —  O pyo  J  cosfl ]  ( 3— 26 ) 

*  Upyd  +  ( (Vj^los  )/vl  ( AR  -  D) 

where : 

Oyof  <^pvo  *  target  hotspot  dispersions  along  Cy  and  epy 
in  the  target  fraoie  of  the  reference  is»ge 

Oy,  Opy  «  target  istage 's  current  dispersions 

ro  ■  sensor -to-target  range  of  reference  isiage 

r  *=  current  sensor-to-target  range 

y  B  target  inertial  velocity  vector 

V  B  sagnitude  of  y 

yj_los  "  cosiponent  of  y  perpendicular  to  the  LOS 
vector  (projected  onto  the  «-a  plane) 

•  2  '2  1/2 

vjjos  •  ■aqnitude  of  Vxlos  “  +01 

«  angle  between  y  and  the  a-a  plane 
AR  -  aspect  ratio  of  reference  IsMge,  Uyg/Opyg 


Since  the  ultimate  goal  is  to  define  the  target  image  in 
terms  of  the  PLIR  plane  coordinates,  a  transformation  from 
the  target  frame  coordinates  to  the  a-B  coordinates  can  be 
expressed  as: 

=  T  X  (3-27) 

j^yJoB  jslnfi  cosej  j^y  target  frame 

where : 

8  angle  between  e^  and  Vxlos  Figure  3.7) 

The  dispersion  matrix  can  be  determined  from  the  following: 

EaB  ■=  I  E.  (3-28) 

where  the  dispersion  Mtrix  is  utilized  in  the  Masurement 
skodel  in  Equation  (3-18). 

3.6  Summary 

This  chapter  has  introduced  several  models  which  depict 
the  real  world  target  tracking  environamnt.  The  processes 
that  have  been  modeled  include:  target  dynamics, 
atmospheric  jitter,  mechanical  bending/vibration,  and 
bacicground  and  FLIR  noise.  Deterministic  target 
trajectories  were  described  to  provide  the  baseline  for 
testing  the  tracker  against  realistic  tracking  scenarios. 
Finally,  a  simulation  space  was  presented  in  order  to  allow 
the  simulation  of  the  tracking  scenario  on  a  digital 
coaq>uter . 


'-a  •j  /--  '>_V_  , 
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IV.  Filter  Models 


4 . 1  Introduction 

This  chapter  describes  the  filter  Models  esiployed  in 
the  HHAP  structure  discussed  in  Chapter  II.  Section  4.2 
discusses  the  reduced-order  dynaMics  Models:  the  Gauss- 
Harkov  acceleration  Model  and  the  constant  turn-rate 
dynaMics  Model,  and  Section  4.3  discusses  the  enhanced- 
correlator/1 inear-MeasureMent  Model . 

4 . 2  DynaMics  Models 

Two  distinct  dynamics  BK>del8  are  presented  for  the 
purpose  of  conparative  analyses  in  the  MMAP.  Both  models 
represent  a  reduced-order,  eight-state  vector  to  estlMate 
the  target's  position,  velocity,  and  acceleration,  and  the 
atMospheric  jitter  in  two  orthogonal  directions  (x,y)  on  the 
FLIR  image  plane.  The  filter  state  vector  for  each  of  the 
two  Models  is  represented  as: 


vhere : 


xt,yt  ■  target  position,  x-  and  y-direction 
Vx,Vy  =  target  velocity,  x-  and  y-direction 
*X'*y  *  target  acceleration,  x-  and  y-direction 
Xafya  ^  ataospheric  jitter,  x-  and  y-direction 

Note  that  the  ataospheric  jitter  is  represented  by  two 
filter  states  rather  than  six  and  that  the  filter  neglects 
the  aodelling  of  the  bend ing/vi brat  ion  states  as  presented 
in  the  truth  model  of  Chapter  III.  The  filter's  atmospheric 
jitter  oKKlel  captures  the  primary  characteristics  of  the 
jitter's  power  spectral  density  while  at  the  same  time 
neglecting  the  high  frequency  effects.  The  omission  of 
modeling  the  bending/vibration  states  within  the  filters  is 
be  presented  as  part  of  the  analysis  on  filter  performance 
in  Chapter  VI . 

Although  the  state  vectors  that  describe  the  two 
dynamics  models  to  be  investigated  are  the  same,  the 
description  of  the  target's  acceleration  process  is  uniquely 
defined  for  each  model.  The  linear  filter  of  Section  4.2.1 
NK>dels  the  target's  inertial  acceleration  as  a  zero-Man, 
first  order  Gauss-Harkov  process,  while  the  nonlinear  filter 
of  Section  4.2.2  models  the  target's  trajectory  as  described 
by  a  series  of  concatenated  constant  turn-rate  segments. 

4.2.1  Gauss-Markov  Acceleration  Model .  This  dynamics 
model,  a  standard  airborne  target  model  not  originating  in 
these  references,  assumes  that  the  target's  inertial 
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acceleration  can  be  described  as  the  output  of  a  first-order 
lag  driven  by  ifhite,  Oaussian  noise.  Developed  by  Millner 
(13:50,51)  and  Kozenchak  (7:75,80)  for  the  PLIR  tracker,  it 
is  well  described  by  the  tisie  invariant,  linear  stochastic 
differential  equation: 

jCf(t)  *  E.£lCf(t)  +  fi£%#f(t)  (4-2) 

where : 

2C£(t)  ■  8-diMnsional  filter  state  vector 

W£(t)  «  4-di»en8ional,  zero-mean,  white  Gaussian  vector 
noise  process  of  strength  Qf . 


0 


0 


0 


I 


where : 


'T'x'^y  =  correlation  tl»e  for  the  x-  and  y-  target 
acceleration  processes 

Tji  ■  correlation  tiise  for  atmospheric  jitter  processes 


2  2 
Ox  /Oy 


variance  and  mean-squared  value  for  the 
X-  and  y-  target  acceleration 


Oj|  «  variance  and  swan-squared  value  for  the 
atsKispherlc  jitter 


The  filter  state  estimate  and  error  covariance  swtrlx 
are  propagated  forward  In  a  dlscrete-tlsie  sense  by  the 
following  set  of  equations  (9:171,172): 


Jlf(ti  +  i  )  -  ♦{(  t)2C£(ti  ) 


(4-6) 


Ef(ti+i~)  -  tf(/lt)Ef(t/)tf'*‘(4t)  +  QLfd  (4-7) 


where : 


2C{(tj[)  =  filter's  estimate  of  the  state  vector 
Pf(ti)  *  filter's  error  covariance  swtrlx 


(ti  ) 


instant  before  the  measurement  is  incorporated 
at  time  t^ 


(t,*, 


instant  after  the  measurement  is  incorporated 
at  tiM  t| 


*f  (^t) 


tlM  invariant  state  transition  matrix 
associated  with  ££  for  propagation  over  the 
sample  period  t: 


^i+1  -  ti 


T.  T. 


tf(ti+i,T)fifa£gf*tf‘(ti+i,T)  dr 


The  ti(^t)  and  Q£(]  matrices  have  been  determined  to  be  of 


the  following  form  (14:47,46): 


0  #15 


^26  0 


If(Jt)  = 


♦46  0 


(4-8) 


♦55  0 


♦66  0 


♦77  0 


where : 


♦l5  *  ■  "^x^^  ~  exp{-4t/Tx)  )  ) 


♦26  ■  7-yl4t  -  Ty(l  -  exp{-4t/Ty}  )  J 


♦3  5  “  ^xll  -  expi-^t/T^l) 


m 
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#45  ■  Tyd  -  expl-4t/Ty)l 
#55  »  exp{-4t/T^) 

#66  “  ®xp{-4t/Ty} 

#77  ■  ^88  •  exp{-4t/7^} 


ftfd  = 


where : 


On 

0 

Ol3 

0 

Ol5 

0 

0 

0 

022 

0 

O24 

0 

026 

0 

Oi3 

0 

O33 

0 

035 

0 

0 

0 

024 

0 

O44 

0 

046 

0 

Ol5 

0 

035 

0 

055 

0 

0 

0 

026 

0 

046 

0 

066 

0 

0 

0 

0 

0 

0 

0 

O77 

0 

0 

0 

0 

0 

0 

0 

2/3)  TxZit^ 

-  2(Tx4tJ 

2 

4  Ty^^(At)ex 

(4-9) 


+  2r^^At  -  Tx^exp{-24t/Txl  +  Tx^) 


Q13  “  +  2  Tx^(4t)exp(-4t/rx)  ♦  '?'x^ 

-  2rx^exp(-zlt/rx)  -  2Ty^ At  +  Tx^expl -2zlt/Tx  1 ) 

0^5  *  ^x  ^x  ( At  )exp(  ~  A  ^  ^  *  *^x  expl— 2^t/’7"xl) 

022  “  ay^( (2/3)  Ty4t^  -  2(Ty4tl^  -  4Ty^ (4t )exp( -4t/Ty ) 

♦  2r^ At  -  Ty^exp( -24t/TyJ  +  Xy^) 


024  *  Vy  (TyAt  1  +  2Ty  (At )exp( -At/ TyJ  + 


y,  ^  y 


-  2  Tv^expC -^t/Tvl  '  2Tv^4t  ♦  Ty^expl -24t/Twl > 


026  ■  Oy^i-2Ty(At)Bxpl~At/Ty]  *  -  Ty^exp[-2At/Ty\} 

O33  ■  ~  4T^^expI -4t/T,j  1  -T^^expt -24t/r^  J } 

2 

O35  “  ®x  <'^x  ~  2T^expl -4t/r,{l  ♦  r^expE -2zlt/T^ ) } 

O44  “  Vy^ilTyAt  -  3Ty^  +  4Ty^eXpl  -^t/Ty)  “  Ty  ^6  X  p  {  ”  2 t  /  Ty  J  } 

2 

Q46  *  Cy  {Ty  -  2TyeXpI  -4t/Ty  1  +  Ty  C X P I  “ 2 ^ t / Ty  )  ) 

O55  ■  -  exp( -22lt/TxH 

066  =  Oy^il  -  expl-2At/Ty]} 

Q77  =  “  ®xpI-24t/TaI  > 

088  =  O77 

The  Multiple  Model  Adaptive  Filter  estinates  are  foraed 
via  a  probabilistically  weighted  averaging  technique. 
Following  each  propagation  cycle,  the  elesiental  filter 

A  A  A  A 

estimates,  xi(ti+i~)  »  xt(ti+i~)  and  X2(ti+i”)  -  yt(ti.H~) 

A  A 

form  the  MMAF  estimates,  xi  ( and  X2(tj^l  ^Hsaf’ 
Analogously,  the  estimates  X|nam£(t^^)  could  also  be  used. 
These  estimates  can  be  considered  as  the  control  signals 
that,  if  driven  to  zero  over  the  next  sample  period,  would 
drive  the  FLIR  optical  centerline  at  the  target. 

4.2.2  Constant  Turn-Rate  Dynamics  Model .  The  constant 
turn-rate  dynamics  model  is  an  alternative  to  the  linear, 
first-order  Gauss  Markov  acceleration  process,  and  it  is 
described  by  a  series  of  concatenated  constant  turn-rate 


•  > 


segments  (19:55,58).  Since  this  model  results  In  a 
nonlinear  dynamics  model,  the  extended  Kalman  filter 
propagation  equations  described  In  Chapter  II  are  utilized. 
The  constant  turn-rate  target  acceleration  process  Is 
modelled  by  the  following  first-order,  nonlinear,  stochastic 
differential  equation: 


■  2 

a(t)  =  -u  y(t>  ♦  w(t)  (4-10) 

whe  r  e : 

Cl)  =  angular  velocity  of  the  constant  turn 
y(t)  X  a(t) 


a(t)  =  2-dimen8ional  target  acceleration  vector; 
coordlnatlzed  In  the  PLIR  image  plane 

v(t)  *  2-dimenBional  target  velocity  vector; 

coordlnatlzed  In  the  PLIR  Image  plane 

w(t)  =  2-dimen8lonal  zero-mean,  white  (Saussian 
vector  noise  process 


This  model  results  In  a  nonlinear  system  state  equation 
of  the  form  of  Equation  (2-1)  when  utilizing  the  state 
variables  defined  in  Equation  (4-1).  The  components  of 
f.(x(t),tl  are  Identical  to  those  of  Equation  (4-3)  with 
exception  to  the  f5  and  fg  components.  These  two  components 
are  functions  of  the  nonlinear  vector  function  defined  in 
Equations  (4-10)  and  (4-11).  The  noise  distribution  matrix, 
Sf,  is  identical  to  Equation  (4-4)  and  the  noise  strength 
matrix,  (If,  can  be  determined  via  an  iterative  tuning 
process  depending  on  the  simulation  scenario  parameters. 
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The  state  estlatate  and  error  covariance  natrlx  are 


propagated  between  sample  periods  by  employing  the  extended 
Kalman  filter  equations  defined  In  Section  2.1.  The 
sensitivity  siatrlx,  of  Equation  (2-11)  Is  described  by 


the  partial  differential 

trlx 

equation: 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

ai 

0 

0 

0 

0 

0 

1 

0 

0 

£  *  * 

(4-12) 

dx 

0 

0 

P53 

F54 

F55 

P56 

0 

0 

0 

0 

P63 

^64 

^65 

^66 

0 

0 

0 

0 

0 

0 

0 

0 

-1/Ta 

0 

0 

0 

0 

0 

0 

0 

0 

-l/'^a 

where : 

A 

4x3^A 

F53  =  — 

—  ( 

-3X3X6 

+  X4X5  + 

1 

A  4X3X4A 

P54  ■  - 


(2x3x5  + 


1 


(-2X4Xe  ♦ 


B 


n3  -  — 


^64 


^65 


A 

-  0x4x5 

b2 

2x4^A 


-  X3X5  ♦ 


2 

4x4^A 


B 


Fee  = 


B‘ 

-2X3X4A 


B' 


A  =  (X3Xg  -  X4X5) 
2  2 

B  =  (X3'^  ♦  X4‘') 


4 . 3  Mea3uren>ent  Model 

The  neasureMent  Model  employed  corresponds  to  the  for* 
of  a  correlation  alqorlthM  developed  by  Rogers  (ie:S3,e3) 
which  provides  sieasureswnts  to  a  linear  Kalnan  filter 
•easureisent  Model.  This  algorlthM  Is  enhanced  over  the 
standard  correlation  algorlthM  by  the  following  (19:58,59): 

1.  The  current  PLIR  data  frame  Is  correlated  with  an 
estimate  of  the  target's  Intensity  function,  an 
adaptively  constructed  template,  as  opposed  to  the 
previous  frame. 

2.  Rather  than  simply  outputting  the  peak  of  the 
correlation  function,  the  enhanced  correlation 
algorithm  outputs  the  center  of  mass  of  that  portion  of 
the  correlation  function  which  Is  greater  than  some 
predetermined  lower  bound,  or  threshhold.  This 
technique  prevents  the  difficulty  of  having  to 
distinguish  global  peaks  from  local  peaks,  as  do  the 
peak-finding  algorithms  of  many  conventional 
correlator  algorithsis.  It  is  also  computationally 
less  burdensome  than  a  peak-finding  method. 


3.  The  PLIR/laeer  pointing  coHnands  are  generated 
according  to  the  Kalaian  filter's  propagated  state 
estlstate  rather  than  to  the  output  of  a  standard 
correlation  algorlthn,  thereby  incorporating  knowledge 
of  dynamics  over  the  ensuing  sample  period. 

The  resultant  outputs,  or  "pseudo-MasureMnts",  of  the 
enhanced  corellation  algorithm  are  the  linear  offsets, 
and  y^.,  of  the  centroid  with  respect  to  the  center  of  the 
f ield-of-view,  which  are  subsequently  incorporated  into  the 
linear  KalMn  filter  update  cycle.  k  description  of  this 
algorithm  is  presented  in  the  next  two  sections  and  a  more 
complete  description  is  discussed  in  (16:52,70). 

4.3.1  Generating  the  Tesu>late .  A  template  of  an  esti¬ 
mate  of  the  target's  intensity  profile  Is  created  over  the 
previous  "N"  centered  target  intensity  functions  (the 
shifting  property  of  Fourier  transforms  is  the  mechanism 
which  centers  the  intensity  functions  on  the  PLIR  plane 
before  averaging).  The  memory  size,  "N",  over  which  the 
intensity  functions  are  averaged  is  dependent  on  how  the 
shape  functions  are  varying  in  time.  Highly  dynamic  inten¬ 
sity  functions  destand  small  "M"  values  while  slowly  varying 
intensity  functions  admit  large  "N"  values. 

Since  a  viable  software  algorithm  concern  for  any 
application  is  the  minimization  of  memory  requirements,  the 
finite  stemory  averaging  is  approximated  through  the  use  of  an 
exponential  smoothing  technique.  The  exponential  smoothing 
technique  resembles  the  properties  of  a  finite  stemory 
filter  (10:33,39);  however,  it  requires  the  storage  in 


Mnory  o£  only  one  FLIR  data  fraine  and  not  "N”  frames 
explicitly.  The  template  for  the  exponential  smoothing 
technique  is  maintained  by: 

A  A 

I(ti)  -  Yl(ti)  ♦  (1  -  Y)I(ti_i)  (4-13) 

where : 

A 

I_(ti)  “  template-generated  "smoothed"  estimate  of  the 
target's  Intensity  function 

I.(tj^)  *  current  PLIR  data  frame's  "raw"  intensity 
function 

Y  =  smoothing  constant;  0  <  Y  <  1 

The  smoothing  constant,  Y,  in  the  exponential  smoothing 
algorithm  functions  as  does  "N"  in  the  finite  memory  filter. 
Large  values  of  Y  emphasize  the  current  data  frame  as  would 
a  corresponding  small  "N"  value.  A  smoothing  constant  of 
Y  =  0.1  had  previously  been  determined  to  suffice  for  this 
application  (18:V-12). 

Figure  4.1  shows  a  bloc)(  diagram  of  the  template 
generation  arrangement  enclosed  within  the  dashed  line.  The 
raw  FLIR  data  array  is  transforiaed  into  the  Fourier  domain 
through  the  use  of  a  Fast  Fourier  Transform  (FFT) 
implemented  using  the  Cooley-Tukey  FFT  technique  (7:18). 

The  data  array  is  then  centered  on  the  FLIR  plane  by  phase 
shifting  it  an  amount  equal  to  the  displacement  (in  the 
original  untrans formed  spatial  coordinates)  of  the  target 
and  atmospheric  position  states  from  the  center  of  the  FLIR 
FOV  by  the  following  shifted  amount: 
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Figure  4.1.  Enhanced  Correlator/Linear  Measurenent  Model 

*8hlft<ti)  =  (4-14) 

yshlft<ti)  -  X2(t/)^£  +  X8(t/)*„f  (4-15) 


The  shifting  property  of  Fourier  transforas,  which  states 
that  a  translational  shift  in  the  spatial  domain  corresponds 
to  a  linear  phase  shift  in  the  frequency  domain  enables  the 
shift  to  take  place  in  the  frequency  domain,  or: 


F(a(x-xahiftfy-y8hift>>  - 

£(  fx/ fy  )exp{  -  j2K(  fx*  Xghl  f  t  ^y*yahift)J 


(4-16) 


where : 


S.<x,y)  B  2-dlBensional  spatial  data  array 
P{ • }  B  Fourier  transform  operator 
g(£x/fy)  =  P{a(x,y)} 


After  centering  the  data  array  by  the  phase  shifting 
technique,  the  data  is  subjected  to  the  exponential 
smoothing  algorithm  of  Equation  (4-13).  The  "smoothed”  data 
is  then  stored  in  the  form  of  a  template  and  correlated  with 
the  subsequent  PLIR  data  frame  which  creates  the  "pseudo¬ 
measurements"  . 

4.3.2  "Pseudo-Measurements"  bv  Enhanced  Correlation . 
The  enhanced  correlation  algorithm  creates  the  "pseudo¬ 
measurements"  in  the  form  of  position  offsets  in  two 
orthogonal  directions.  These  offsets  represent  the  distance 
between  the  center  of  the  FLIR  FOV  and  the  centroid  of  the 
IR  iskage.  The  current  PLIR  data  frante  is  correlated  in  the 
Fourier  domain  with  the  template  generated  from  the  previous 
sample  time,  thus  resulting  in  a  cross  correlation.  This 
cross  correlation  is  performed  by  computing  the  Inverse  Past 
Fourier  Transform  (IFFT)  of  the  following: 


F{g.(x,y)*l(x,y)  I  =  G(fx,£y)lt 


(4-17) 


where : 


a( X, y ) *!.( x, y )  =  cross  correlation  of  g.(x,y)  and  l.(x,y) 

a(x,y)  =  measured  target  Intensity  function, 
current  PLIR  data  frame 

l.(x,y)  B  expected  target  Intensity  function, 
generated  In  the  form  of  a  template 


g(fx/ty)  -  P{q(x,y)) 

* 

t  <fx»ty)  “  complex  conjugate  of  Pli.(x,y)) 


After  performing  the  IPPT,  a  threshholding  is  performed 
on  the  correlation  function,  g.<x,y)*l.(x,y),  that  effectively 
zeros  out  any  value  that  is  less  than  30%  of  the  function's 
maximum  value  (14:52,55).  The  resultant  function  corres¬ 
ponds  to  the  relative  displacement,  or  offset,  of  the 
center-of-intensity  between  the  current  PLIR  data  frame  and 
the  template.  These  offsets  are  due  to  the  effects  of 
target  dynamics,  atmospheric  jitter,  and  measurement  noise 
(bending  effects  are  ignored  in  the  filters,  however,  the 
variances  of  the  measurement  noises  can  be  increased  during 
filter  tuning  to  account  for  this),  or: 


’^offset  “  ’^t  ’^a  •'x 
^offset  =  yt  ^  ya  ''y 


(4-18) 


(4-19) 


Equations  (4-18)  and  (4-19)  can  be  represented  in  state 


space  form  by: 


Z.(ti)  =  Hfx^(ti)  +  yf(ti) 


(4-20) 


where 


£(ti)  =  Ixoffset  /offset*  *  measured  in  pixels 


j5.f(ti)  =  filter  state  vector  of  Equation  (4-1) 


) 

(4-; 


VfCtj)  =  2-dimensional,  discrete-time,  zero-mean,  white 
Gaussian  measurement  corruption  noise  measured 
in  pixels,  o£  covariance  ^ 


The  measurement  noise,  Vfiti),  encompasses  the 
spatially  correlated  background  noise,  the  PLIR  sensor 
error,  errors  due  to  ignoring  the  vibration/bending  effects, 
and  errors  due  to  the  PPT/IPPT  processes.  The  covariance 
matrix  associated  with  this  error  (except  for  the  impact  of 
the  ignored  bending  effects)  has  been  determined  empirically 
to  be  (16:63,68)  : 


R  = 


0.00436  0 

0  0.00598 


(4-22) 


4 . 4  Summary 

This  chapter  has  formulated  the  filter  and  measurement 
models  for  the  MMAF's  elemental  filters.  One  filter  model 
represents  the  target's  acceleration  as  a  first  order  Gauss- 
Markov  process,  while  the  other  filter  itodel  treats  the 
target's  trajectory  as  a  series  of  constant  turn-rate 
segments.  Both  dynamics  models  utilize  a  common  linear 
measurement  model  that  generates  "psuedo-measurements"  via  a 
cross-correlation  technique  between  a  template  generated 
from  the  previous  data  frame  and  the  current  data  frame.  An 
advantage  to  the  linear  measurement  model  is  the  ability  to 
employ  the  linear  Kalman  filter  update  equations  (9:117,118) 
rather  than  algorithms  which  are  more  computationally 


burdensome . 


This  chapter  embodies  the  concepts  presented  in  the 


preceding  three  chapters  and  describes  the  tracking  algor- 
itha  eaployed  in  this  research  endeavor.  Following  the 
discussion  on  the  overall  tracking  algorltha,  the  aechanlsas 
exercised  In  evaluating  the  algorithm  are  also  discussed. 

5.2  Overview  o£  the  Tracking  Algorithm 

A  Bayesian  Multiple  Model  Adaptive  Filter  including  a 
bank  of  five  Independently  operating  Kalman  filters  provides 
the  basis  for  the  tracking  algorithm.  The  processed  fields- 
of-vlew  for  each  filter,  along  with  the  assumed  target's 
dynamics  in  the  FLIR  plane  directions,  appear  in  Table  5.1. 
Note  that  filters  II,  12,  and  13  are  represented  by  square 
f lelds-of-vlew,  whereas  filters  |4  and  15  are  rectangular 
f leld-of-vlew  filters. 


Table  5.1 

MMAF  Elemental  Filters 


Filter  # 

FOV  (Dlxels) 

Assumed 

x-dvnamlcs 

Assumed 

y-dynamlcs 

1 

8X8 

benign 

benign 

2 

24  X  24 

20  g 

20  g 

3 

8X8 

10  g 

10  g 

4 

24  X  8 

20  g 

benign 

5 

8  X  24 

benign 

20  g 

The  five  elemental 


5.2.1  Gauss-Markov  Model  Tracker . 

filters  In  the  Gauss-Harkov  MMAF  (GH  HHAP)  are  based  upon 

the  Gauss-Harkov  acceleration  model  of  Section  4.2.1.  Each 

filter  was  tuned  with  the  appropriate  direction-dependent 

2 

acceleration  variance,  a  ,  and  time  constant,  r,  (14:68). 
These  parameters  are  shown  in  Table  5.2  for  each  of  the  five 
filters.  Notice  that  the  parameters  pertaining  to  x-  and  y- 
directions  for  filters  #1,  12,  and  13  are  the  same,  whereas 
they  differ  for  filters  14  and  15  due  to  the  rectangular 
geometry  of  the  POV  of  these  filters.  Filters  14  and  15  are 
constructed  with  a  larger  f ield-of -view  size  in  the  direction 
in  which  harsh  maneuvers  are  expected.  The  filter  driving 


Table  5.2 

Gauss-Harkov  Filter  Parameter  Values 


Filter  1 

Tx 

2 

2 

Oy 

Ta 

At 

1 

4.0 

250 

4.0 

250 

.0707 

.2 

1/30 

2 

1.5 

10000 

1.5 

10000 

M 

N 

R 

3 

1.5 

2000 

1.5 

2000 

N 

R 

R 

4 

1.5 

10000 

4.0 

250 

N 

R 

R 

5 

4.0 

250 

1.5 

10000 

M 

R 

II 

Units : 

"^X  ^  "^y/ 

Ta,  z^t 

:  seconds 

a 

2  2 
X  '  ''y 

2  4 

:  pixels  /seconds 

2 

<^a 

2 

:  pixels 

noise  strength,  "0",  is  not  shown  but  can  easily  be  coaputed 
froa  the  following  relationship: 

flli  =  2a^/T  (5-1) 

where : 

2 

a  the  filter's  assumed  variance  for  that  target 

acceleration  process  or  jitter  position  process 

'<' =  correlation  tiM  of  target  acceleration  or  of 
jitter 

The  paraaeters  shown  in  Table  5.2  for  filters  12  and  13 
%irare  determined  by  tuning  a  single,  large  POV  and  small  POV 
Kalman  filter,  respectively.  Using  a  trajectory  12  target 
maneuver  (reference  Section  3.4.1),  filter  12  was  tuned  to 
20-g  dynamics  and  filter  13  was  tuned  to  10-g  dynamics. 
Parameter  values  for  filters  14  and  IS  are  chosen  such  that 
the  highly-varying  dynamics  coincide  with  the  appropriate 
channel  direction  of  the  rectangular  POV.  The  jinking 
maneuvers  correspond  to  the  elongated  direction  of  the  POV 
and  the  short  end  of  the  POV  is  characterized  by  the  benign 
dynamics  of  filter  II.  Also  note  that  the  sample  period. 

At,  for  past  research  efforts  has  been  30  Hz.  Section  6.3 
addresses  the  potential  for  implementing  the  tracking 
algorithm  at  a  50  Hz  rate. 

5.2.2  Constant  Turn-Rate  Model  Tracker .  The  MHAP 
structure  of  the  constant  turn-cate  (CTR)  model  tracker  is 
of  the  same  configuration  as  shown  in  Table  5.1.  The  key 
difference  between  the  CTR  and  GM  model  trackers  is  the 
strength  of  the  driving  white  noise  terms,  "0", 


as  shovm  in 


Table  5.3  fox  the  CTR  tracker  (the  computational  burden  of 
the  CTR  model  is  unnecessary  for  the  benign  dynamics  of 
filter  II).  Notice  that  the  "Q"  values  do  not  correspond  to 
the  relation  described  in  Equation  (5-1).  Previous  research 
has  shown  that  a  CTR  filter  is  considered  to  be  adequately 
tuned  to  a  trajectory  12  maneuver  when  the  peak  aiean-error 
is  approximately  1.5  tiates  the  filter  computed  error 
standard  deviation  or  root-mean-squared  (RHS)  value  at  a 
given  time,  rather  than  when  the  filter's  computed  RMS 
errors  match  the  actual  RHS  errors  connitted  throughout  the 
simulation,  as  evaluated  over  a  10-run  Monte  Carlo 
simulation  (10:69,71).  The  consequence  of  this  relationship 
is  that  the  CTR  filter  overestimates  its  own  errors  more 
than  does  the  OM  filter. 

Table  5.3 


CTR  Filter  Parameter  Values 


Filter 

1  Ox 

Qy 

2 

15000 

15000 

3 

4000 

4000 

4 

15000 

125 

5 

125 

15000 

Units : 

2  5 

pixels  /seconds 

The  f ields-o£-vlew  for  all  the  filters  are  represented 
as  8  X  8  data  arrays  In  order  to  exploit  the  single  generic 
correlation  algorithm.  This  is  accomplished  by  redefining 
the  pixel  size  in  filters  12,  14,  and  15  for  correlation 
processing  (filters  11  and  13  already  are  parameterized  in 
an  8  X  8  array).  Filter  I2's  POV  is  represented  as  an  8  X  8 
data  array  by  partitioning  the  24  X  24  array  into  an  8  X  8 
array  of  3  X  3  elements,  where  each  3X3  element  consists 
of  an  average  of  the  nine  "small"  (20  micro-radians  on  a 
side)  elements  from  the  larger  POV.  This  essentially 
simulates  a  data  array  of  a  "large"  (60  micro-radian)  pixel 
array.  The  redefined  pixel  for  filter  14  is  averaged  over  a 
3X1  array  of  "small"  pixels  due  to  the  rectangular 
geometry  and  measures  60  micro-radians  x  20  micro-radians. 

By  similar  argument,  filter  IS's  redefined  pixel  measures  20 
micro-radians  x  60  micro-radians.  The  POVs  and  unit  pixel 
sizes  (shaded  regions)  for  correlation  processing  are  shorn 
in  Figures  5.1  and  5.2. 

5. 4  Filter  Parameters 

At  Initialization  of  the  tracking  simulation,  accurate 
a  priori  knowledge  of  velocity  is  assumed  and  position  is 
acquired  by  an  algorithm  developed  by  Tobin  (19:130,139). 

The  position  states,  x^  and  X2,  sre  initialized  by 
positioning  the  target's  center  of  mass  in  the  center  of  the 
PLIR  POV.  The  velocity  states,  X3  and  X4,  are  initialized 
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dependent  upon  the  trajectory  selection  and  the  acceleration 
states,  X5  and  X5,  are  Initialized  by  computing  the  delta 
velocities  at  time  t,,  and  time  tj^  and  dividing  by  the  sample 
period.  The  atBK>spheric  jitter  states,  X7  and  xg,  are 
initially  zero.  The  error  covariance  matrix,  P(to),  is 
initially  reflected  as: 
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The  conditional  probabilities. 

Pl» 

P2/  P3»  P4 

,  and 

are  initialized 

to  0.96, 

0 

. 

0 

0.01 

,  0. 

01,  and  0.01  for 

filters  11 

thru 

15,  respectively. 

The 

initial  "heavy" 

iMighting  of  the  conditional  probability  on  filter  II  is  due 
to  the  "best-Btatching"  characteristics  of  filter  II 's  benign 
dynamics  to  the  target's  straight  and  level  trajectory  for 
time  less  than  two  seconds.  In  actual  practice,  initial 
acquisition  can  be  presumed  to  occur  on  a  benign  target. 

In  the  event  of  loss  of  track  due  to  gross  errors  in 


any  of  the  filter's  target  position  estimates,  the  tracking 
algorithm  executes  a  reacquisition  routine.  This  routine 


resets  the  state  vector  and  covariance  natrlx  of  the 


divergent  filter  as  a  linear  conbinatlon  of  the  state 
vectors  and  covariance  Matrices  of  the  nondivergent  filters. 
The  conditional  probabilities  of  the  nondivergent  filters 
are  then  scaled  to  sun  up  to  one  prior  to  evaluating  the 
state  vector  and  error  covariance  matrix  of  the  divergent 
filters  by  : 

-  N 

JSd  »  E  PnXn  (5-3) 

n*l 

Ed  =  E  PntEn  *  (Jin  “  Jid^Jin  “  3id»  *  (5-4) 

n  =  l 

where : 

JSd'Ed  ~  state  vector  and  error  covariance  matrix  of 
each  divergent  filter 

Jin'En  state  vector  and  error  covariance  matrix  of 
each  nondivergent  filter 

N  =  number  of  nondivergent  filters 

Pn  ®  hypothesis  conditional  probability  of  the  n-th 
nondivergent  filter.  The  p,^  values  are  scaled 
such  that:  Pi  P2  •••  Pn  ”  1*0 

5.5  Tracking  Algorithm  Statistics 

A  Monte  Carlo  simulation  provides  the  mechanism  for 
evaluating  the  tracking  algorithms'  performance.  Based  upon 
the  conclusions  of  previous  research  efforts  that  the  sample 
statistics  of  ten  Monte  (3arlo  simulation  runs  provides 
sufficient  convergence  to  the  actual  statistics  from  an 
infinite  number  of  runs  (19:60),  the  statistics  from  ten 


Monte  Carlo  runs  will  act  as  the  basis  upon  «rhich  tracker 
performance  is  evaluated. 


The  sample  mean  errors  of  the  tracking  algorithm's 
estiskates  are  computed  by  : 

_  M 

®xd^*-iJ  *  ®xdn^*-i^ 

n=l 

N 

=  l/N  E  IXdn(ti)  -  Xdnf(ti)l  (5-5 

n  =  l 

where : 

Bxd(^i)  ”  sample  mean  error  of  the  x  target  position 
estimate  at  tiste  tj[,  averaged  over  N  runs 

exdn^^i)  ~  error  in  the  MHAP  x-position  estimate  at  t^^ 
during  simulation  n 

A 

^nf(^i)  ~  HMAF  estimate  of  target's  x-position  at  t^ 
during  simulation  n 

^dn(^i)  ”  truth  model  value  of  the  x-position  of  the 
target  at  t^  during  simulation  n 

N  >  number  of  Monte  Carlo  runs 


The  sample  variance  of  the  error,  defined  in  terms  of  the 
same  paraswters  as  the  sample  mean  errors,  is  described  by: 

t^xd^tti)  =  (1/(N-1)J  E  exdn^<ti)  -  {N/{N-1 )  ti )  (5-6 

n“l 


Two  error  parameters  are  of  particular  concern  when 
evaluating  the  tracking  algorithm.  The  first  is  the  error 
committed  in  estismting  the  target's  position  in  both  the  x 
and  y-  PLIR  plane  directions.  This  parameter  serves  as  the 
fundamental  evaluation  mechanism  for  tracker  performance. 
The  second  parameter  is  the  estimation  error  committed  in 


locating  the  target's  centroid  of  the  FLIR  plane.  Because 
the  target's  centroid  location  is  paraaount  for  centering  It 
on  the  PLIR  plane  for  the  template  generation/  this  error 
parameter  supplies  Information  with  respect  to  the  filters' 
ability  to  reconstruct  the  target's  shape  function  and 
centroid  location. 

5.6  Performance  Plots 

Eight  performance  plots  are  utilized  In  this  research 
effort  to  evaluate  filter  performance.  The  performance 
plots  for  both  the  x-  and  y-dlrectlon  in  the  PLIR  plane  are 
as  follovrs  : 

1.  True  x-positlon  RHS  error  vs.  filter -computed 
x-position  RHS  error 

2.  True  y-posltlon  RHS  error  vs.  filter-computed 
y-posltlon  RHS  error 

3.  Hean  x-target  position  error,  ±  one  standard 
deviation,  plotted  at  tj^-  for  all  1 

4.  Hean  y-target  position  error,  ±  one  standard 
deviation,  plotted  at  t|^-  for  all  1 

5.  Hean  x-target  position  error,  ±  one  standard 
deviation,  plotted  at  tj'f  for  all  1 

6.  Hean  y-target  position  error,  f  one  standard 
deviation,  plotted  at  t^-t-  for  all  1 

7.  Hean  x-centrold  error,  ±  one  standard  deviation, 
plotted  at  t^t  for  all  1 

8.  Hean  y-centrold  error,  1.  one  standard  deviation, 
plotted  at  t^'t-  for  all  1 


Performance  plots  il  and  #2  provide  Inforauitlon  with 
respect  to  the  tuning  process,  plots  13  thru  16  provide 


tractclnq  perforstance  evaluation,  and  plots  17  and  18  provide 
Information  regarding  the  effectiveness  of  the  target 
intensity  template  identification. 


Illustrations  of  plots  82,  §4,  16,  and  18  are  shown  in 
'figures  5.3  thru  5.6,  respectively,  for  the  Oauss-Markov 
MMAP  algorithm  tracking  a  10-g,  trajectory  12  target 
maneuver  (note  that  the  expected  level  of  bending  is 
included  in  the  truth  model  but  not  in  the  filter  model). 
Bach  of  these  plots  portrays  the  evidence  of  the  target 
maneuver  at  t  =  2.0  seconds. 

Bach  performance  plot  is  labelled  with  a  plot 
designation  code  that  identifies  the  specific  simulation 
scenario  for  which  the  data  is  referrenced.  This  code  is  in 
the  following  format: 

field  1/field  2/field  3/field  4/field  5/field  6 

where  fields  15  and  16  are  optional,  and  the  fields  contain 
the  following  information: 


field 

1  ; 

OH 

; 

Gauss-Markov  filter  model 

CTR 

constant  turn-rate  filter  model 

field 

2  : 

HMAP 

multiple  model  adaptive  filter 

SMG-S 

single  filter  with  small  (8x8)  POV, 
used  to  establish  the  benchmark 

field 

3  : 

T(i) 

• 

identifies  type  trajectory  (11,  12, 
13,  14,  15,  or  16) 

field 

4  : 

(l)-G 

: 

identifies  the  magnitude  (in  g's) 
of  the  target  maneuver 

,  -  .  -  .  -  ,  -  .  -  .  -  .  •  .  •  -  -  ..  >■  -  •  _  »•  -  ^  ,•'  J’  \S  w*  ^  „  '  A  «  -  -  ■ 


field  5  ;  (Optional) 

T=  ’^i,'^2  ’  identifies  the  assumed  filter 

correlation  time  pairs  when  unlike 
the  nominal  (4,1.5) 

50  Hz  :  identifies  a  simulation  with  the 
50  Hz  sample  time  incorporated 


QB(I)  :  identifies  the  level  of  bending 
phenomena  included  in  the  truth 
model 

Zq  =  :  identifies  when  the  initial 

z-coordinate  is  different  from  the 
nominal  (20,000  meters) 

field  6  :  KP  =  :  identifies  the  different  pixel 

proportionality  constant  size  from 
nominal  (20  virad) 

R(4-)  :  identifies  when  the  measurement 

variance,  R,  is  different  from  that 
identified  in  Equation  (4-22) 

4* 

For  example,  the  plot  designation  codes  appearing  in 
Figures  5.3  thru  5.6  reflect: 


■/  i 


GM/MMAP/T2/10-G/QB1 

This  designation  implies  that  the  simulation  involved  the 
Gauss-Narkov  MHAF  tracking  model  with  the  expected  level  of 
bending  included,  against  a  trajectory  12  target  maneuver 
with  the  target  pulling  a  10-g  turn. 


Greyscales  are  symbols  used  to  chararterir^ 


ranges  on  FLIR  plane  images  and  th'  filter 
sysR>ols  are  shown  in  Table  5.4.  The  41  e.-  .  • 


M-MM  C19  n  MULTIPLE  HOKL  MOflPTtVE  TtflCKlNO  HLQORITHN  MAIHST 
aiRMMHE  TMMCTSIU)  AIR  FORCE  INST  OP  TECH 
HRIOMT-FATTERSOH  AFl  OH  SCHOOL  OF  EMQIMEERIMO 
UNCLASSIFIED  T  A  LEEHEV  DEC  07  AFIT/QE/EHQ/OTD-I?  F/Q  17/5. 1  ML 


which  is  a  pictorial  view  of  greyscale  symbols,  illustrates 
the  adaptive  identification  of  the  target's  intensity  shape 
function  in  the  form  of  a  tes^late  on  a  24-by-24  grid  of 
pixels.  The  exastple  shown  in  Figure  5.7  is  the  result  of  a 
coBiparison  between  a  template  and  a  noise-corrupted  PLIR 
■masurement  array.  The  spread,  or  arrangement  of  the 
greyscale  symbols  indicate  the  dispersion  of  the  hotspot 
about  its  centroid  and  also  the  location  of  this  intensity 
peak  relative  to  the  center  of  the  PLIR  POV. 


123456789012345678901234 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

_ +  +  +  - 

16 

17 

18 

19 

20 

21 

22 

23 

24 

Figure  5.7.  Noise  Corrupted  PLIR  Measurement  Array 


5.6  SuwBwry 

This  chapter  has  blended  the  ideas  and  concepts 
discussed  in  the  previous  three  chapters  into  a  viable 
multi-filter  tracking  algorithm.  As  a  basis  of  comparison, 
two  distinct  tracker  models  have  been  proposed  -  a  Oauss- 
Markov  acceleration  model  tracker  and  a  constant  turn-rate 
■K>del  tracker,  both  of  which  utilize  the  same  FOV  processing 
algorithms.  These  algor ithsis,  although  normally  implemented 
optically,  represent  the  FLlR-generated  information  in  the 
form  of  digitally  generated  intensities  that  can  be  por¬ 
trayed  via  greyscale  diagrams,  %rhich  serve  as  the  templates 
as  discussed  in  Section  4.3.1.  Finally,  the  tracker  statis¬ 
tics,  through  «rhich  the  evaluations  in  Chapter  VI  will  be 
analyzed,  %rere  presented. 


This  chapter  presents  the  results  o£  the  various 
analyses  for  those  topics  outlined  in  Section  1.3.  The 
evaluation  of  the  benefits  of  simulating  the  tracking 
algorithms  at  a  50  Hz  rate  versus  a  30  Hz  rate  are  discussed 
in  Section  6.2.  After  the  effects  of  the  bending/vibration 
states  are  analyzed  utilizing  a  single  filter  for  both  the 
Gauss-Markov  model  tracker  and  the  constant  turn-rate  model 
tracker  in  Section  6.3,  the  next  section  then  analyzes  the 
effects  of  bending/vibration  effects  with  the  Gauss-Markov 
multiple  model  adaptive  filter  algorithm.  A  sensitivity 
study  of  various  parameters  for  which  the  filters  are  not 
retuned  (other  than  ^ust  altering  that  specific  parameter  in 
the  filter)  is  presented  in  Section  6.5,  and  a  sensitvlty 
analysis  of  two  additional  target  trajectories  is  presented 
in  Section  6.6. 

6.2  50.  Hi  Implementation  Analysis 

As  suggested  in  Section  1.3.2,  the  performance  gains 
attributable  to  a  50  Hz  Implementation  are  Investigated. 

This  assessment  is  performed  by  first  establishing  a  bench¬ 
mark  of  perforamnce  for  both  the  Gauss-Markov  and  constant 
turn-rate  tracker  «K>dels  simulated  at  the  presently- 
configured  sample  rate  of  30  Hz.  After  the  benchmark  has 
been  established,  the  appropriate  soft%#are  changes  are 
implemented  to  accommodate  the  faster  sampling  frequency. 


and  an  Identical  simulation  scenario  as  conducted  with  the 
benchmark  runs  is  performed,  indicating  the  improvement  or 
degradation  as  a  result  of  sampling  at  a  higher  frequency. 

As  mentioned  in  Section  1.3,  Tobin's  results  (13)  form 
the  basis  of  the  tracking  performance  benchmark  for  the 
analyses  to  be  presented  in  this  chapter.  This  benchmark  is 
defined  as  the  performance  of  a  single  filter  which  is 
artificially  "told”  about  certain  dynamic  parameters  at  the 
onset  of  a  maneuver.  The  filter  employed  for  the  benchmark 
simulations  is  a  small  POV  (8x8)  filter  for  a  trajectory 
2,  10-g  maneuver.  The  changing  parameter  for  this  case  is 
the  strength  of  the  filter  driving  noise,  "Q",  which  is 
increased  prior  to  the  onset  of  a  maneuver  so  that  the 
filter -computed  RMS  error  reaches  a  steady-state  value  by 
the  time  the  maneuver  occurs,  thus  allowing  the  filter 
bandwidth  to  open  by  the  time  the  maneuver  is  initiated. 
These  simulations  are  conducted  for  both  the  Gauss-Harkov 
tracker  and  the  constant  turn-rate  tracker  smdels. 

6.2.1  Gauss-Markov  Mode  1 .  The  performance  plots  for 
the  benchmark  run  and  for  the  case  of  the  50  Hz  simulation 
are  shown  in  Figures  C-1  through  C-8,  and  C-9  through  C-16, 
respectively.  Notice  that  each  set  of  perforMnce  plots,  as 
is  true  of  all  subsequent  simulations,  is  composed  of  eight 
different  types  of  plots.  The  first  two,  as  in  Figures  C-1 
and  C-2,  show  the  f i Iter -computed  estimation  error  standard 
deviation  compared  to  the  actual  RMS  error  (in  pixels)  for 


both  the  X-  and  y-directions/  respectively.  These  two  plots 
indicate  the  tuning  adequacy  of  the  £ilter(s)  by  showing  a 
coa4>arlson  between  the  actual  and  the  fllter-coaiputed  RMS 
values.  Plots  three  and  four,  as  in  Figures  C-3  and  C-4, 
Identify  the  actual  estlauitlon  nean  error  t  one  standard 
deviation  at  tlM  »  ti  for  the  x-  and  y-directions, 
respectively.  These  parameters  are  used  as  a  basis  for 
control  generation.  Plots  five  and  six,  as  In  Figures  C-5 
and  C-6,  Identify  the  same  Information  as  plots  three  and 
four,  but  at  time  t^^.  These  parameters  provide  the  best 
possible  filter(s)  estimation  accuracy.  The  initiation  of 
the  target  maneuver  at  two  seconds  Into  the  simulation  is 
evident  in  Figures  C-4  and  C-6.  Plots  seven  and  eight,  as 
8ho%m  in  Figures  C-7  and  C-6,  show  the  actual  centroid 
position  estimation  error  mean  t  one  standard  deviation  at 
time  <=  tj^^  for  the  x-  and  y-directions,  respectively. 

These  parameters  provide  Information  regarding  the  adequacy 
of  iMge  centering  to  aid  in  the  construction  of  the  target 
image  template. 

In  both  the  30  Hz  and  the  50  Hz  case,  the  dynamics 
driving  noise  strength, "Q" ,  is  increased  at  frame  155  for 
the  30  Hz  simulation  and  at  frame  #195  for  the  50  Hz 
simulation,  allowing  for  a  five-frame  transient  prior  to 
target  maneuver  initiation.  Tables  6.1  and  6.2  show  the 
tesiporally  averaged  filter  error  statistics  (mean  and 
standard  deviation,  or  "sigma"),  the  filter's  peak-mean  y- 
positlon  estimate  errors  and  the  recovery  time  for  both 


cases  (x-posltlon  Is  just  as  Important;  however,  the  maneu¬ 
vers  evaluated  here  exhibit  harshly  changing  y-direction 
characteristics).  The  two  time-averaged  intervals,  0.5  to 
2.0  and  3.5  to  5.0  seconds,  are  chosen  to  allow  a  0.5  second 
minimum  for  transients  to  die  out.  These  intervals  are 
selected  to  show  comparable  statistical  data  during  the 
benign  region  and  during  the  maneuver.  The  peak-mean  error 
is  defined  as  the  greatest  mean-error  occurring  during  a 
given  simulation  and  the  recovery  time  is  defined  as  the 
time  it  takes  for  the  filter  to  recover,  or  reach  steady 
state,  from  the  maneuver.  Both  of  these  parameters  are 
precise  only  to  the  degree  that  the  values  are  discernible 
from  the  plots  (Figures  C-4,  C-$,  C-12,  and  C-14),  and  are 
considered  accurate  to  within  five  percent.  In  these  and 
subsequent  tables,  the  heading  designations  correspond  to 
the  convention  established  in  Section  5.6. 

In  comparing  the  temporally  averaged  mean  and  one  sigma 
values  for  the  two  tracking  simulations.  Tables  6.1  and  6.2, 
reveals  a  slight  improvement  in  the  oiean  error  associated 
with  estimating  the  position  and  an  improvement  in  the  true 
standard  deviation  about  that  mean  error,  which  is  expected 
since  the  algorithm  samples  the  data  at  a  faster  rate  and 
thus  allows  less  time  for  error  growth.  This  improvement  in 
error  standard  deviation  about  the  siean  can  be  seen  in 
Figures  C-3  through  C-6  and  C-11  through  C-14  by  comparing 
the  relative  displacement  of  the  one  sigma  curves  relative 


Table  6.1 


Single  Filter  Benchmark; 
GM/SNG-S/T2/10-G 


Temporally  Averaged 

TIm 

Interval 

Error  Parameter 

(mean  /  1  sigma) 

(0.5 

2.0] 

(3.5 

,  5.01 

Xerr (ti") 

-.1623 

/ 

.4089 

.2073 

/  .4912 

yert<ti') 

.0136 

/ 

.4241 

-.2588 

/  .4941 

*er r  ^  ^1  ^ 

-.1044 

/ 

.3765 

.1479 

/  .4290 

yerr^^l  ^ 

.0031 

/ 

.3876 

-.0998 

/  .4240 

x-centgrf (tj^) 

.0068 

/ 

.1226 

.0007 

/  .1546 

y-centgrr 

-.0083 

/ 

.0706 

.0577 

/  .0714 

y  peak- 

mean  error 

(ti') 

«  -2.0  pixels 

y  peak- 

Man  error 

(t/) 

>  -1.4  pixels 

recovery  tlM 

.  7  seconds 

Table  6.2 

Single  Filter  Simulation; 
GM/SNG-S/T2/10-G//50  Hz 


Temporally 

Averaged 

TIm 

Interval 

Error  ParaMter 

(Man  / 

1  sigma) 

(0.5 

9 

2.0) 

(3.5 

,  5.0) 

Xerr<ti" 

) 

-.0385 

/ 

.3779 

.1110 

/  .4557 

yerr  (*^l' 

) 

.0125 

/ 

.3787 

-.2352 

/  .4397 

’‘err  ^^1 

) 

-.0041 

/ 

.3568 

.0776 

/  .4172 

yerr 

) 

.0096 

/ 

.3522 

-.1486 

/  .3973 

x-centerr  <^1*) 

.0400 

/ 

.0889 

.0504 

/  .0911 

y-centgrr  <^1^) 

.0169 

/ 

.0795 

.0719 

/  .0815 

y  peak- 

Man  error 

(ti') 

■  -2.0  pixels 

y  peak- 

Man  error 

>  -1.3  pixels 

recovery  tlM 

*  .  6  seconds 

to  the  mean  curve  In  the  two  sets  of  plots  (also  note  the 
difference  in  scales  between  the  two  sets  of  plots). 

6.2.2  Constant  Turn-Rate  Model .  The  performance  plots 
for  the  two  simulations,  the  benchmark  run  at  30  Hz  and  the 
50  Hz  run,  are  sho%#n  In  Figures  C-17  through  C-24  and  C-25 
through  C-32,  respectively.  Notice  that  the  filter  tuning 
for  the  CTR  model  is  much  more  conservative  than  is  the  case 
with  the  GM  model.  Since  the  filters  with  the  CTR  dynamics 
model  exhibit  a  longer  filter-computed  covariance  matrix 
transient  than  do  the  GM  filters,  "Q"  was  increased  at  frame 
•35  (frame  1175  in  the  50  Hz  case)  to  allow  the  transient  to 
reach  steady  state  and  the  filter  bandwidth  to  open 
appropriately  prior  to  the  maneuver  initiation.  Tables  6.3 
and  6.4  show  the  statistical  data  for  the  two  simulations, 
which  show  results  that  are  comparable  to  those  of  the 
Gauss-Markov  tracker.  In  comparing  the  filter  estimated 
error  curves  for  both  the  simulations.  Figures  C-25  and  C-26 
reveal  that  the  relative  displacement  of  the  filters' 
estimated  errors  between  time  t^  and  t^^  is  larger  than 
what  is  observed  in  the  30  Hz  case  of  Figures  C-17  and  C-18. 
This  occurs  due  to  changing  the  continuous-time  dynamics 
noise  strength,  ”0”,  for  the  simulations  to  allow  the 
discrete-time  noise  covariance,  "Q(j",  to  remain  the  same. 

The  "0"  relationship  is  described  by  the  following  approxi¬ 
mation  inherent  to  the  CTR  model: 


%'  -.1 


■  »--  Tii  'T.  W-,  )r^  tTfc  Wc  K-i  w- 


Table  6.3 

Single  Filter  Benchmark; 
CTR/3NG-S/T2/10-G 


Temporally  Averaged 
Error  Parameter 
(mean  /  1  sigma) 

(0.5  , 

Time 

2.0) 

Interval 

(3.5 

,  5.0] 

Xerr (ti") 

.4106  / 

.3905 

.2343 

/  .3881 

yerr<ti~) 

-.0223  / 

.3822 

.3372 

/  .4232 

Xerr  ^  ^1  ^ 

.3406  / 

.3715 

.1872 

/  .3588 

yerr^^^l  ^ 

-.0242  / 

.3570 

.2987 

/  .3829 

x-cent^rr 

.0009  / 

.1145 

-.0203 

/  .1163 

y-centerr 

-.0156  / 

.0708 

.1347 

/  .0764 

y  peak- 

mean  error 

(ti~) 

=  -6.0  pixels 

y  peak- 

mean  error 

(ti^ 

=  -4.0  pixels 

recovery  time 

=  1.3  seconds 

Table 

6.4 

Single  Filter  Simulation; 
CTR/SNG-S/T2/10-G//50  Hz 

Temporally  Averaged 
Error  Parameter 
(mean  /  1  sigma) 

(0.5  , 

Time 

2.0] 

Interval 

(3.5 

,  5.0] 

Xerr^^l  ^ 

.5956  / 

.3160 

.1049 

/  .3985 

ye  r  r  ^  ^  1  ^ 

.0164  / 

.3069 

.3702 

/  .3618 

Xerr^^i  ^ 

.5461  / 

.3058 

.0744 

/  .3805 

yerrtti'") 

.0149  / 

.2943 

.3389 

/  .3435 

x-centerr 

.0084  / 

.0826 

-.0402 

/  .0919 

y-centerr 

.0232  / 

.0771 

.1962 

/  .0799 

y  peak- 

mean  error 

(ti") 

=  -6.0  pixels 

y  peak- 

mean  error 

(ti"^) 

=  -5.0  pixels 

recovery  time 

«  1.2  seconds 
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This  relationship  effectively  decreases  the  discrete-tine 


noise  strength  as  the  sample  frequency  is  increased, 
provided  the  continuous-time  noise  strength  is  unchanged. 
Although  conceptually  it  otakes  sense  to  keep  "Q"  constant 
rather  than  "Qf)",  "Q"  was  increased  to  maintain  the 
equivalent  discrete-time  filter  driving  noise  for  this  simu¬ 
lation  (to  Hiaintain  a  direct  comparison  between  the  two 
algorithms  at  the  different  sampling  frequency),  thus  reali¬ 
zing  a  higher  slope  between  t^  and  tj^^  for  the  filter- 
computed  estimation  errors.  Since  processing  at  a  50  Hz 
greatly  Increases  (73%  increase)  computer  processing  time, 
the  benefit  of  decreasing  the  uncertainty  of  the  position 
estimate  on  the  order  of  a  tenth  of  a  pixel,  as  compared  to 
the  30  Hz  version,  is  not  warranted.  Therefore,  in  an 
effort  to  conserve  computer  processing  time,  subsequent 
simulations  are  processed  with  a  sample  period  of  30  Hz. 

6 . 3  Single  Filter  Bend Inq/Vibrat ion  Analysis 

This  section  addresses  the  impact  of  including  the 
bending/vibration  states  of  a  large  space  structure  into  the 
tracking  algorithm's  truth  model  without  modifying  the 
tracker's  filter  model.  Thus,  this  is  a  robustness  study  of 
the  original  filter  design.  The  two  benchmark  scenarios, 
processed  at  30  Hz,  discussed  in  Section  6.2,  form  the 
baseline  of  comparison  for  both  the  Gauss-Harkov  tracking 
model  and  the  constant  turn-rate  tracking  model. 


6.3.1. 


Gauss-Markov  Model  Tracker .  The  performance 
plots  for  the  case  of  Including  the  bending/vibration  states 
are  shovm  In  Figures  D-1  through  D-8.  Notice  from  Figures 
C-1,  C-2,  D-1,  and  D-2,  that  the  filter  estimates  for  the 
simulation  with  the  bending/vibration  states  included  and 
without  additional  tuning  underestimate  the  actual  error. 
Although  the  true  RMS  value  of  x-  and  y-posltlon  has 
increased  with  the  inclusion  of  vibration  and  bending,  the 
filter-computed  error  does  not  appear  to  change.  This 
implies  that  the  filter  is  more  confident  of  its  outputs 
than  it  should  be.  Although  the  filter  can  be  improved  by 
retuning,  the  robustness  of  the  original  filter  is  the  prime 
concern  under  consideration.  In  comparing  the  time-averaged 
statistics  of  Tables  6.1  and  6.5,  it  is  evident  that, 
although  the  mean  error  between  the  two  cases  does  not  seem 
to  vary  much,  the  one  sigma  values  about  the  mean  have 
Increased  by  approximately  2/10  to  3/10  pixels  with  the 
inclusion  of  the  bending/vibration.  This  relationship  can 
also  be  seen  in  Figures  C-3  through  C-6  and  Figures  D-3 
through  D-6.  This  increase  in  the  uncertainty  with  respect 
to  the  position  variables  is  due  partly  to  the  filter  not 
estimating  the  bending  states  and  partly  to  the  enhanced 
correlation  algorithm  not  considering  the  effects  of  bending 
when  establishing  the  offsets  from  the  center  of  the  field- 
of-view  in  constructing  the  image.  Recall  the  measurement 
equations  for  the  two  directions: 


22  =  Yd  +  ya  +  Yb  ♦  V2 


(6-3) 


where  these  measurements  represent  the  location  of  the 
target  centroid  as  indicated  by  the  correlator.  However, 
the  filter  is  not  provided  any  information  with  respect  to 
the  two  bending  states,  X5  and  yj^,  nor  is  it  tuned  with 
additional  measurement  uncertainty  noise  associated  with  V]^ 
and  V2.  No  appreciable  change  in  estimating  the  location  of 
the  centroid  appears  between  the  two  simulations. 

Table  6.5 

Single  Filter  With  Bending/Vibration; 

GM/SNG-S/T2/10-G/QB1 


Temporally  Averaged 

Time 

Interval 

Error  Parameter 
(mean  /  1  sigma) 

(0.5 

,  2 

.0] 

[3.5 

,  5.0) 

^e  r  r  ^  ^  1  ^ 

-.1449 

/  . 

6297 

.1536 

/  .7625 

Yerr (ti") 

.0610 

/  . 

6407 

-.1588 

/  .6698 

*^e  r  r  ^  ^  i  ^ 

-.0849 

/  . 

5947 

.0922 

/  .7133 

Yerr <ti^) 

.0474 

/  . 

6026 

-.0015 

/  .6096 

x-cente^f  ( t|^  ) 

-.0057 

/  . 

1195 

-.0333 

/  .1296 

y-centgrr ^^i^J 

-.0055 

/  . 

0758 

.0617 

/  .0774 

y  peak- 

mean  error 

(ti') 

=  -2.0  pixels 

y  peak- 

mean  error 

(ti*) 

=  -1.2  pixels 

recovery  time 

=  .7  seconds 

*,  .  -w  ^  ■ 


6.3.2  Constant  Turn-Rate  Model  Tracker .  The 
performance  plots  for  the  case  where  the  bendlng/vlbratlon 
states  are  included  using  the  constant  turn-rate  model 
tracker  are  shown  in  Figures  D-9  through  D-16,  and  the 
statistical  data  appears  in  Table  6.6.  As  was  evident  in 
the  Gauss-Harkov  case,  although  the  time  averaged  x-  and  y- 
position  errors  are  roughly  of  the  same  order,  the  1-sigma 
deviation  for  the  bending/vibration  case  tend  to  be  .2  to 
pixels  greater  than  for  the  benchmark  case  in  Table  6.3. 
This  implies  that  the  inclusion  of  the  bending/vibration 
phenomena  increases  the  uncertainty  of  the  actual  position 


Table  6.6 

Single  Filter  With  Bending/Vibration; 
CTR/SNG-S/T2/10-O/QB1 


Temporally  Averaged 

Time 

Interval 

Error  Parameter 

(mean  /  1  sigma) 

10.5 

,  2. 01 

[3.5 

,  5.0) 

Xgr r  ^ ^i  ) 

.4511 

/  .6175 

.2088 

/  .7312 

yerr^*-i  ^ 

.0098 

/  .6179 

.4558 

/  .6266 

Xgrr  ^ ^i  ^ 

.3800 

/  .5877 

.1636 

/  .7000 

yert^*^i  ^ 

.0109 

/  .5884 

.4162 

/  .5872 

x-centgr j  ) 

.0022 

/  .1088 

-.0129 

/  .1109 

y-centgrr  (ti"*^) 

.0025 

/  .0647 

.1584 

/  .0622 

y  peak- 

mean  error  (t^  ) 

=  -6.0  pixels 

y  peak- 

mean  error  (t^*) 

=  -4.0  pixels 

recovery  time  =  1.3  seconds 


of  the  target  by  2/10  to  3/10  of  a  pixel,  or,  the  impact  of 
nodellng  the  bending  phenomena  has  a  the  same  impact  on 
the  constant  turn-rate  tracker  model  as  was  evident  on  the 
Gauss-Markov  model  tracker.  Again,  as  was  the  case  for  the 
Gauss-Markov  model  tracker,  the  time  averaged  statistics 
pertaining  to  the  location  of  the  centroid  are  not  affected 
by  the  inclusion  of  the  vibration  and  bending  phenomena. 

The  effects  of  performing  the  simulation  at  the  faster 
sampling  time  (50  Hz)  and  the  effects  due  to  modeling 
bending  and  vibraion  have  shown  comparable  results  to  both 
the  Gauss-Markov  model  tracker  and  the  constant  turn-rate 
model  tracker  (performance  data  not  shown).  Although  the 
constant  turn-rate  model  tracker  has  exhibited  better 
performance  than  the  Gauss-Markov  model  tracker  at  close 
ranges,  the  robustness  issues  to  be  addressed  do  not  concern 
close  range  scenarios;  therefore,  subsequent  simulations 
utilize  the  GM  MMAF  tracking  algorithm  since  past  studies 
have  shown  it  to  be  preferable  to  the  CTR  model  tracking 
algorithm  at  moderate  to  long  ranges. 

6 . 4  MMAF  Bending/ Vi brat  ion  Anal vs  is 

An  analysis  comparable  to  that  of  the  single  filter 
case.  Section  6.3,  is  performed  with  the  GM  MMAF  tracking 
algorithm.  In  addition  to  the  two  scenarios  analyzed  in  the 
single  filter  case  (excluding,  then  including,  the  bending 
and  vibration  phenomena),  two  additional  scenarios  are 
considered.  The  four  scenarios,  which  are  distinguishable 
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by  the  Magnitude  o£  the  zero  frequency  PSD  value  of  bending 

2 

noise  (see  Appendix  A  for  the  full  development  of  /  which 

is  in  fact  that  PSD  value  at  zero  frequency),  are  addressed 

2 

in  subsequent  sections:  1)  K5  =  0,  or  effects  of  bending/ 
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vibration  states  are  excluded;  2)  K5  =  5  X  10  ,  includes 

the  anticipated  effect  of  bending/vibration  without  filter 

2  -12 

retuning;  3)  K5  =  5  X  10  ,  Increases  the  effect  of 

bending/vibration  states  by  an  order  of  magnitude  above  the 
anticipated  nominal  values,  again  without  filter  retuning; 
and  4)  same  as  (3)  plus  the  addition  of  filter  tuning. 

These  analyses  are  performed  with  the  Gauss-Markov, 
multiple-model  adaptive  filter  trac)cing  algorithm  against  a 
trajectory  12,  10-g  maneuver. 

6.4.1  MMAF  Benchmark  Simulation.  Figures  E-1  through 
E-6  show  the  benchmark  of  performance  simulation  for  the 
case  where  the  bending/vibration  states  are  excluded  from 
the  truth  model.  These  results  are  the  same  as  those 
determined  by  Tobin  (19:150-157).  The  statistics  for  the 
HMAF  benchmark  are  shown  in  Table  6.7  and  the  dominant 
elemental  filters  are  reflected  in  Table  6.8.  These  two 
sets  of  information  serve  as  the  basis  of  performance 
comparison  for  subsequent  analyses  involving  the  Gauss- 
Harkov,  Multiple  Model  Adaptive  Filter,  10-G  pull-up 
maneuver,  trajectory  #2  simulations. 

Notice  that  the  reacquisition  algorithm  discussed  in 
Section  5.4  is  accomplished  at  frames  174  and  #75,  as  shown 
in  Table  6.6.  The  state  estimate  of  filters  II,  13,  14,  and 
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Table  6.7 


Multiple  Model  Adaptive  Filter  Benchiaark; 
OM/MMAF/T2/10-G 


Temporally  Averaged 

Time 

Interval 

Error  Parameter 
(mean  /  1  sigma) 

(0.5 

/ 

2.01 

[3.5 

/  5 

.01 

Xerr <ti~) 

.0024 

/ 

.4351 

.2394 

/  . 

5155 

yerr^^i  ^ 

.0001 

/ 

.4327 

.0968 

/  . 

7136 

^r r  ^  ^i  ) 

.0116 

/ 

.3852 

.2010 

/  . 

4303 

Yerr  (ti'*’) 

-.0090 

/ 

.3950 

.2311 

/  . 

6094 

x-centgjr 

-.0033 

/ 

.1123 

.0852 

/  . 

1787 

y-centg^r (ti^) 

-.0159 

/ 

.0634 

.3991 

/  . 

2781 

y  peak- 

mean  error 

(ti~) 

»  -1.7  pixels 

y  peak- 

mean  error 

(t/) 

=  -.7  pixels 

recovery  time 

-  .  6  seconds 

Table  6.8 

Dominant  Elemental  Filters; 
GM/MMAF/T2/10-G 


Interval 

Dominant 

( Frames ) 

Filter (s) 

Comments 

1-65 

11 

Tracking  of  benign  trajectory; 
maneuver  begins  at  frame  60 

66  -  73 

12,  15 

Y-direction  maneuver  recognized 

74  -  75 

12 

Other  filters  "lose  lock" 

76  -  120 

12,  14 

Wide  FOV  tracking 

121  -  150 

•  3 

10-g  maneuver  recognized 

15  are  all  reset  to  natch  filter  12*8  state  estimate  since 
filter  12  Is  the  only  filter  that  doesn't  "lose  lock".  At 
frame  174,  a  good  estimate  of  the  target's  y-accelerat ion  is 
predicted  by  filter  12,  which  is  then  transferred  to  the 
other  filters  during  the  reaquisition  cycle.  This  transfer 
effectively  resets  the  conditional  mean  value  of  the  Gauss- 
Harkov  acceleration  processes  in  the  remaining  filters,  thus 
providing  each  filter  an  accurate  estimate  of  actual  y- 
direction  target  dynamics.  It  is  reasonable,  therefore,  to 
stretch  the  POV  in  the  x-dlrection  (weighting  filter  14 
appropriately)  during  this  period  since  most  of  the  changes 
in  the  target's  acceleration  now  occur  in  the  x-dlrection. 

6.4.2  Effects  of  Bend ino/Vibrat ion.  This  section 

addresses  the  Impact  of  including  the  bending/vibration 

states  in  the  simulation  truth  model  but  not  telling  the 

filter  any  additional  inforntation.  Recall  from  Appendix  A 

that  the  strength  of  the  bending/vibration  white  noise,  Q|,, 

2 

is  expressed  in  terais  of  the  parameter  and  that  the 
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baseline  value  for  K5  is  5  X  10  (equal  to  the  zero- 
frequency  PSD  value).  Figures  E-9  through  E-16  show  the 
performance  plots  for  this  simulation.  Figures  E-9  and  B-10 
show  that  the  actual  error  exceeds  the  filter-computed  error 
at  random  intervals  throughout  the  simulation,  thus  implying 
that  the  filter  is  oblivious  to  the  changes  being  contrib¬ 
uted  from  the  bending/vibration  phenomena,  as  expected.  As 
was  the  case  for  the  single  filter  analysis,  the  MHAF 
appears  to  be  more  confident  of  its  outputs  than  it  should 


Jl 


Temporally  Averaged 

Time 

Interval 

Error  Parameter 
(mean  /  1  sigma) 

10.5 

f 

2.01 

(3. 

5  , 

5.0) 

Xerr<ti‘) 

.0410 

/ 

.6374 

.1837 

/  . 

7821 

Ve  r  r  ^  ^  i  ^ 

.0582 

/ 

.6463 

.0944 

/  . 

7983 

*er r  ^ ^i  ^ 

.0498 

/ 

.5926 

.1395 

/  . 

7224 

Ferr 

.0439 

/ 

.6076 

.2306 

/  . 

6912 

x-centg^r  <  t^i^ ) 

.0010 

/ 

.1332 

.0451 

/  . 

1471 

y-centgrr  (ti"*^) 

-.0081 

/ 

.0718 

.3143 

/  . 

2415 

y  peak- 

mean  error 

(ti') 

-  -1.6  pixels 

y  peak- 

mean  error 

iti*) 

*  -.7  pixels 

recovery  time 

»  .6  seconds 

Table  6.10 

Dominant  Eleoiental  Filter 
OM/MMAP/T2/10-G/QB1 


Interval 

Dominant 

( Frames ) 

Filter (s) 

Comments 

1-66 

•  1 

Tracking  of  benign  trajectory; 
maneuver  begins  at  frame  60 

67  -  72 

12,  15 

Y-direction  ataneuver  recognized 

73 

•  2 

Other  filters  "lose  lock" 

74  -  138 

12,  #4 

Vide  FOV  tracking,  intermittent 
attempts  to  lock  on  x-direction 
filter  (14) 

139  -  150 

13,  14 

10-g  maneuver  recognized,  harsh 
changes  in  x-direction  followed 

be.  Table  6.9  sho%rs  the  tiine  averaged  statistics  for  this 
simulation.  A  comparison  of  the  x-  and  y-direction  error 
statistics  of  Table  6.9  to  those  of  the  benchmark  case  of 
Table  6.7  reveal  that,  although  the  estimate  of  the  mean 
error  In  both  cases,  with  and  without  bendlng/vibratlon,  are 
very  similar,  the  uncertainty  about  that  Man  is  increased 
by  an  average  of  0.2  pixels  when  the  bendlng/vibratlon 
phenomena  Is  introduced.  This  is  a  direct  correspondence 
with  that  found  In  the  single  filter  benchmark  case.  This 
Increase  In  the  uncertainty  can  also  be  seen  by  comparing 
Figures  E-11  through  B-14  to  the  benchmark  equivalent  plots 
of  Figures  B-3  through  B-6.  The  estimate  of  the  x-  and  y- 
centrold  errors  and  their  associated  one  sigma  values  for 
both  cases  show  no  apparent  dissimilarities  between  the  two, 
as  was  the  case  with  the  single  filter  analysis  In  Section 
6.3.  The  dominant  elemental  filters  for  this  simulation  are 
shown  In  Table  6.10.  As  compared  to  the  benchmark.  In  Table 
6.8,  the  two  cases  reveal  some  differences.  First  of  all, 
this  tracking  scenario  tracks  the  benign  portion  of  the 
simulation  for  one  frame  longer  than  did  the  benchmark, 
otherwise  Interpreted  as  an  additional  one  frame  delay 
before  a  sianeuver  was  recognized.  Secondly,  and  also  sur¬ 
prisingly,  the  simulation  "lost  track”  for  only  one  frame 
(rather  than  two  as  did  the  benchmark)  shortly  after  recog¬ 
nizing  the  y-jlnk  ntaneuver.  And  thirdly.  It  takes  the 
algorithm  longer  to  recognize  the  10-g  maneuver  when  the 
bending  phenomena  Is  siodeled. 


section  addresses  the  concern  of  what  Magnitude  of  filter 
perforaance  degradation  can  be  expected  if  the  phenomena  due 
to  bending/vibration  is  an  order  of  swignitude  greater  than 
expected,  or  =  5  X  10  Figures  R-17  through  E-24 

show  the  perforawince  data  and  Table  6.11  reflects  the 
statistics  for  this  simulation.  Figures  B-17  and  E-16 
reveal  the  adverse  impact  as  a  result  of  the  increased 
effect  of  the  bending/vibration  phenoatera.  Although  the 
actual  RMS  error  has  increased  significantly,  the  filter- 
coatputed  RHS  error  ignores  these  effects.  This  is 
comparable  to  the  case  of  a  finely-tuned  filter  for  a 
scenario  which  is  really  not  present,  which  is  expected 
since  the  bending  phenomena  is  not  BK>deled  in  the  filter. 

The  temporally  averaged  statistics  of  Table  6.11  show  a 
significant  difference  in  the  one  sigaia  values  about  the 
Oman  x-  and  y-error  parameters,  although  the  aiean-error 
estiaates  had  not  changed  significantly.  The  standard  devi¬ 
ation  of  the  errors  has  increased  by  a  factor  of  four  to 
five  times  beyond  what  was  seen  in  the  benchaiark  simulation. 
This  wide  spread  of  the  one  sigata  error  is  clearly  visible 
in  Figures  E-19  through  E-22.  Since  the  average  uncertainty 
associated  with  the  x-  and  y-errot  paraaMters  lie  in  the 
range  from  1.5  to  2.0  pixels,  filter  tuning  is  appropriate 
in  an  attempt  to  force  the  filter  to  "recognize"  the 


changing,  dynamic  environment  which  It  is  ignoring,  again. 


A 


Table  6.11 


HHAF  with  Increased  Bend Inq/Vibratlon 
GM/MMAF/T2/10 'J/0B2 


Temporally  Averaged 

Time 

Interval 

Error  Parameter 
(mean  /  1  sigma) 

(0.5 

9 

2.0] 

(3.5 

,  5.0] 

*er r  ^  ^i  ^ 

.0927 

/ 

1.539 

.1294 

/  1.976 

Yerr (ti") 

.1610 

/ 

1.515 

.2304 

/  1.625 

Xerr 

.1004 

/ 

1.915 

.0741 

/  1.915 

Verr^^i  ) 

.1397 

/ 

1.449 

.3446 

/  1.528 

x-centgrr  <ti^) 

-.0345 

/ 

.1493 

.0053 

/  .1570 

y-centerr  <ti^) 

.0062 

/ 

.0932 

.1862 

/  .1976 

y  peak- 

mean  error 

(ti‘) 

=  -1.7  pixels 

y  Peal«- 

mean  error 

(ti^ 

=  -.7  pixels 

recovery  time 

B  .6  seconds 

which  is  expected  due  to  the  original  tuning.  Even  though 
the  filter  grossly  underest inates  the  error  variances  in  its 
estinates  of  the  target's  x-  and  y-position,  the  x-  and  y- 
centroid  errors  do  not  appear  to  be  affected  as  compared  to 
the  benchmark  simulation.  This  shows  that,  although  the 
state  estiskates  are  corrupted  via  the  reduced-order  filter 
model,  there  is  little  impact  on  centering  the  image  for 
averaging . 

6.4.4  Increased  Bend ino/Vi brat  ion  With  Tuning .  This 


section  addresses  the  attainable  performance  when  the  order 
of  magnituoe  of  bending/vibration  phenomena  such  as  shown  in 


Section  6.4.3  Is  evident  in  the  truth  nodel  but  the  filter 


is  allovred  to  be  retuned  for  best  estimation  accuracy.  In 
selecting  the  appropriate  filter  tuning  parameters,  the 
filter  "Q"  is  disregarded  since  the  dynamics  models  have  not 
varied  from  the  baseline  simulation.  The  filter  "R"  value, 
however,  should  vary  due  to  the  fact  that  the  bending  and 
vibration  phenomena  was  not  modeled  in  the  filter  (recall 
Equations  (6-2)  and  (6-3))  and,  in  fact,  the  "R"  is  expected 
to  increase,  to  account  for  the  variance  (RMS)  contribution 
of  xj,  and  y^  to  and  22/  respectively.  An  initial 
approximation  as  to  what  "R"  value  to  choose  can  be  computed 
from  the  adaptive  estimation  relation  (10:122): 


B.(ti)  -  E  -  H(t j)P(t j")H^(t j)  1  (6-2) 

W  j=i-N+l 


The  first  approximation  of  "R"  revealed  that  it  should 
be  Increased  almost  two  orders  of  magnitude  greater  than  the 
empirically  determined  values  identified  in  Equation  (4-22) 
without  bending  effects  considered.  This  increase,  however, 
only  improved  tracking  performance  slightly  better  than  what 
had  been  shown  without  the  additional  tuning.  A  closer 
examination  of  the  mean-squared  value  contribution  associ¬ 
ated  with  the  bending  states  (over  a  single  Monte  Carlo 
simulation)  revealed  that  an  increase  of  four  orders  of 
magnitude,  rather  than  two,  over  the  baseline  "R"  values  was 
appropriate.  Filter  tuning  for  the  "best"  values  resulted 
in  the  following,  which  is  the  result  of  adding  4|.  to  the 


values  identified  in  Equation  (4-22)  (the  value  of  4  was 
derived  from  taking  the  average  contribution  of  the  variance 
from  and  over  a  single  Monte  Carlo  run): 


R  = 


4.00436  0 

0  4.00598 


(6-3) 


The  corresponding  perforinance  plots  are  shown  in 
Figures  E-25  through  B-32  and  the  statistical  data  appears 
in  Table  6.12.  Notice  that  the  filter  exhibits  a  delayed 
response  to  the  maneuver  by  approximately  one-third  of  a 
second  as  illustrated  by  the  filter-computed  RMS  error 
curves  of  Figures  E-25  and  E-26.  In  addition,  these  plots 
show  the  effect  of  the  increased  measurement  noise  "masking” 
the  filter-computed  RMS  values,  implying  that  the  filter  is 
tuned  to  place  less  trust  in  the  measurements  and  to  rely 
iBorQ  on  the  internal  filter  dynamics  models.  The  Increase  in 
measurement  noise  also  implies  that  there  lies  a  greater 
uncertainty  associated  with  the  measurements,  which  Is  Indi¬ 
cated  In  Figures  E-27  through  E-30  by  the  larger  standard 
deviation  as  compared  to  the  previous  case  without  the  addi¬ 
tional  tuning,  as  depicted  In  Figures  E-18  through  E-24. 

This  relationship  can  also  be  seen  by  comparing  the  x-  and 
y-statistlcs  shown  In  Table  6.12  with  the  corresponding  data 
in  Table  6.11.  In  addition  to  the  Increased  standard  devia¬ 


tion,  the  y-channel  also  displays  a  bias  of  approximately 
-4.0  pixels,  beginning  at  the  onset  of  the  maneuver  and 


Table  6.12 


MMAP  with  Increased  Bending  and  Tuned 
OM/MMAP/T2/10-O/QB2/R+4 


Temporally  Averaged 

Time 

Interval 

Error  Parameter 
(mean  /  1  sigma) 

(0.5 

i 

2.01 

(3.5 

i 

5.01 

*err  <ti~) 

-.0660 

/ 

1.515 

.5191 

/ 

2.102 

yerr  (<:i~> 

.0200 

/ 

1.509 

-1.392 

/ 

1.581 

>«etr 

-.0506 

/ 

1.462 

.4398 

/ 

2.049 

Yerr 

.0136 

/ 

1.456 

-3.703 

/ 

1.514 

x-centetr  <^1*) 

-.1880 

/ 

.6455 

.4274 

/ 

.5502 

y-centgrr (ti^) 

-.0932 

/ 

.5986 

-4.015 

/ 

.4384 

y  peak-mean  error  (tj  )  =  -6.5  pixels 
y  peak-mean  error  (tj'*’)  »  -4.5  pixels 
recovery  time  -  .6  seconds 


continuing  throughout  the  duration  of  the  simulation.  This 
phenomenon  Is  also  exhibited  by  the  estimated  y-centroid 
error  as  Illustrated  In  Figure  B-32,  with  an  Initial  bias  of 
about  -3.0  pixels  and  Increasing  (negatively)  In  time. 

These  biases  in  the  y-channel  are  a  result  of  the  Increased 
smasureswnt  variances  In  the  filter-assumed  model.  Examina¬ 
tion  of  the  dominant  elemental  filters  for  this  simulation 
revealed  that  filter  12,  the  large  POV  filter,  remained 
dominant  throughout  after  the  tlem  of  the  maneuver 
Intlatlon.  This  Implies  that  the  Increased  measuresmnt 
variance  inhibited  the  abllitly  of  the  HMAP  algorithm  to 
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distinguish  between  the  "good*  and  "bad”  filters,  and  thus 
the  increased  bias  in  the  y-channel  is  due  to  filter  12*8 
poor  estiinate  of  y-channel  dynamics.  This  can  also  explain 
the  exceptionally  high  time-averaged  y-centroid  error  para¬ 
meter  value  shown  in  Table  6.12  between  3.5  and  5.0  seconds 
As  a  result  of  the  proposed  tuning  method.  Figures  B-24  and 
B-25  indicate  that  better  performance  is  attained  with  the 
higher  "R"  values;  however.  Figures  B-26  through  B-32 
indicate  the  converse  (especially  in  the  y-channel). 

As  a  further  investigation  into  this  concept  of  the 
unexpected  bias  contribution,  the  discrete-time  measurement 
noise  matrix  was  changed  to  reflect  the  values  of  1.5  where 
the  square  root  of  this  is  the  standard  deviation  contribu¬ 
tion  from  the  bending  states  over  a  ten  run  Honte  (2arlo 
simulation  (it  makes  more  sense  to  alter  the  measurement 
matrix  by  the  variance;  however,  the  intent  of  this  simula¬ 
tion  is  to  show  a  trend,  rather  than  a  physical  realiza¬ 
bility).  The  performance  plots  for  this  simulation  are 
shown  in  Figures  B-33  through  E-40.  Figures  E-33  and  B-34 
reflect  that  "better"  tracking  as  compared  to  the  case  with 
out  tuning  is  effected;  however,  a  similar  occurrence  to 
that  of  the  previous  tuned  case  is  evident  -  a  biasing 
effect  appears  after  the  maneuver  initiation  in  the  y- 
channel .  Again,  the  dominant  filter  after  the  maneuver  was 
filter  12,  and,  thus,  this  shows  that  too  large  of  a 
measurement  variance  was  being  utilized  in  order  for  the 
algorithm  to  adequately  distinguish  the  correct  model. 
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6.4.5.  Impact  of  Bend Inq/Vl brat  Ion  Phenomena .  As  has 
been  shown,  the  effect  of  modeling  bending  and  vibration  in 
the  truth  model  but  not  providing  the  filter  such  additional 
information  does  not  impede  the  Gauss-Markov  MHAP  tracking 
algorithm  from  accurately  pointing  and  tracking,  provided 
the  level  of  bending  and  vibration  is  on  the  order  of  that 
expected.  However,  if  the  level  of  bending  and  vibration  is 
greater  than  that  which  is  expected,  tracking  can  be  main¬ 
tained  but  accurate  pointing  for  the  purpose  of  focusing  a 
weapon  such  as  a  laser  is  severely  degraded.  The  two  tuning 
exercises  associated  with  the  increased  bending  level  have 
implied  that  the  result  of  increasing  the  measurement 

t 

variance  to  account  for  the  unmodeled  effects  impedes  the 
MMAF  algorithm's  ability  to  choose  the  appropriate  filter 
model  in  the  bank. 


6 . 5  Scenario  Sensitivity  Ana 1 vs i s 

This  section  addresses  the  sensitivity  of  the  tracking 
algorithm  to  variations  of  pixel  size,  target  type,  target 
range,  and  noise  characteristics.  The  simulations  performed 
for  this  analysis  do  not  include  the  effects  of  bending  and 
vibration  in  the  truth  model;  in  addition,  software 
modifications  have  been  introduced  at  this  point  to 
accommodate  the  description  of  two  additional  trajectories 
which  are  introduced  in  the  following  section.  These 
simulations  are  performed  using  the  Gauss-Markov,  HHAP 


tracking  algorithm  employing  a  trajectory  12  configuration. 
Table  6.13  shows  the  scenarios  considered  for  this  analysis. 

The  target  type  parameters,  or  correlation  time  pairs, 
of  8/3,  4/1.5,  4/. 8,  and  4/. 4  are  perceived  to  represent  the 
dynamics  of  a  long-range  missile,  a  bomber  aircraft,  a  high- 
performance  fighter  aircraft,  and  an  air-to-air  missile, 
respectively.  These  time  pairs  identify  the  correlation 
times  assumed  by  the  filter  for  both  the  benign  and  the 
dynamic  case;  for  the  benchmark,  the  benign  T=  4  seconds 
and  the  dynamic  t  =  1.5  seconds.  The  pixel  size  corresponds 
to  the  angular  f ield-of -view  of  a  single  pixel  and  the  range 
parameter  corresponds  to  the  target’s  initial  condition  at  a 
location  perpendicular  to  the  FLIR  FOV.  With  respect  to 
pixel  size,  those  simulations  that  vary  from  the  nominal 
(nominal  =  20  4rad  on  a  side)  necessitate  changing  the 
inertial  x-  and  y-initial  conditions  in  order  to  provide 
the  tracker  the  identical  initial  angular  orientation  of  the 
target  as  seen  by  the  FLIR  FOV  as  in  previous  studies.  This 
relationship  with  respect  to  elevation  and  azimuth  is 
depicted  in  Figure  2.6. 

The  noise  parameters,  Qr2'  respective 

continuous-time  noise  strengths  (nominal  values  determined 
by  computing  Equation  (5-1)  on  the  values  depicted  in  Table 
5.2)  for  the  dynamics  driving  noises  for  the  models  upon 
which  the  filters  are  based.  Recall  from  Equation  (5-1) 
that  "0"  is  expressed  as  a  function  of  the  acceleration 
variance  (or  mean  squared  value)  and  correlation  time  for 


Table  6.13 


Scenario  Sensitivity  Parameters 


t 

Target  Type 
(seconds ) 

Pixel  Size 
( radians ) 

Range 

(meters) 

Noise 

(Q) 

Maneuver 

1. 

8/3 

200  nrad 

2000000 

Qri 

2-G 

2. 

8/3 

200  nrad 

2000000 

Qri 

10-G 

3. 

4/1.5 

20  )lrad 

20000 

Qri 

10-G 

4. 

4/1.5 

20  4rad 

20000 

Qri 

20-G 

5. 

4/.  8 

20  4rad 

20000 

Qri 

10-G 

6. 

4/. 8 

20  Mrad 

20000 

Qri 

20-G 

7. 

4/.  4 

20  Mrad 

20000 

Qri 

10-G 

8. 

4/.  4 

20  Mrad 

20000 

Qri 

20-G 

9  . 

4/.  4 

2  mrad 

20000 

Qr2 

10-G 

10. 

4/.  4 

2  mrad 

20000 

Qr2 

20-G 

the  first-order  Gauss-Harkov  filters.  The  parameter  Qpi  in 
Table  6.12  implies  that,  for  those  simulations,  the  zero 
frequency  PSD  (and  low  frequency)  value  for  "Q”  must  be 
equivalent  to  that  of  the  benchmark  simulation  in  order  to 
retain  a  viable  performance  analysis  comparison,  or,  the 
mean  squared  value  of  acceleration  (where  the  energy  corres¬ 
ponds  to  the  area  under  the  PSD  curve  of  Figure  6.1)  remains 
the  same  between  the  scenarios  to  be  compared.  Therefore, 
for  those  simulations  where  the  target  type  is  different 
from  that  of  the  benchmark  simulation  (  *r =  4/1.5  seconds), 
the  corresponding  filter  variance  must  also  be  adjusted  to 
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retain  this  relationship.  Figure  6.1  illustrates  the  rela¬ 
tionship  between  "Q"  and  the  PSD  curve  £or  a  first-order 
Gauss-Markov  process.  The  noise  parameter  Q{^2  excep¬ 

tion  to  the  above  discussion  based  upon  the  intuitive  notion 
that,  as  the  pixel  size  decreases  in  magnitude  without  a 
corresponding  increase  in  range  magnitude,  then  a  similar 
proportional  increase  of  the  filter's  driving  noise  strength 
must  occur  in  order  to  compensate  for  the  "apparent”  changes 
in  target  dynamics  due  to  decreased  pixel  size.  Note  that 
adjusting  the  baseline  ”Q”  value,  Qri,  is  not  necessary  for 
the  first  scenario  listed  in  Table  6.13  as  compared  to  the 
nominal  since  the  pixel  size  not  only  decreased  by  two 
orders  of  magnitude  from  the  nominal  but  the  initial 
(nominal  «  20,000  meters)  range  coordinate  also  increased  by 


Figure  6.1  PSD  Plot  of  a  First-Order  Gauss-Markov  Process 


two  orders  of  magnitude.  This  corresponds  to  the  argument 
that,  as  the  pixel  size  decreases/increases  in  magnitude  and 
the  range  to  target  increases/decreases  by  a  proportional 
degree  of  magnitude,  the  driving  noise  strength  necessary  in 
the  filter  to  depict  the  same  magnitude  of  target  maneuvers 
is  unchanged.  In  addition,  this  differs  from  Tobin's  "Q-vs- 
Range"  function  (19:126)  described  by: 

Qr2  =  (r2/ri)^  Orl  (6-4) 

where : 

0^2.  =  strength  of  driving  noise  for  a  target  at  range  rl 
Qr2  =  strength  of  driving  noise  for  a  target  at  rabge  r2 

where  this  relationship  is  based  upon  targets  at  different 
ranges  while  maintaining  the  equivalent  pixel  size. 

Rather  than  discussing  all  ten  of  the  simulation 
scenarios  listed  in  Table  6.13,  scenarios  1,  4,  5,  8,  and  9 
will  be  discussed  since  these  simulations  encompass  the 
different  variations  of  parameter  values.  Scenario  12  is 
the  same  as  II  but  with  the  turn  maneuver  occurring  at  a 
higher  g-level.  Scenario  13  is  the  same  as  the  benchmark 
run  discussed  in  Section  6.4.1  and  scenarios  16,  17,  and  110 
are  similar  to  #5,  #8,  and  19,  respectively,  with  the 
exception  of  the  different  g-level  maneuver. 

6.5.1  Ranqe/Pixel  Size  Sensitivity.  Scenario  II 
considers  the  case  where  the  pixel  size  is  decreased  by  two 
orders  of  magnitude  and  the  range-to-target  is  increased  by 
two  orders  of  magnitude  (nominal  range  =  20000  meters). 


This  scenario  is  performed  for  a  target  represented  by  a 
"benign"  correlation  time  of  8  seconds  and  a  "dynamic" 
correlation  time  of  3  seconds  throughout  the  2-g  turn.  The 
performance  plots  for  this  simulation  are  found  in  Figures 
F-1  through  F-8  and  the  time  averaged  statistics  are 
compiled  in  Table  6.14.  The  statistics  for  this  simulation 
are  very  similar  to  those  of  the  benchmark  case  as  depicted 
in  Table  6.7  with  one  unique  exception:  the  y-peak-mean 
error  is  approximately  1/2  that  of  the  benchmark  run,  or  the 
filter  does  a  better  job  estimating  the  target's  position  at 
the  time  immediately  following  the  harsh  maneuver  when  the 

Table  6.14 

MHAF  Scenario  tl  Statistics 
GM/MMAF/T2/2-G/TS8, 3/Zo=2000000 


Temporally 

Averaged 

Time 

Interval 

Error  Parameter 

(mean  / 

1  sigma) 

[0.5 

/ 

2.0) 

(3.5 

,  5.0] 

Xgrr  ^ ^i 

) 

.1572 

/ 

.5273 

.1141 

/  .5398 

ye  r  r  ^  ^  i 

) 

-.0258 

/ 

.3543 

-.0426 

/  .3760 

Xgrr  ^ ^i 

) 

.1360 

/ 

.4230 

.1156 

/  .4331 

ye  r  r  ^  ^  i 

) 

-.0223 

/ 

.3114 

-.0259 

/  .3354 

x-centerr 

.1223 

/ 

.1035 

.0566 

/  .1345 

y-centerr 

.0049 

/ 

.0589 

-.0043 

/  .0817 

y  peak- 

mean  error 

(ti‘) 

=  -.9  pixels 

y  peak- 

mean  error 

(tl*) 

=  -.6  pixels 

recovery  time 

=  .  7  seconds 

s  %  s  s 


'  V^’ 
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target  is  positioned  at  great  distances  from  the  sensor. 
Figure  F-2  shows  this  phenomenon,  as  compared  to  the  bench¬ 
mark  case  in  Figure  E-2,  where  the  actual  RMS  error  changes 
drastically  for  the  time  between  2.0  and  3.0  seconds  and  the 
filter  does  a  good  job  of  following  the  changes.  Also  of 
particular  interest  are  the  different  character ist ics  asso¬ 
ciated  with  the  estimated  x-  and  y-centroid  positions  as 
shown  in  Figures  F-7  and  F-8  as  compared  to  Figures  E-7  and 
E-8  (note  the  difference  in  scales).  The  estimated  x-plus 
centroid  position  exhibits  a  positive  bias,  or  overestimate, 
until  the  time  when  the  maneuver  is  initiated,  shown  in 
Figure  F-7,  as  compared  to  the  relatively  stable  estimation 
shown  in  Figure  E-7.  On  the  contrary,  the  estimated  y-plus 
centroid  position  exhibits  a  much  more  stable  estimation, 
shown  in  Figure  F-8,  as  compared  to  that  of  Figure  E-8, 
which  exhibits  a  bias  for  time  greater  than  3.5  seconds. 

6.5.2  Hiqh-q  Maneuver  Sensitivity.  Scenario  14 
represents  a  target  of  correlation  time  pair  of  4/1.5 
seconds  pulling  a  20-g  maneuver  at  nominal  range  and  pixel 
size  and  utilizing  the  baseline  values  of  Performance 

plots  for  this  simulation  are  shown  in  Figures  P-9  through 
F-16  and  the  statistical  data  is  reflected  in  Table  6.15. 

As  had  been  demonstrated  in  a  previous  research  effort  by 
Tobin  (19:107,109),  the  filter  performed  rather  well  against 
this  scenario,  which  serves  as  the  benchmark  to  measure  the 
effectiveness  of  subsequent  20-g  scenario  simulations.  Note 
how  fast  the  filter  recovers  from  the  maneuver,  as  can  be 
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Temporally  Averaged 

Time 

Interval 

Error  Parameter 

(mean  /  1  sigma) 

10.5 

9 

2.01 

[3.5 

,  5.0] 

^rr  ^ ^i  ) 

-.0456 

/ 

.4261 

.5565 

/  .5519 

yerr (ti") 

o 

o 

• 

/ 

.3577 

.4192 

/  .6163 

^r r  ^ ^i  ) 

-.0354 

/ 

.3707 

.4294 

/  .4377 

Ye  r  r  ^  ^  i  ^ 

-.0137 

/ 

.  3173 

.5792 

/  .5144 

x-centgrr  ( ti"*^ ) 

-.0006 

/ 

.0948 

.2694 

/  .1697 

y-centerr  < ) 

-.0015 

/ 

.0536 

.6919 

/  .2719 

y  peak- 

mean  error 

(ti") 

=  -2.3  pixels 

y  peak- 

mean  error 

(t/) 

=  -1.1  pixels 

recovery  time 

»  . 5  seconds 

seen  by  comparing  Figure  F-2  to  E“2,  thus  implying  that  the 
filters  are  "better"  tuned  for  20-g  maneuvers  than  they  are 
for  10-g  maneuvers.  This  is  also  evident  In  Figures  F-6  and 
E-6  in  the  estimation  of  the  y-plus  position  and  Figures  F-8 
and  E-8  in  estimating  the  y-plus  centroid  position. 

6.5.3  Medium  Correlation  Time  Target .  Scenario  #5 
represents  a  target  of  correlation  time  pair  4/. 8  seconds  at 
nominal  pixel  size,  range,  and  noise  strength  pulling  a  10-G 
turn.  The  variances  associated  with  elemental  filters  12, 
•3,  14  x-dynamics,  and  15  y-dynamics  were  adjusted  to  5333, 
1066,  5033,  and  5033,  respectively.  This  adjustment  retains 
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the  equivalent  mean  squared  value  of  acceleration  that  was 
available  to  the  filters  in  the  benchmark  runs  so  that  a 
one-to-one  comparison  is  effected  (pictorially  illustrated 
in  Figure  6.1).  The  performance  plots  for  this  simulation 
appear  in  Figures  F-17  through  F-24  and  the  statistics  are 
shown  in  Table  6.16.  Evident  in  both  the  statistical  data 
and  the  performance  plots  are  characteristics  comparable  in 
performance  to  the  benchmark,  with  the  exception  that  the  x- 
channel  estimates  shown  in  Figures  F-19,  F-21,  and  F-23  (as 
compared  to  Figures  E-3,  E-5,  and  E-7,  also  note  difference 
in  scales)  for  both  position  and  centroid  position  exhibit  a 
time -increasing  bias  after  roughly  0.5  seconds  after  the 


Table  6.16 

MHAF  Scenario  15  Statistics 
OM/HMAF/T2/10-G/r=4, .8 


Temporally  Averaged 

Time 

Interval 

Error  Parameter 
(mean  /  1  sigma) 

[0.5 

9 

2.0) 

[3.5 

,  5 

.0) 

’‘err^^i  J 

-.0516 

/ 

.4243 

.  2955 

/  . 

5216 

Yerr  (ti~) 

-.0097 

/ 

.3574 

-  .0735 

/  . 

6532 

’^err^^i  ^ 

-.0374 

/ 

.  3699 

.  2476 

/  . 

4196 

Yerr 

-.0125 

/ 

.  3174 

.0912 

/  . 

5883 

x-centepj.  ( ti^ ) 

.0010 

/ 

.0978 

.1419 

/  . 

1456 

y-centgrr  <  ^ ) 

-.0005 

/ 

.0541 

.  3575 

/  . 

3041 

y  peak- 

mean  error 

(ti‘) 

=  -1.7  pixels 

y  peak- 

mean  error 

<ti*> 

=  -.7  pixels 

recovery  time  = 


.  5  seconds 


maneuver  has  been  initiated.  The  probable  cause  of  the  bias 
is  the  fact  that  the  filters  were  previously  tuned  for  a 
target  correlation  time  of  1.5  seconds,  thus,  the  strength 
of  the  driving  white  noise,  "0",  was  predicated  on  this  type 
of  correlation  time.  This  is  comparable  to  optimally  tuning 
the  Kalman  filters  to  track  bomber  type  aircraft  but 
performing  the  tracking  simulation  against  fighter  type 
aircraft.  One  means  to  handle  this  phenomenon  would  be  to 
provide  adaptive  estimation  of  parameters  such  as  the  filter 
"Q",  or  other  parameters  where  appropriate,  and  adjust  the 
parameters  to  react  to  the  changing  environment. 

6.5.4  Fast  Correlation  Time  Target .  Scenario  18  rep¬ 
resents  a  target  with  correlation  time  pair  of  4/. 4  seconds 
at  nominal  pixel  size,  range,  and  noise  strength  pulling  a 
20-G  turn.  Again,  the  variances  related  with  elemental 
filters  12,  #3,  14  x-dynamics,  and  15  y-dynamlcs  were 
adjusted  to  2667,  533,  2667,  and  2667,  respectively,  to 
account  for  the  matching  of  the  mean  squared  value  of  accel¬ 
eration  between  this  simulation  and  the  benchmark  (see 
Figure  6.1).  For  performance  comparison  purposes,  this 
simulation  is  measured  against  the  performance  exhibited  in 
the  Scenario  #4,  the  20-g  benchmark.  The  performance  plots 
for  this  simulation  appear  in  Figures  F-25  through  F-32  and 
the  statistical  data  is  shown  in  Table  6.17.  In  comparing 
the  performance  data.  Figures  F-27,  F-29  and  F-31,  as 
compared  to  Figures  F-11,  F-13,  and  F-15  (note  difference 


Table  6.17 


HHAP  Scenario  16  Statistics 
OM/MMAF/T2/20-G/t=4, . 4 


Temporally 

Averaged 

Time 

I nterval 

Error  Parameter 

( mean  / 

1  sigma) 

10.5 

2.01 

13.5 

,  5.0) 

^r r  ^ ^i 

) 

-  .0679 

/ 

.4213 

1.175 

/  .5201 

yerr 

) 

- .0089 

/ 

.3557 

-.8953 

/  .3924 

’‘err  ^  ^i 

) 

- .0467 

/ 

.3702 

.  8846 

/  .4199 

Verr 

) 

-.0119 

/ 

.3154 

-.6203 

/  .3609 

x-cent^pp (ti* ) 

.0013 

/ 

.0970 

.  5576 

/  .1525 

y-centerj.  (ti"^ ) 

.0004 

/ 

.0528 

.0665 

/  .1045 

y  peak- 

mean  error 

(ti') 

=  -3.0  pixels 

y  peak- 

mean  error 

=  -1.5  pixels 

recovery  time 

=  .  7  seconds 

in  scales),  the  bias  in  the  x-channel  as  seen  in  the 
scenario  15  is  prevalent.  In  fact,  this  bias  is  larger  than 
seen  previously,  primarily  due  to  the  greater  mismatch  of 
the  correlation  times  as  compared  to  the  benchmar)( .  Again, 
a  parameter  estimation  technique  would  seem  appropriate  to 
compensate  for  the  large  bias  errors. 

6.5.5  Pixel  Size/Noise  Strength  Relationship . 

Scenario  19  represents  a  target  of  correlation  time  pair  of 
4/. 4  seconds  at  nominal  range  and  pixel  size  of  two 
mllliradians  pulling  a  10-G  maneuver.  The  noise  strength 
for  this  simulation  was  decreased  by  two  orders  of  magnitude 
for  each  filter  to  correspond  to  the  same  increased  pixel 


%A 


size  from  nominal.  The  performance  plots  appear  in  Figures 
P-33  through  P-40  and  the  statistical  data  is  depicted  in 
Table  6.18.  Of  particular  interest  is  the  high  x-position 
estimation  error  for  time  less  than  1.2  seconds  as  compared 
to  the  benchmark  simulation.  Although  this  is  evident  in 
Figure  F-33,  it  is  more  clearly  visible  in  the  mean  error 
estimate  shown  in  Figures  F-35  and  F-37.  The  mean  value  of 
the  error  in  the  filter's  prediction  estimate  has  increased 
by  almost  two  orders  of  magnitude  compared  to  what  had  been 
experienced  in  the  benchmark  scenario.  This  is  felt  to  be 
due  to  the  fact  that,  although  the  pixel  size  had  Increased 


Table  6.18 

^ MMAF  Scenario  19  Statistics 

GM/MMAF/T2/10-G/r=4, .4/KPT=,002 


Temporally  Averaged 

Time 

Interval 

Error  Parameter 
(mean  /  1  sigma) 

10.5 

9 

2.01 

(3.5 

9 

5.0) 

r  <  ^i  J 

.2263 

/ 

.4067 

.0753 

/ 

.  3103 

yerr (ti”) 

-.0440 

/ 

.3241 

-.0932 

/ 

.  3087 

^e  r  t  ^  ^  i  ^ 

.1616 

/ 

.  2964 

.0756 

/ 

.2916 

ye  r  r  <  ^  i  ^ 

-.0422 

/ 

.2964 

-.0761 

/ 

.2889 

x-centefp(ti  ) 

.0070 

/ 

.1021 

.0132 

/ 

.0982 

y-centerr 

- .0156 

/ 

.0634 

-.0160 

/ 

.0632 

y  peak- 

mean  error 

(ti") 

=  -.3  pixels 

y  peak- 

mean  error 

(ti*) 

=  -.2  pixels 

recovery  time 


1.1  seconds 


by  two  orders  of  magnitude,  the  initial  position  and 
velocity  parameters  associated  with  the  target  truth  model 
had  not  changed;  therefore,  the  target  "appears"  to  be 
moving  less  harshly  at  the  onset  of  the  simulation  to  the 
larger  pixel  size  than  it  had  for  the  nominal  case.  Another 
interesting  facet  about  this  simulation  is  the  stablity 
associated  with  tracking  the  estimated  y-position.  The 
peak-mean  error  estimate  after  the  onset  of  the  turning 
maneuver  had  Improved  to  l/6th  of  that  estimated  by  the 
benchmark  run  given  In  Table  6.7  (as  expressed  In  pixels, 
even  though  the  physical  pixel  sizes  themselves  have  changed 
between  the  two  cases). 

Target  Trajectory  Sensitivity 

This  section  analyzes  the  sensitivity  of  the  tracking 
algorithm  against  two  target  trajectories,  15  and  16,  for 
which  the  algorithm  had  not  been  previously  tuned.  Unlike 
the  sensitivity  studies  in  the  last  section,  the  effects  of 
bendlng/vlbratlon  are  included  in  the  truth  model  for  these 
analyses.  The  trajectory  sensitivity  analyses  are  conducted 
with  the  GM  MMAF  tracking  algorithm  against  targets  pulling 
a  10-g  turn. 

6.6.1  Sensitivity  to  Trajectory  t5.  The  first 
trajectory  to  be  analyzed  is  trajectory  15.  As  shown  In 
Figure  2.5,  trajectory  15  Is  accomplished  by  rotating 
trajectory  12  by  45°,  thus  reorienting  the  truth  model 


trajectory  generation  In  space  rather  than  redefining  the 


orientation  of  the  Inertial  coordinate  system.  The  perfor¬ 
mance  plots  are  shown  In  Figures  G-1  through  G-8,  the 
statistical  data  appears  In  Table  6.1<>^  and  the  dominant 
filters  throughout  the  tracking  simulation  are  depicted  In 
Table  6.20.  The  statistical  data  shows  that  both  the  time 
averaged  siean  and  one  sigma  deviations  have  Increased  for 
all  time  as  compared  to  the  benchmark.  However,  a  direct 
comparison  of  x-channel  or  y-channel  results  for  the  two 
cases  doesn't  make  total  sense,  since  the  trajectory  charac¬ 
teristics  are  now  rotated  with  respect  to  the  tracker's  x- 
and  y-directions .  Also,  reviewing  the  time  histories  of  y- 
mlnus  position  and  y-plus  position  show  that  the  filter 
initially  (for  time  less  than  .8  seconds),  has  difficulty  in 
estimating  the  actual  y-position,  which  is  depicted  by  the 
mean  error  estimates  fo  Figures  G-3  and  G-5.  This  phenomena 
had  only  been  experienced  at  the  onset  of  the  turn  maneuver 
in  previous  scenario  simulations.  In  addition,  the  data 
enumerated  in  Table  6.20  shows  that  the  MHAF  algorithm 
does  not  consistently  identify  a  dominant  filter  for  any  one 
significant  stage  of  the  scenario.  This  is  due  to  the 
orientation  of  the  target  and  the  orientation  of  the 
velocity  vector  not  dominating  any  one  direction  with 
respect  to  the  f ield-of-vlew,  especially  prior  to  the 
maneuver  initiation,  but  equally  represented  in  both  the 
FLIR  X-  and  y-dlrectlons . 


Temporally  Averaged 

Time 

Interval 

Error  Parameter 

(mean  /  1  sigma) 

(0.5 

»  2 

.01 

(3.5 

5.0) 

*err^^i  ^ 

.0770 

/  . 

6824 

.1464 

/ 

.8269 

yerr^^i  ^ 

.0407 

/  . 

6911 

.1999 

/ 

.6552 

*e  r  r  ^  ^  1  ^ 

.0730 

/  . 

6302 

.0565 

/ 

.7491 

Ve  r  r  ^  ^  1  ) 

.0534 

/  . 

6106 

.2284 

/ 

.  5897 

x-centgrr  J 

-.0083 

/  . 

1048 

-.0748 

/ 

.1238 

y-centgrr 

.0420 

/  . 

0802 

.1353 

/ 

.1789 

y  peak- 

mean  error 

(ti*) 

=  -1.4  pixels 

y  peak- 

mean  error 

<ti*) 

=  -1.1  pixels 

recovery 

time 

=  .  6  seconds 

Frame 


Dominant  Filter(s) 


1 

-  14 

•  1 

15 

-  16 

•  1,  15 

17 

-  69 

11,  12,  13 

Remarks 

Acquisition  phase 

Attempts  to  follow  a 
y-dlxection  maneuver 

Intermittent  behavior 
between  y-direction 
and  10-g  maneuver 


70  -  150 


12,  14 


Attempts  to  follow  a 
x-direction  maneuver 


6.6. 2  Raroplnq  Acceleration  at  Turn  Initiation . 
Trajectory  16  is  of  the  same  configuration  as  trajectory  12 
with  the  exception  that  when  the  pull-up  maneuver  occurs 
(time  =  2.0  seconds)/  the  acceleration  increases  as  a  ramp 
function  and  attains  its  maximum  after  S  sample  periods  have 
expired/  rather  than  occurring  as  a  step  change.  This  type 
of  acceleration  profile  is  more  realistic/  less  harsh  than 
the  artificial  step  increase  in  acceleration.  Figure  6.2 
shows  the  initial  acceleration  profiles  at  commencement  of 
the  pull-up  maneuver  for  both  trajectory  12  and  16.  Notice 
that  the  greatest  difference  between  the  two  trajectories 
occurs  in  the  x-direction  acceleration  profile.  The 
trajectory  12  accelerations  for  the  x-  and  y-directions  are 
expressed  as: 


a^  =  -tc/V,^qCos  (u>t ) 


(6-3) 


ay  =  wVjjjggS in  (wt ) 


(6-4) 


where : 


‘^mag  =  magnitude  of  velocity  at  t  =  0 

CJ  =  turn-rate  of  the  target 

t  =  time,  measured  from  maneuver  initiation 

~  ^simulation  “  2  sec  (Note:  applicable  for  the 
maneuver  to  begin  at  tgimui^tion  =  2.0  sec) 


To  effect  a  ramping  acceleration,  Eguations  (6-3)  and  (6-4) 
are  premultiplied  by  a  dimensionless  ramp  function  of  the 
following  form  for  the  first  5  sample  periods  after  maneuver 


initiation: 


t 


where  the  6  Is  present  to  account  tor  the  ramping  effect  to 
attain  its  maximum  after  5  sample  periods  (each  sample  is 
l/30th  of  a  second).  Thus,  multiplying  Equations  (6-3)  and 
(6-4)  by  Equation  (6-5)  and  Integrating  the  outcome  produces 


•  1 
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TIME  2.0 
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b)  TRAJECTORY  #6 

UNITS  :  TIME 
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ACCELERATION  -  METERS/SECONDS 

Figure  6.2  Trajectory  12  and  16  Acceleration  Profiles 
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the  resultant  x-  and  y-direction  Inertial  velocity 
relationships  described  by: 

=  Vmag{6ltco8(/Jt )  -  ( l/(j)8in (a;t )  1  -  1)  (6-6) 

~  ^'^mag  ( (l/tL>)  (cos(6ut)  -  11  +  tsin(a;t)l  (6-7) 

and 

'^x  “  '^mag  @  t  =  0.0 
Vy  =  0.0  «lt  =  0.0 

Integrating  Equations  (6-6)  and  (6-7)  provide  the  following 
X-  and  y-direction  position  relationships: 

^t  ■  6V^g(  2/tL»^  (cos (Celt )  -  11  +  tAc>sin(Ct>t )  -  t)  +  xto  (6-8) 

yt  =  6Vn^g{  (  2/tc;^sin(Wt )  ]  -  t/coll  +  cos(tJt)n  +  yto  (6-9) 

where : 

Xfco  =  true  x-position  at  time  =  0.0 
yto  =  true  y-position  at  time  =  0.0 

Upon  implementing  these  relationships  into  the  GM  MMAF 
tracking  algorithm,  a  slight  enhancement  in  y-direction 
performance  is  expected  for  the  simulation  beyond  time  »  2.0 
seconds  as  compared  to  the  benchmark  data  with  bending 
included  appearing  in  Figures  E-9  through  E-16.  The  perfor¬ 
mance  plots  appear  in  Figures  G-9  through  G-16  and  the 
statistical  data  appears  in  Table  6.21.  Upon  comparing  the 
perfromance  plots  for  the  two  cases,  no  apparent  differences 
are  discernible.  As  the  data  in  Table  6.21  reveals,  the 
mean  errors  in  the  estimates  of  y-position  and  y-centroid 
error  have  improved  both  at  time  equal  to  t^  and  t^^  as 


compared  to  the  benchmark  including  bending  and  vibration  as 
shown  in  Table  6.9.  In  addition,  when  comparing  the  one 
sigma  values  for  all  time,  this  simulation  (for  which  the 
acceleration  begins  as  a  ramp)  shows  a  slight  improvement 
over  the  case  where  the  acceleration  acts  li)(e  a  step  in  the 
y-direction,  which  is  expected.  One  additional  note  is  that 
the  dominant  elemental  filters  for  this  simulation  (not 
shown)  are  basically  the  same  as  Table  6.10. 

Table  6.21 

Trajectory  16  Statistical  Data 
GM/MMAF/T6/10-G/QB1 


Temporally  Averaged 

Time 

Interval 

Error  Parameter 
(mean  /  1  sigma) 

[0.5 

/ 

2.0) 

[3.5 

,  5.0) 

Xq].].  ( t^  ) 

.0427 

/ 

.6309 

.1878 

/  .7882 

ye  r  r  ^  ^  i  ^ 

.0586 

/ 

.6439 

.0397 

/  .7693 

^rr  ^ ^i  ^ 

.0510 

/ 

.5869 

.1431 

/  .7275 

Yerr  <ti'^) 

.0439 

/ 

.6052 

.1803 

/  .6683 

x-centgp j ( t£  ^ ) 

.0025 

/ 

.1318 

.0468 

/  .1463 

y-centerr 

-.0091 

/ 

.0706 

.2814 

/  .2204 

y  peak- 

mean  error 

(ti") 

=  -2.2  pixels 

y  peak- 

mean  error 

(ti^, 

*  -.8  pixels 

recovery  time 

=  .6  seconds 

This  chapter  has  analyzed  several  aspects  pertaining  to 
airborne  tracking  systems  simulations.  First/  the 
potential  of  operating  at  a  50  Hz  sample  frequency  versus  a 
30  Hz  rate  was  investigated  for  both  the  Gauss-Harkov 
tracker  and  the  constant  turn-rate  model  tracker,  utilizing 
a  single  filter  configuration.  Secondly,  both  tracker 
models  were  evaluated  with  respect  to  adding  the  bending  and 
vibration  effects  of  a  large  space  structure  to  the  truth 
model  without  changing  the  basic  structure  of  the  filter, 
thus  representing  a  reduced-order  filter.  After  evaluating 
these  effects  against  the  single  filter  configuration,  the 
same  analysis  was  performed  with  the  Gauss-Markov  MMAF 
tracking  algorithm.  In  addition,  the  magnitude  of  the 
bending  was  increased  and  filter  tuning  was  performed  to 
compensate  for  this  unmodeled  effect.  After  the  bending  and 
vibration  analysis,  a  sensitivity  study  was  performed  in 
order  to  evaluate  the  effects  of  changing  pixel  size,  target 
type,  or  correlation  time,  driving  noise  strength,  and 
sensor-to-target  range.  This  scenario  sensitivity  study  was 
then  followed  by  a  sensitivity  study  of  two  target  trajectory 
profiles  for  which  the  filters  had  not  been  tuned.  The 
conclusions  drawn  from  these  analyses  are  summarized  in  the 
following  chapter. 


VII  . 


Assessments  and  Recommendations 


7 . 1  Introduction 

This  chapter  assesses  the  impact  of  the  simulation 
analyses  discussed  in  the  previous  chapter  and  suggests 
research  areas  requiring  further  study.  Section  7.2 
examines  the  simulations  discussed  in  Chapter  VI  and 
evaluates  each  for  their  effect  on  the  PLIR  tracking 
algorithms.  Section  7.3  addresses  the  preliminary  studies 
performed  on  investigating  a  rotating  rectangular  field- 
of-view  concept  and  Section  7.4  enumerates  recommendations 
for  further  study. 

7 . 2  Assessments 

7.2.1  Practicality  of  the  50  Hz  Sample  Frequency.  The 
potential  of  implementing  the  tracking  algorithm  at  a  50  HZ 
rate  rather  than  at  a  30  Hz  rate  was  investigated.  The 
impact  of  processing  at  the  50  Hz  rate  results  in  slight 
improvements  in  estimation:  both  the  mean  error  and  the 
variance  about  the  mean  error  in  both  the  x-  and  y- 
directions  decrease  by  a  small  amount  (average  6\  decrease). 
However,  the  computer  loading  penalty  paid  as  a  result  of 
the  more  frequent  sampling  greatly  outweighs  the  potential 
benefits  obtainable.  As  an  example,  the  computer  processing 
time  associated  with  a  ten  Monte  Carlo  simulation  utilizing 
the  single  filter  Gauss-Harkov  tracking  model  results  in  a 
73\  increase  for  the  50  HZ  version  over  that  which  had  been 
experienced  with  the  30  Hz  version.  One  would  expect  a  6T% 
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Increase  since  50/30  =  1.67;  the  remainder  of  the  Increase 
is  due,  perhaps,  to  the  overhead  associated  with  the 
statistical  evaluations  when  performed  on  the  CYBER.  As  a 
consequence,  the  30  Hz  sampling  frequency  Is  sufficient  for 
the  tracking  algorithm  to  provide  good  performance  and  at 
the  same  time  minimize  computer  processing  time. 

7.2.2  Impact  of  Hardware  Bendlng/Vlbratlon .  The  issue 
related  to  the  repurcussions  due  to  the  bendlng/vlbrat ion 
phenomena  addresses  two  concerns.  The  first  concern 
considers  the  effect  on  the  filter  algorithm  as  the  expected 
degree  of  bending  and  vibration  of  the  space  structure  (see 
Appendix  A)  is  considered  In  the  truth  model  while  the 
filter  Is  not  provided  Information  regarding  this  phenomena. 
This  is  a  robustness  concern. 

Sections  6.3.1  and  6.3.2  discuss  this  Impact  with 
respect  to  single  filter  tracking  algorithms  based  on  Gauss- 
Harkov  acceleration  and  constant  turn  rate  models, 
respectively,  and  Section  6.4.2  discusses  this  issue  with 
respect  to  Gauss-Markov  MMAF  algorithm.  In  all  three 
instances,  the  inclusion  of  the  bending  and  vibration 
phenomena  displays  little  or  no  effect  to  the  mean  position 
errors  as  compared  to  the  same  scenario  without  the 
bend ing/vlbrat ion  states.  However,  the  error  standard 
deviations  associated  with  the  error  reflected  Increases  on 
the  order  of  0.2  to  0.3  pixels.  Only  in  the  case  of  the  y- 
posltlon  estimate  mean  error  and  one  standard  deviation  with 


the  CTR  filter  did  the  total  RHS  error  exceed  one  pixel, 
thus  distorting  the  tracking  ability  of  the  algorithm.  Both 
Gauss-Harkov  algorithms,  the  single  filter  and  the  HMAP, 
maintained  favorable  tracking  ability  with  the  expected 
vibration  and  bending  levels  included  in  the  truth  model. 

The  second  concern  addresses  the  effect  of  increasing 
the  vibration/bending  phenomena  by  an  order  of  magnitude; 
this  was  presented  in  Section  6.4.3  without  additional 
tuning.  This  simulation  showed  that,  although  the  filter's 
mean  errors  did  not  change  significantly,  the  one  sigma 
values  increased  by  a  factor  of  four  or  five,  to  between  1.5 
to  almost  2.0  pixels.  The  desired  tracking  accuracy  to  less 
than  a  pixel  in  resolution  is  severely  impeded  by  the  large 
standard  deviation  values.  Thus,  to  maintain  adequate 
resolution  and  tracking  performance  capability,  the  degree 
of  bending  and  vibration  due  to  the  flexure  of  the  structure 
should  not  exceed  that  discussed  in  Appendix  A,  provided  no 
filter  remodeling  has  compensated  for  this  effect.  If 
retuning  is  allowed,  as  was  presented  in  Section  6.4.4,  it 
was  shown  that  the  MMAF's  ability  to  select  the  appropriate 
filter  model  for  accurate  tracking  and  pointing  was  severely 
impeded  by  the  measurement  variance  compensation. 

7.2.3  Filter  Sensitivity  to  Varying  Parameters .  This 
section  summarizes  the  relationships  observed  for  varying 
selected  parameters  in  the  tracking  simulation  as  discussed 
in  Section  6.5.  The  first  scenario,  scenario  11,  shows  that 
increasing  target  range  by  two  orders  of  magnitude  from 


nominal,  where  the  nominal  is  20,000  meters,  results  in 
comparable  performance  (mean  t  1  sigma  values)  to  the 
benchmark  scenario  when  the  pixel  size  is  decreased  by  two 
orders  of  magnitude  from  nominal  (20  )irad  on  a  side).  This 
concept  provides  a  viable  comparison  to  any  other  tracking 
algorithm  where  either  the  range  or  the  pixel  size  is 
different  from  that  which  is  employed  in  this  algorithm. 

The  second  scenario,  scenario  14,  consists  of  the  same 
parameters  as  the  benchmark  scenario  with  the  exception  that 
the  target  pulls  a  20-g  turn  instead  of  a  10-g  turn.  This 
scenario  is  performed  with  the  intent  of  establishing  the 
benchmark  of  performance  for  a  subsequent  20-g  simulation 
and  demonstrating  the  capability  of  the  tracker  to  maintain 
good  performance  against  a  vehicle  pulling  a  20-g  turn. 

This  performance  data  proved  comparable  to  previous  efforts. 

The  third  and  fourth  scenarios,  scenarios  15  and  17, 
reflect  different  realizations  of  filter-assumed  correlation 
times  than  had  previously  been  assumed.  Both  showed  that, 
when  tracking  a  target  represented  by  a  correlation  time 
other  than  that  for  which  the  filter  was  tuned,  tracking 
performance  was  not  severely  impacted.  However,  both  cases 
exhibited  a  time-increasing  bias  in  the  x-channel  commencing 


at  the  time  of  the  maneuver  initiation.  Also,  the  bias 
tends  to  worsen  as  the  correlation  time  gets  faster,  i.e., 
represented  by  a  smaller  t.  Also  recall  that  the  mean 
squared  value  of  acceleration  with  respect  to  the  driving 
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white  noise  remained  equivalent  to  that  in  the  benchmark 
case  where  the  slower  correlation  time  had  been  assumed. 

This  gives  rise  to  the  conclusion  that,  given  a  bank  of 
filters  based  upon  a  particular  set  of  assumed  correlation 
time  parameters,  if  a  time-increasing  bias  in  the  mean  error 
is  prevalent,  a  possible  way  to  counteract  this  effect  would 
be  to  increase  the  strength  of  the  filter  driving  white 
noise,  possibly  via  an  adaptive  parameter  estimation  method. 

The  fifth  scenario,  scenario  19,  demonstrated  the  rela¬ 
tionship  between  pixel  size  and  filter  "Q",  given  that  all 
other  parameters  don't  change.  This  demonstration  showed 
that,  as  the  pixel  size  decreases,  the  filter  white  noise 
strength  should  decrease  by  a  corresponding  amount  (both 
decrease  or  both  increase)  if  the  same  basic  estimation 
characteristics  are  to  be  maintained. 

7.2.4  Target  Trajectory  Results .  The  tracking 
algorithm's  performance  against  two  target  trajectories 
not  previously  considered  showed  that  the  algorithm 
adequately  tracked  the  targets.  However,  in  both  cci.jes,  the 
RHS  value  of  the  error  statistics  approached  the  total  of 
one  pixel,  the  bound  upon  which  good  tracking  is  predicated. 

7.3.  Investigations  of  Rotating  Rectangular  Field-of-View. 

Preliminary  investigations  as  to  the  feasibility  of 


implementing  a  rotating  rectangular  f ield-of -view  filter 
have  been  conducted.  The  intent  of  these  investigations  was 
to  determine  a  method  by  which  a  rectangular  POV  filter 


could  be  rotated  so  as  to  align  the  "elongated"  side  with 
the  best  estimate  of  the  acceleration  vector.  The  reason 
behind  desiring  a  rotating  filter  is  the  possibility  of  the 
tracking  algorithm  to  perform  its  function  with  four  filters 
rather  than  five  (or  even  more  if  one  were  to  allow  a  larger 
f ield-of-view  dimension  to  exist  in  directions  other  than 
pure  azimuth  and  elevation),  as  had  been  done  previously. 

The  first  instinct  in  pursuing  this  task  was  to  utilize  the 
estimated  acceleration  vector,  as  determined  from  the 
multiple  model  adaptive  estimator  subroutine  embedded  in  the 
software,  for  orienting  the  longer  side  of  a  rectangular 
f ield-of-view.  Analysis  of  the  data  provided  from  a  single 
Monte  Carlo  simulation  for  the  Gauss-Harkov  tracking 
algorithm  against  a  trajectory  12,  10-g  maneuver  revealed 
that  utilizing  the  estimated  acceleration  was  not  a  good 
choice.  The  reason  for  this  is  that  the  estimated  accelera¬ 
tion  vector  is  a  noisy  estimate  that  can  change  quite 
rapidly,  and  did  in  fact  change  sign  (direction)  from  frame 
to  frame  at  one  and  one-half  seconds  into  the  scenario. 

This  is  where  the  estimated  acceleration  vector  tended  to 
oscillate  between  the  first  and  second  quadrant  (of  a 
Cartesian  coordinate  system)  for  a  period  of  17  frames.  In 
addition,  from  a  f rame-to-f rame  evaluation  of  the  respective 
angular  separation  of  sequential  estimated  acceleration 
vectors,  rotations  of  30  )irad  to  135.31  mrad  would  be 
required  in  order  for  the  POV  to  Rkaintain  alignment  of  its 
"elongated"  side  with  this  vector.  The  inconsistency 


associated  with  estimating  the  acceleration  vector  and  the 
corresponding  large  rotation  angle  (135  nrad  0.75  deg) 
necessary  to  rotate  an  array  of  picture  elements  compelled 
further  consideration  of  the  rotational  concept  with  respect 
to  following  the  estimated  velocity  vector. 

Initial  evaluation  of  the  estimated  velocity  vector  as 
a  candidate  by  which  to  command  rotations  showed  much  more 
stable  characteristics  than  did  the  estimated  acceleration 
vector.  Oscillations  from  frame  to  frame  were  not  apparent 
and  angular  difference  from  frame  to  frame  ranged  from  as 
small  as  10  )irad  to  as  large  as  7.38  mrad.  Rather  than 
attempt  to  align  the  "elongated”  axis  with  the  estimated 
velocity  vector  (as  was  the  intent  if  estimated  acceleration 
were  used),  aligning  the  "short"  axis  of  the  rectangular  POV 
with  the  velocity  vector  direction  is  considered  since  this 
maintains  the  acceleration  vector  along  the  "elongated”  side 
(assuming  that  acceleration  is  predominantly  orthogonal  to 
velocity) . 

Several  rotational  methods  were  attempted  using 
different  combinations  of  the  filter  states,  the  truth  model 
states,  and  the  hotspot  intensity  pea)(s,  i.e.  rotate  filter 
states,  truth  model  states,  and  hotspot  intensity  pea)cs  all 
in  same  direction  or  rotate  filter  states  and  truth  model 
states  in  one  direction  and  rotate  hotspot  intensity  pea)cs 
in  the  opposite  direction.  All  evaluations  of  these 
rotations  provided  similar  results.  The  traclcing  algorithm 


perforned  well  initially  and  then  tended  to  "lose  track”  for 
the  majority  of  the  simulations.  This  problem  stemmed  from 
the  fact  that,  although  the  states  from  the  current  frame 
are  rotated,  the  rotated  image  array  generated  from  the 
current  frame  is  correlated  with  the  ”smoothed”  image 
(template)  from  the  previous  frame  which  is  not  rotated  with 
respect  to  the  current  frame.  This  correlation  of  the 
present  data  with  the  previous  image  (smoothed)  is  the 
mechanism  by  which  the  filter  measurements  are  extracted. 

In  effect,  what  has  been  shown  is  that,  if  the  currently 
generated  rotated  image  is  correlated  with  an  unrotated 
image  (or  otherwise  interpreted  as  delaying  the  rotation  by 
a  sample  period),  the  performance  capability  of  the  Gauss- 
Harkov  tracking  algorithm  degrades  significantly.  This  is, 
in  effect,  a  suboptimal  Implementation  of  the  rotating  POV 
method.  One  means  of  compensating  for  the  time  lag  would  be 
to  rotate  the  template  (previously  generated  image)  by  the 
same  amount,  and  at  the  same  time,  as  the  filter  estimates. 
However,  if  the  PLIR  image  plane  is  physically  rotated  by  an 
angle  0,  the  image  in  that  plane  appears  to  rotate  by  -0; 
therefore,  the  simulation  of  the  true  image  rotation  by  0 
corresponds  to  the  tracker's  x  and  y  axes  rotated  with 
respect  to  the  original  orientation  and  the  reconstructed 
image  requiring  a  rotation  by  -0. 

K  method  of  implementing  a  rotation  to  a  spatial  image 
array  can  be  found  in  (6).  This  method  transforms  an  N  X  N 
x-y  spatial  array  to  a  2N  X  N/2  R-theta  plane  array  as 


depicted  by  Figure  7.1.  Each  row  In  the  R-Theta  plane 
corresponds  to  an  angular  displacement  in  the  x-y  plane  from 
the  center  of  the  plane;  each  row  (or  each  elemental 
increment  in  the  horizontal  direction)  equates  to  an 
increment  of  x/N  radians  if  an  image  is  considered  to  cover 
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360°  In  the  x-y  plane  (synnetry  can  be  exploited  if  the  x-y 
plane  lies  in  the  Fourier  domain,  which  is  presently  not 
considered,  and,  in  this  case,  an  image  can  be  considered  to 
cover  180°  in  the  mapped  space  (x-y  plane)  and  accompanied 
by  its  complex  conjugate). 

Rotations  are  simply  performed  in  the  R-Theta  plane  by 
shifting  rows  up  or  down  by  the  number  of  incremental 
angles  that  correspond  to  the  desired  angle  of  rotation. 

One  would  expect  not  to  be  able  to  rotate  by  exactly  a 
desired  amount,  but  to  the  "closest"  incremental  rotation 
possible  based  upon  the  number  of  "N"  discrete  quantities. 
Since  the  image  array  in  the  GH  tracicing  algorithm  is 
represented  by  a  24  X  24  array,  each  row,  if  this  method 
were  to  be  employed,  would  correspond  to  an  increment  of 
r/48  radians.  To  implement  such  a  rotational  algorithm  for 
this  tracl^er,  preliminary  investigations  have  disclosed  that 
angle  increments  in  Theta  on  the  order  of  1  to  3  mrads  would 
be  necessary  for  the  filters  to  be  able  to  maintain  "lock" 
on  a  f rame-to-f rame  basis.  One  topic  discussed  in  the 
following  section  is  a  proposed  implementation  technique 
that  may  warrant  further  investigation. 


7 . 4  Recommendations  for  Further  Study 

The  following  recomendations  are  suggested  for  further 
study  in  enhancing  the  FLIR  tracking  system: 

1.  Further  investigation  of  the  effects  of  the 


bending/vibration  phenomenon  may  be  necessary.  If  the  level 


of  bending/vibration  is  actually  greater  than  that  which  is 
expected  (see  Appendix  A),  additional  filter  tuning  should 
be  performed  (taking  into  consideration  the  bias  effect)  or 
the  possibility  of  modeling  the  bending  effects  in  the 
filter  model  should  be  conducted.  Representing  the  bending 
effects  in  the  filter  model  via  reduced  order  modeling 
should  be  considered. 

2.  Implementing  the  rotating  rectangular  f ield-of-view 
should  be  Investigated.  One  possible  Bieans  of  rotating  the 
image  template  would  be  to  reconstruct  the  image  array 
consisting  of  a  24  X  24  intensity  data  array  into  a  96  X  96 
intensity  data  array  (effectively  transforming  each  elesient 
in  the  24  X  24  array  into  4  elements  in  the  96  X  96  array) 
and  perform  the  R-Theta  transformation  previously  described. 
This  redefined  image  could  then  be  rotated  in  increments  of 
11/192,  or  3.27  mrad .  Although  the  96  X  96  array  definition 
provides  a  good  starting  point  for  possible  implementation, 
it  is  certainly  not  the  only  possible  choice.  The  image 
could  then  be  inverse  transformed  back  to  the  x-y  spatial 
domain  and  redefined  to  a  24  X  24  intensity  data  array  to  be 
correlated  with  the  currently  constructed  image. 

3.  The  problem  of  initial  target  acquisition  should  be 
investigated.  Since  accurate  velocity  and  acceleration 
estisetes  are  not  available  a  priori  (as  was  assumed  in  this 
effort),  they  should  be  determined  on-line.  One  potential 
solution  would  be  to  estimate  velocity  and  acceleration  by 


the  method:  =  4x/4t  and  a,(  =  Av'^/At,  where  the  x-centroid 

and  y-centrold  locations  at  adjacent  sample  times  could  be 
used  to  determine  ^x  and  ^y.  Also,  At  should  be  made  as 
small  as  computation  time  will  allow  and  still  minimize  the 
time-to-acquisition.  Another  potential  solution  would  be  to 
employ  the  inverse  covariance  Kalokan  filter  propagation  and 
update  forms  described  in  Maybeck  (9)  to  estimate  initial 
position,  velocity,  and  acceleration  variables  on  the  basis 
of  a  batch  of  N  measurements  (no  a  priori  information 
assumed,  i.e.,  ^  =  0.),  and  then  switch  back  to  the 

conventional  form  once  P  ^  becomes  nonsingular. 

4.  Still  requiring  investigation  for  this  tracking 
algorithm  is  an  indepth  target/decoy  sensitivity  analysis. 
Sensitivity  to  the  following  parameters  should  be  pursued: 
(1)  decoy  intensity,  shape,  and  size,  (2)  separation  from 
the  target  during  decoy  ignition,  and  (3)  tracking  time 
prior  to  decoy  release.  Image  and/or  filter  residual 
characteristics  should  be  analyzed  in  order  to  suggest  ways 
in  which  decoys  may  be  distinguished  from  targets. 


The  mechanical  bendlng/vibrat Ion  model  used  in  this 
study  was  based  on  a  report  prepared  for  APtfL  by  the  staff 
of  R  &  D  Associates  and  of  Cambridge  Research  (15).  This 
model  is  by  no  means  the  only  possible  structural  configura¬ 
tion  for  a  space-based  structure;  however,  it  does  provide 
the  necessary  mechanical  considerations  felt  necessary  to 
describe  the  effects  of  bending  and  vibration  at  several 
locations  about  a  typical  large  structure.  The  predominant 
sources  of  random  excitation  in  a  space  based  structure  are 
caused  by  (15:4-1): 

1.  Coolant  flow  through  the  mirrors,  the  resonator 
power  management  equipment  and  the  associated 
supply  systems; 

2.  Fluctuating  combustion  chamber  pressures; 

3.  Fluctuating  pressure  generated  by  the  laser  exhaust 
flows; 

4.  Conditioning  and  transfer  of  laser  reactants;  and 

5.  Steady  state  operation  of  the  control  moment  gyros. 

The  power  spectral  density  (PSD)  curve  of  Figure  2.2,  which 
is  shown  on  the  following  page  as  Figure  A.l,  was  determined 
to  be  the  predominant  contributor  to  the  bending/vibration 
effects  on  a  tracker's  line  of  sight,  with  the  response  in 
both  the  X-  and  y-  directions  being  similar.  Although  this 
PSD  response  represents  the  disturbance  seen  at  the 
secondary  mirror  LOS,  a  similar  response  with  approximately 
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a  60  dB  isolation  between  the  two  is  expected  at  the 
tracker's  LOS  (5).  This  research  effort  approximates  the 
response  of  the  system  due  to  the  PSD  curve  of  Figure  A.l  to 
be  represented  by  a  second  order  Harkov  process  described  by 
the  transfer  function: 


G(s) 


2  2 
s  2^ti>|,s  4- 


(A-1) 


where : 


bending/vibration  mode  natural  frequency', 
3.14159264  rad. /sec. 


.5  Hz. 


=  gain  adjusted  to  obtain  the  desired  root  mean 
square  (RMS)  bending/vibration  output,  c'bv* 
where  oi^v  is  defined  later  in  Equation  (A-23) 

=  damping  ratio,  =  .15  (17:359) 


The  roots  of  the  character istic  equation  can  be  determined 
by  solving  the  quadratic  equation  formed  by  setting  the 
denominator  of  Equation  (A-1)  to  zero,  resulting  in: 


‘1,2  =  <^b  ^  i‘^b  =  -.47124  +  J3. 10605 


(A-2) 


The  corresponding  state  space  representation  in  standard 
observable  form  for  both  directions  would  appear  as: 


Fll  Pl2 

F21  F22 


Fll  ^12 
F21  F22 


Xbi  0  0 

Xb2  i5b  0 

+ 

Ybl  0  0 

yb2  0  Gb 


(A-3 


and : 


’  V’  N.'  O  ^  ’ 


where : 


Fii  =  0 

Fi2  =  1 

F21  =  “CL>n^  '  -9.647547  rad  .  ^/sec  .  ^ 

F22  =  ~  "*94248  rad. /sec. 

Gb  =  =  9.647547Kb 

Equations  (A-3)  and  (A-4)  are  of  the  form: 

* 

X  =  Fx  +  Gw  (i 

and 

2.  “  Hx  (. 


The  state  transition  matrix,  can  be  derived  by 
employing  the  following  relationship: 


tb^t) 


L  ^UsI-EbJ 


At 


where : 

L  =  the  Laplace  operator 
I  =  the  identity  matrix 
At  =  the  sample  period,  t^^i  -  t^ 


( 


Performing  the  operation  defined  by  Equation  (A-7)  on  £5 
results  in  the  following  form  for  the  bending  phenomena 


state  transition  matrix: 


♦bl  ♦b2  0  0 

♦b3  *54  0  0 

0  0  ♦bl  ^52 

0  0  ♦bS  ♦b4 


(A-8) 


where 


=  exp(-Ob4t )  (co8(6Jb4t)  +  (Ob/‘*^b)8ln(tJ54t )  1 

=  .9945787985 

#b2  =  exp(-ab4t)  (  (l/tJb)sln(tJb/3t)  1 
=  .03275523095 

♦b3  “  exp( -OTb/it )  { -  ( 1  (Cb/Wb)  J  ^ 

=  -.1040812213 

♦b4  =  exp(-Ob^t)  Icos(ajbi4t)  -  (CTb/‘*^b)®i"<‘^b'^^  )  J 
=  .9637076486 

A  discrete-time  model  of  the  propagation  of  Xb  is  given  by; 


21b<ti  +  i)  »  it^b(*^i>  ♦  Kbd 


(A-9) 


where : 


Wbd  =  /  ♦b*'^l  +  l/ >Sb(7')ab<'^>‘*'^ 


snd  Wb/j  has  the  following  statistics: 


(A-10) 


EIwbd(<^i)J  *  0 
ElsaxJ  <^i  )Sfbd^<^j  )  j  =  Qbd^ij 


(A-11) 


(A-12) 


ti  +  1 


tb(ti+i/T)ab<'r)abC^)ab^(T)tb’^(ti+i,T)dT  (a-i3) 


Qbdl  Qbd2  0  0 

Qbd3  ttbd4  0  0 

0  0  Qbdl  Qbd2 

0  0  Qbd3  Qbd4 


where : 


Qbdl  = 


Qbd2 


Qbd4 


(A-14) 


ti  +  1 


Qb^*b2^<i'r  =  -3254294713  Kb^ 


■i  +  1 


Qb^*b2*b4d'r  =  Qbd3  =  -01038495095 


1+1 


Gb^*b4^<J'r  •  .3055407229 


ti 


Note  that  the  elements  of  Q^xl  ^^ve  been  expressed  in  terms 
of  the  power  spectral  density  gain  parameter,  K^.  Although 
current  technology  assessments  on  large  space  structure 
mechanical  disturbances  indicate  radial  displacement  levels 
consistent  with  those  of  Figure  A.l,  this  parasmter  can 
easily  be  varied  to  analyze  the  level  of  disturbance  that 
would  degrade  filter  performance. 

For  the  simulation  implementation,  the  discrete-time 
propagation  equation  was  Implemented  in  the  following  form: 

2£b(ti+i)  -  ♦bJLb<*^i>  *  e/ttbd  Wn<ti)  (A-15) 


where : 


y  Qbd  =  Cholesky  square  root  natrlx  of 

>  vector  of  independent,  unit-variance, 
gaussian  noises 

In  order  to  deteraine  the  value  of  Ki^,  the  form  of  the 
covariance  matrix  must  be  knovm.  If  uncorrelatedness 
between  the  x-  and  y-  directions  are  assumed,  then  the 
covariance  matrix  appears  as: 

Eb  = 


Pbl 

Pb2 

0 

0 

Pb3 

Pb4 

0 

0 

0 

0 

Pbl 

Pb2 

0 

0 

Pb3 

Pb4 

Now,  if  the  tisie  derivative  of  the  covariance  matrix  is  set 
equal  to  zero,  the  steady  state  value  of  K5  can  be  derived 
by  solving  the  following  relationship: 

Eb  =  EbEb  +  EbEb’*'  +  SbilbSb’'’  *  0  (A-17) 

which  yields  a  4  X  4  block  diagonal  matrix  where  each  2X2 
block  is  of  the  form: 


0 


2Pb2 

Pb4  *  f’bsPbl  ♦  fb4Pb3 


Pb4  +  ^ba^bl  +  Pb4**b2 
2Pb3Pb2  ♦  2Fb4Pb4  ♦  ^b^ 


<A-18) 


This  implies  that: 


Two  equations  with  three  unknovm  values  (Pbl/  ^b2' 


remain : 


Pb4  +  Pb3Pbl  =  0 


2Fb4Pb4  ♦  Gb  =  0 


(A-20) 


(A-21) 


Another  equation  relating  one  or  more  of  the  unknowns,  Gb, 
Pbl/  Fb4/  Is  necessary  In  order  to  accommodate  a  solution 
to  the  above  relationship.  The  variance  of  the  output 
provides  such  a  relationship: 


Pz2  =  Etzz  > 

=  E(H)<bllbVl 


ElXbi  1  0 


0  BlYbl  J 


(A-22) 


Making  use  of  the  fact  that 


EtXbi^l  =  ECybi^J  =  Pbl  = 


(A-23) 


and  Obv  equal  to  the  area  under  the  PSD  curve  which  is 

approximated  by  the  following  equation  (a  geometrical 
approximation  to  the  area  of  Figure  A.l): 


( .5)(Kb^)  +  (l/2)( .2)(14Kb^) 


where : 


+  (l/2)(700)(15Kb  ) 
5251. 9(Kb^)  rad.^/sec.^ 


=  PSD  magnitude  at  zero  frequency 


(A-24) 


The  first  term  in  this  approximation  is  represented  by  the 
area  in  Figure  A.l  shaded  by  horizontal  lines,  the  second 
term  by  the  area  shaded  by  vertical  lines,  and  the  third 
term  by  the  diagonal  lines  (Note:  levels  below  lE-19  are 
ignored).  Equation  (A-24)  is  used  to  solve  Equation  (A-20) 
for  P54  in  terms  of  the  zero  frequency  PSD  value,  K5: 

Pb4  “  t^n^fl^bv^  “  5.103E+O4(K5^)  rad.^/sec.^  (A-25) 


The  derivations  thus  far  have  defined  all  of  the  parameters 
necessary  to  perform  a  time  propagation  of  the  states  of  the 
bending/vibration  model.  In  order  to  provide  measurement 
update  relationships  compatible  with  the  tracking  algorithm, 
the  measurement  variables  should  be  expressed  in  terms  of 
pixels.  Since  the  propagation  variables  are  expressed  in 
terms  of  radians,  a  simple  conversion  from  radians  to  pixels 
is  necessitated,  with  the  resulting  output  relationship 
appearing  as: 


Ifb 


0  0 


0 


Sb 


I 


s 


1 


(A-26) 


Get  the  initial  hotspot  offsets  from 
the  target  center  of  mass 


Begin  tracking  the  target  for  the 
number  of  frames  read  in  from 
the  input  data  file 


Define  the  Gaussian  peak  locations 
for  target  image  simulation  based 
on  centroid  position 


Get  the  measurement  noise  array  and 
define  a  measurement  noise  array  for 
each  filter’s  f leld-of-view  size 


Calculate  the  measurement  noise 
for  each  filter 


Determine  the  measurement  array 
for  each  filter 


Pad  the  measurement  arrays  with  zeros 
for  those  elements  greater  than  the 
filter's  f ield-of-view  (measurement 
array  is  dimensioned  64x64) 


Print  the  filter's  template  and/or 
print  the  data  array 


Propagate  the  filter  estimates  forward 
in  time  by  one  sample  period 


Calculate  the  multiple  model  adaptive 
filter  estimate  prior  to  the  measurement 


Apply  control  by  redefining  location  of 
the  f ield-of-view  for  each  filter 
prior  to  the  next  measurement 


Test  for  loss-of-track ,  or  when  the 
estimate  of  target  position  is  within 
one  pixel  from  the  edge  of  the  FOV, 
and  perform  reaquisition  if  necessary 


Propagate  the  truth  model  states  forward 
in  time  by  one  sample  period 


Return  to  "Define  the  location  of 
the  f ield-of-view  for  each  filter 
with  respect  to  the  initial  position" 
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Figure  D-16.  CTR/SNG-S/T2/10-G/QB1 


Figure  E-8. 


Figure  E-9.  GM/ttlAF/T2/10-G/QBl 
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Figure  E-Z2.  GM/>fWF/T2/10-G/QB2 


Figure  E-28.  GM/WAF/TZ/10-G/QB2/R+4 


Figure  E-30.  GM/MMAF/T2/10-G/QB2/R+4 
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APPENDIX  G 


MULTIPLE  MODEL  ADAPTIVE  FILTER  SIMULATIONS 
TARGET  TRAJECTORY  SENSITIVITY  ANALYSIS 
(REFERENCED  PROM  SECTION  6.6) 
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